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PREFACE 


Upon  completion  of  my  graduate  work  at  the  Naval 
Postgraduate  School  I  will  be  assigned  as  an  Assistant 
Professor  in  the  Department  of  Mathematical  Sciences,  at  the 
United  States  Military  Academy.  While  there  I  will  teach  the 
cadets  calculus.  I  will  incorporate  this  package  on  integral 
calculus  that  I  have  developed  into  the  cadets'  instruction. 

As  an  instructor  I  will  see  this  integral  calculus  package 
used  by  the  cadets.  I  will  get  feedback  from  them  on  its  use 
and  its  value  in  the  instruction.  I  will  then  use  the 
feedback  to  make  improvements  to  the  package  and  to  expand  its 
capabilities.  It  is  my  vision  that  such  an  instructional  aid 
will  improve  the  overall  quality  and  effectiveness  of  the 
calculus  instruction  at  the  Military  Academy. 
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ABSTRACT 


One  of  the  most  difficult  concepts  in  the  area  of  integral 
multivariate  calculus  is  finding  the  limits  of  integration. 
This  thesis  describes  an  interactive  computer  program  designed 
to  help  students  understand  this  important  concept.  The 
program  shows  how  a  given  domain  is  plotted  and  teaches  how  to 
find  the  limits  of  integration  when  evaluating  two-dimensional 
integrals. 

The  program  allows  the  user  to  enter  any  known  information 
about  a  region  and  then  evaluates  the  integral.  The  region  is 
plotted,  the  limits  of  integration  are  given  along  with  the 
area  of  the  region.  The  program  handles  cartesian  and  polar 
coordinate,  two-dimensional  integral  problems.  This  program 
could  be  used  independently  by  the  student  and/or  used  by  the 
calculus  instructor  in  the  classroom. 
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THESIS  DISCLAIMER 


The  reader  is  cautioned  that  computer  programs  developed 
in  this  research  may  not  have  been  exercised  for  all  cases  of 
interest.  While  effort  has  been  made,  within  the  time 
available,  to  ensure  that  the  programs  are  free  of 
computational  and  logic  errors,  they  cannot  be  considered 
validated.  Any  application  of  these  programs  without 
additional  verification  is  at  the  risk  of  the  user. 
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I .  INTRODUCTION 


In  the  area  of  integral  multivariate  calculus  one  of  the 
most  difficult  topics  is  deciding  on  the  limits  of  integra¬ 
tion.  This  package  shows  the  calculus  student  how  a  given 
domain  is  plotted  and  teaches  the  student  how  to  find  the 
limits  of  integration  when  evaluating  two-dimensional 
integrals. 

Integral  calculus  is  the  mathematics  we  use,  for  example, 
to  find  lengths  and  areas  of  irregular  shapes;  to  calculate 
average  values  of  functions;  and  to  predict  future  values  in 
population  sizes  and  future  costs  of  living.  The  development 
of  integral  calculus  starts  from  the  calculation  of  areas. 
The  limits  used  to  define  areas  are  special  cases  of  a 
definite  integral. 

The  focus  of  the  package  is  on  teaching  the  student  how  to 
correctly  find  the  limits  of  integration  for  a  region,  while 
for  typical  calculus  problems  the  limits  of  integration  are 
determined  analytically,  it  is  quite  helpful  to  see  the 
regions  of  interest.  It  then  is  important  to  be  able  to 
accurately  graph  the  regions.  The  calculus  student  may  have 
difficulty  with  curve  sketching  especially  in  polar 
coordinates.  This  package  allows  the  student  to  visualize  the 
region  of  interest  enabling  the  student  to  better  understand 
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the  problem.  It  is  essential  that  the  student  be  able  to  see 
accurate  plots. 

This  package  was  developed  using  the  Mathematica  software. 
Mathematics  uses  an  adaptive  sampling  algorithm  for  plotting 
to  determine  when  and  where  to  sample  a  function.  A  section 
with  large  curvature  is  sampled  more  frequently  than  a  flat 
section  of  a  function.  This  sampling  algorithm  coupled  with 
a  set  prescribed  number  of  equally  spaced  samples  produces  the 
quite  reasonable  plots  we  want.  This  graphing  capability  is 
one  of  the  greatest  attributes  of  this  package. 

This  package,  which  teaches  the  student  how  to  find  the 
limits  of  integration,  has  two  parts.  The  first  part,  which 
is  found  in  Appendix  A,  consists  of  a  lab  book.  This  lab  book 
addresses  two-dimensional  integrals  and  walks  the  student 
through  numerous  example  problems.  The  two-dimensional  case 
considers  regions  given  in  the  cartesian  coordinate  system  as 
well  as  regions  in  the  polar  coordinate  system.  In  all 
examples  the  student  is  shown  the  plot  of  the  functions  and 
then  is  shown  how  the  limits  of  integration  are  determined. 
The  integrals  are  set  up  and  evaluated  and  the  area  of  the 
region  is  given. 

The  second  part  of  this  package  is  an  interactive  program 
written  in  Mathematica 'a  programming  language.  With  this 
program  the  student  can  solve  two-dimensional  integral 
problems,  given  in  both  the  cartesian  coordinate  system  and 
the  polar  coordinate  system.  The  student  enters  the  functions 
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which  define  a  region's  boundaries  and  any  other  known 
information  about  the  region.  The  curves  are  then  plotted. 
The  limits  of  integration  are  determined  and  the  area  of  the 
region  is  calculated.  The  program  output  includes  the  limits 
of  integration,  the  curves  and  the  calculated  area  of  the 
regions.  The  program  also  can  be  used  to  display  any  of  the 
plots  for  the  example  problems  found  in  the  lab  book. 

The  program's  greatest  value,  however,  is  in  allowing  the 
user  to  input  his  own  problems.  The  program  is  designed  to 
handle  those  types  of  integral  calculus  problems  found  in  any 
introductory  level  calculus  text  book. 

No  knowledge  of  Mathematics  is  required  to  use  the 
interactive  program.  The  tutorial,  found  in  Appendix  B,  tells 
the  user  how  to  load  the  program  and  how  to  use  it.  The 
tutorial  takes  the  user  through  three  example  problems.  The 
first  is  an  integral  problem  given  in  cartesian  coordinates. 
The  second  problem  is  in  polar  coordinates  and  finds  the  area 
between  a  polar  curve  and  the  origin.  The  final  problem  is 
also  in  polar  coordinates  and  finds  the  area  between  two  polar 
curves.  These  three  problems  are  simple  problems  that  show 
how  the  program  is  used.  The  tutorial,  while  not  compre¬ 
hensive,  gives  the  user  some  idea  of  the  types  of  problems  the 
program  handles.  It  demonstrates  how  the  program  prompts  the 
user  for  input  and  how  that  input  is  entered.  The  program  is 
capable  of  much  more  complicated  problems  than  those  in  the 
tutorial.  In  general,  once  the  package  is  loaded  and  the 
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program  is  started  the  user  simply  follows  the  instructions 
that  the  program  displays. 

Due  to  the  variety  of  ways  in  which  a  region  can  be 
defined  or  in  which  a  problem  can  be  worded  the  program  must 
be  versatile.  The  program  prompts  the  user  in  such  a  way  as 
to  Swlicit  the  necessary  information  about  the  region 
regardless  of  how  the  problem  is  stated.  As  mentioned 
previously  the  program  is  written  to  handle  the  types  of 
problems  found  in  calculus  text  books.  These  problems  are 
generally  well  defined  such  that  the  program  has  no  difficulty 
in  determining  the  region  of  interest.  On  the  other  hand  it 
is  quite  possible  to  contrive  problems  which  are  not  well 
defined  and  for  which  the  program  is  unable  to  determine  the 
region.  An  example  of  this  type  is  when  there  is  not  enough 
information  provided.  It  is  also  possible  to  have  a  problem 
which  is  seemingly  well  defined  for  which  the  program  does  not 
work.  An  example  of  this  type  of  problem  is  in  cartesian 
coordinates  when  the  functions  intersect  and  do  not  define  a 
single  unique  region.  For  instance  the  functions  y=12-x2, 
y=x  and  y=0  intersect  to  form  three  separate  regions  (see 
Figure  1,  page  5). 

If  such  a  problem  is  entered,  the  program  produces  either 
no  output,  bad  output,  or  an  error  message.  There  is  still  a 
way  in  which  the  program  can  be  used  to  solve  these  problems. 
For  problems  given  in  cartesian  coordinates  there  is  a 
cartesian  plot-only  option.  If  this  option  is  selected  the 
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user  can  input  all  the  functions  for  specified  ranges  of  the 
variables.  The  functions  will  then  be  displayed  and  the  plot 
can  be  used  to  estimate  the  limits  of  integration.  With  the 
estimated  limits  of  integration  obtained  from  the  plot  <?f  the 
functions  the  user  can  run  the  program  again  with  this  new 
information.  This  plot-only  option  can  be  used  also  if  the 
user  just  wants  to  see  what  a  particular  function  looks  like. 


y 


In  general  the  program  does  not  have  any  difficulty  with 
the  plots  of  the  functions  given  in  polar  coordinates.  If  the 
limits  of  integration  are  not  known  in  advance  the  program 
plots  the  polar  functions  from  Q=0  to  Q=2rr  radians  where  Q 
represents  the  polar  angle  theta.  For  the  two  options, 
finding  the  area  between  a  polar  curve  and  the  origin  and 
finding  the  area  between  two  polar  curves,  this  feature  acts 
as  a  plot-only  option.  The  values  of  Q  for  which  the  polar 
functions  intersect  are  displayed  when  finding  the  area 
between  two  curves.  The  values  of  Q  for  which  the  polar 
function  is  egual  to  zero  are  displayed  when  finding  the  area 
between  a  polar  curve  and  the  origin. 
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However  if  the  user  only  vunts  to  plot  the  polar  functions 
without  solving  the  integral  problem,  he  can,  both  for 
plotting  the  area  between  a  polar  curve  and  the  origin  and  for 
plotting  the  area  between  two  polar  curves. 
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II.  MOTIVATION 


Integral  calculus  has  a  wide  variety  of  applications  in 
mathematics,  physics,  engineering,  economics,  medicine,  and 
other  fields.  Integrals  are  used,  for  example,  to  compute 
areas,  volumes,  lengths  of  curves,  forces  of  fluids,  work, 
weight,  centers  of  gravity,  and  probabilities,  and  to 
determine  functions  from  their  rates  of  change.  In  these  many 
applications  often  the  most  difficult  part  is  in  properly 
setting  up  the  integrals  to  be  evaluated.  More  exactly  the 
difficulty  lies  in  determining  the  limits  of  integration  for 
these  integrals. 

Let  R  denote  the  region  in  the  xy-plane  that  is  bounded  on 
the  top  by  the  graph  y=f(x)  and  on  the  bottom  by  the  graph 
y=g(x)  and  that  extends  from  x=a  on  the  left  to  x=b  on  the 
right.  Double  integrals  over  R  may  be  expressed  as  iterated 
integrals  by  the  formula 

x=b  y=f(x)  -t 

JJ  h(x,y)dA  =  J  |  J  h ( x , y ) dy  Idx.  (1) 

R  x=a  L  y=g(x) 

To  evaluate  the  iterative  integral  we  start  at  the  innermost 
integral  and  work  toward  the  outside. 

The  most  difficulty  with  evaluating  two-dimensional 
integrals  given  in  cartesian  coordinates  is  in  determining 
g(x)  ,  f(x),  a,  and  b  which  are  limits  of  integration.  The 
region  R  may  have  either  changing  upper  or  lower  boundaries  in 
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which  case  the  region  R  is  composed  of  subregions.  The  double 
integral  over  R  then  becomes  equal  to  the  summation  of  double 
integrals,  that  are  over  each  of  the  subregions.  Then  for 
each  subregion  the  limits  of  integration  must  be  determined. 
This  situation  is  more  complicated,  making  determining  the 
limits  of  integration  difficult. 

The  focus  of  this  thesis  is  teaching  how  to  find  the 
limits.  As  such,  all  the  two-dimensional  integral  problems  in 
this  thesis  are  to  find  the  area  of  a  given  region.  That  is 
h(x,y)=l  and  so  we  have  for  integrals  given  in  cartesian 
coordinates 


Area= 


x=b  y=f (x) 

=  I  [  J  dy  1  dx* 

x=a  >-  y=g  (x)  J 


(2) 


After  evaluating  the  innermost  integral  the  equation  becomes 


x=b 

Area=  J  [ f (x) -g(x) ]dx. 
x=a 


(3) 


The  functions  y=f(x)  and  y=g(x)  can  be  found  by  graphing 
the  given  functions.  If  the  left  limit,  x=a,  and  the  right 
limit,  x=b,  are  not  known  they  are  determined  analytically. 
The  functions  y=f(x)  and  y=g(x)  are  equated  and  solved  for  x. 
Similarly  in  polar  coordinates  we  have 


rr  rQ=b  r  rr=r2(Q)  l 

Area=  Jj  rdA  =  j  If  rdr  dQ 

R  Q=a  L  r=rl(Q)  4 


(4) 
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where  both  rl  and  r2  are  functions  of  Q.  For  finding  the 
area  between  a  polar  curve  and  the  origin,  rl-0  and  we  have 


0-b  r=r2 (Q) 

Area*  J  [  J  rdr  J  dQ. 

Q*a  0 


If  we  evaluate  the  innermost  integral  we  get 


-Q=b  2 

Area  =  J  1  r2  dQ. 
Q=a  2 


(5) 


(6) 


For  finding  the  area  between  two  polar  curves  we  evaluate 
the  innermost  integral  in  equation  (4)  above  to  get 

Q=b  2  2 

Area=  J  1  [r2  -  rl  ]  dQ.  (7) 

Q=a  ^ 

Again  the  most  difficulty  with  two-dimensional  integrals 
given  in  polar  coordinates  is  setting  up  the  integrals.  More 
specifically  the  difficulty  is  in  determining  the  elements  a, 
b,  rl,  and  r2  which  are  the  limits  of  integration.  In  polar 
coordinate  integral  problems  the  graphs  of  the  polar  curves 
are  used  to  find  the  limits  of  integration.  In  finding  the 
area  between  a  polar  curve  and  the  origin  the  limits  a  and  b 
(if  not  known)  are  determined  by  solving  the  equation  r2=0  for 
r.  In  finding  the  area  between  two  polar  curves  the  limits  a 
and  b  (if  not  known)  are  determined  by  equating  rl  and  r2  and 
solving  for  r. 


The  calculus  student  must  do  numerous  example  problems  in 
order  to  become  proficient  and  confident  in  properly  setting 
up  an  integral.  The  more  problems  a  student  sees  and  does, 
the  more  confident  the  student  becomes.  This  package  provides 
numerous  worked  solutions  to  problems.  The  program  provides 
the  solutions  to  as  many  integral  problems  as  the  student  may 
wish  to  work.  The  program  can  be  used  with  whatever  calculus 
textbook  the  student  may  be  using. 

The  program  gives  solutions  that  include  not  only  the 
numerical  area,  but  also  the  graph  of  the  region  and  the 
limits  of  integration.  The  students  then  can  use  this 
information  to  identify  their  problem  areas.  The  students  may 
find  that  they  are  getting  the  right  limits  of  integration 
when  they  work  the  problems  but  that  their  solutions  are  not 
correct.  Then  it  is  likely  that  they  are  making  mistakes  when 
evaluating  the  integral.  In  this  case  the  students  see  that 
they  should  review  their  integration  techniques.  They  may  be 
finding  incorrect  limits  of  integration.  Then  the  graph  of 
the  region  may  be  helpful  in  determining  where  the  errors  are 
being  made.  At  any  rate  the  program  provides  a  great  deal  of 
useful  information.  Most  textbooks  provide  just  the  value  of 
the  integral. 

The  lab  book  contains  numerous  worked  out  example 
problems.  It  will  complement  the  few  worked  problems  that  are 
likely  to  be  found  in  most  textbooks.  The  fact  that  the 
package  is  adaptable  to  any  calculus  textbook  is  very  useful. 
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The  simple-to-use  program  gives  the  students  an  extremely 
efficient  and  effective  way  of  learning  how  to  find  limits  of 
integration.  The  use  of  mathematical  software  can  also  make 
learning  calculus  more  enjoyable.  Typically  students  are 
first  introduced  to  mathematical  software  in  a  higher  level 
course.  By  using  this  package  in  their  calculus  instruction 
the  students  are  introduced  to  mathematical  software  early 
during  their  education.  In  this  age  of  computers  this  becomes 
more  and  more  important  all  the  time. 

At  the  Naval  Postgraduate  School  there  is  currently  no 
introduction  to  mathematical  software  included  in  the  syllabus 
for  the  instruction  of  calculus.  While  it  is  true  that  the 
textbook  being  used  references  numerous  examples  and  applica¬ 
tions  of  different  mathematical  software,  typically  no  use  of 
mathematical  software  is  included  in  the  course  instruction. 

By  introducing  mathematical  software  in  the  instruction  of 
calculus  the  students  can  be  exposed  to  some  of  the  very 
powerful  uses  and  applications  of  mathematical  software.  With 
the  students'  introduction  to  mathematical  software  in  their 
calculus  instruction  via  this  package  a  major  step  is  being 
taken  in  changing  current  methods  of  instruction.  This 
package  could  be  used  at  the  United  States  Military  Academy  as 
well  as  other  colleges  and  universities. 
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III.  USB  OF  PACKAGE  IN  CALCULUS  INSTRUCTION 


This  section  addresses  incorporating  the  integral  package 
into  the  calculus  instruction  on  two-dimensional  integrals. 
The  integral  package  is  not  intended  to  be  used  by  itself  to 
teach  the  calculus  student  how  to  find  limits  of  integration 
and  calculate  two-dimensional  integrals.  Rather  it  should 
complement  the  classroom  instruction  and  the  homework  problems 
on  this  topic. 

The  package  should  be  introduced  to  the  students  after 
they  have  received  the  classroom  instruction  on  two- 
dimensional  integrals  given  in  cartesian  coordinates.  The 
students  should  also  first  attempt  to  work  a  number  of 
assigned  homework  problems  that  require  finding  the  area  of  a 
region  given  in  cartesian  coordinates.  After  the  students  see 
the  solutions  to  these  problems,  they  could  be  shown  the 
package . 

The  ideal  way  to  introduce  the  package  is  in  a  computer 
lab  in  which  all  students  can  have  hands-on  use  at  a  terminal. 
The  instructor  could  either  talk  the  students  through  turning 
on  Mathematica,  loading  the  package  and  using  it  or  let  the 
students  simply  follow  the  instructions  in  the  tutorial  found 
in  Appendix  B.  At  this  point  in  the  calculus  instruction  the 
students  would  only  be  ready  for  the  first  example  problem  in 
the  tutorial.  In  this  problem  the  students  would  use  the 
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program  to  find  the  area  of  a  region  given  in  cartesian 
coordinates. 

The  instructor  could  have  other  problems  for  the  students 
to  solve.  These  could  be  problems  given  in  cartesian 
coordinates  which  come  directly  out  of  the  calculus  textbook 
being  used.  The  greater  variety  of  problems  which  the 
students  solve  using  the  program  will  demonstrate  the 
program's  many  capabilities  to  solve  integral  problems  given 
in  cartesian  coordinates. 

Once  the  students  are  shown  how  to  use  the  integral 
program  part  of  the  package  they  should  be  introduced  to  the 
lab  book  which  contains  numerous  worked  problems.  These 
problems  show  the  students  step-by-step  how  to  find  the  limits 
of  integration  and  evaluate  the  integral.  These  example 
problems  will  supplement  the  examples  worked  in  the  classroom 
instruction  and  those  found  in  the  calculus  textbook  being 
used.  These  detailed  solutions  will  further  reinforce  the 
procedures  established  in  the  classroom  for  finding  the  limits 
of  integration  and  setting  up  the  proper  integrals. 

The  students  at  this  point  would  not  have  received  any 
classroom  instruction  on  integrals  involving  polar 
coordinates.  As  such  the  students  should  only  be  told  of  the 
program's  ability  to  handle  such  integral  problems.  The 
students  should  first  be  taught  how  to  graph  in  polar 
coordinates.  Once  the  students  have  received  this  classroom 
instruction  on  graphing  and  on  integrals  in  polar  coordinates 
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and  had  the  opportunity  to  work  some  of  these  problems  they 
should  be  shown  the  polar  coordinate  capabilities  of  the 
integral  program.  For  polar  coordinates  the  students  should 
be  shown  how  the  program  can  find  the  area  between  a  polar 
curve  and  the  origin  or  to  find  the  area  between  two  polar 
curves.  Again  the  tutorial  is  best  suited  for  this,  with  one 
example  of  each  of  these  problems. 

For  all  integral  problems  the  student  should  first  work 
the  problems  without  the  use  of  the  program.  The  program 
should  be  used  by  the  students  to  check  their  work  and  to 
identify  any  errors.  All  solutions  given  by  the  program 
include  the  plot  of  the  functions  and  the  region.  The 
students  can  then  see  graphically  presented  the  limits  of 
integration. 

It  is  important  that  the  student  always  first  attempt  to 
do  the  problems  on  their  own.  The  program  is  not  intended  to 
do  the  student's  homework  but  rather  it  is  to  assist  the 
student  in  learning.  The  goal  is  that  the  students  can,  by 
themselves,  find  the  limits  of  integration  regardless  how  a 
problem  is  presented.  The  program  is  an  educational  aid  and 
should  be  used  as  such.  This  integral  package  properly 
incorporated  into  the  calculus  instruction  can  make  learning 
to  find  limits  of  integration  for  two-dimensional  integrals 
easier  and  fun.  This  integral  package  complements  the 
classroom  instruction  on  finding  the  limits  of  integration  for 
two-dimensional  integrals. 
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IV.  WHY  KATHEKATICA 


While  there  are  many  types  of  mathematical  software  that 
I  could  have  used,  I  selected  Kathematiea.  Since  I  have  not 
used  any  other  mathematical  software  extensively  I  am  unable 
to  comment  on  the  relative  advantages  or  disadvantages  of  one 
software  over  another.  One  of  my  primary  reasons  for 
selecting  Kathematiea  is  its  availability  at  both  the  Naval 
Postgraduate  School  and  the  United  States  Military  Academy. 

Kathematiea  also  had  a  number  of  important  attributes 
which  lead  me  to  select  it.  One  such  attribute  of  the 
software  is  its  sophisticated  graphing  capabilities. 
Producing  quality  plots  was  essential  to  the  package. 
Kathematiea  enabled  me  to  use  a  variety  of  graphing 
capabilities  and  options.  Not  only  did  they  plot  the 
functions  but  also  provided  necessary  information  about  the 
regions.  This  information  was  used  to  find  limits  of 
integration,  determine  upper  and  lower  curves,  left  and  right 
curves,  or  inner  and  outer  curves,  and  finally  the  area  of  the 
region. 

A  second  attribute  of  the  software  was  its  programming 
capability.  Kathematiea  has  an  excellent  pattern  matching 
feature  which  was  used  extensively.  I  also  found  it  easy  to 
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program  and  debug  in  Mathematic*.  With  Professor  David 
Canright's  expertise  in  Mathematic*  which  was  quite 
accessible,  I  felt  confident  in  using  it. 

Another  advantage  or  rather  attribute  of  Mathematic*  is 
the  ability  for  students  to  use  the  package  and  other 
Mathematic*  capabilities  at  the  same  time.  In  using 
Mathematic*  I  was  also  able  to  take  advantage  of  many  of 
Mathematics'*  built-in  functions. 

Those  platforms  that  support  the  program  are  IBM 
compatible  personal  computers  with  the  Windows  environment, 
the  Unix  workstation  with  X-Windows,  and  Macintosh  computers. 
The  program  operation  described  herein  is  with  the  notebook 
interface  and  the  Windows  version  of  Mathematic*.  The  details 
differ  for  the  non-notebook  interface  found  with  the  Unix 
version  of  Mathematic*. 
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V.  STRUCTURE  OF  THE  PROGRAM 


The  integral  program  is  organized  into  ten  Mathematics 
packages  which  have  nine  major  subprograms  and  numerous  other 
defined  functions.  A  package  is  simply  a  file  containing 
Mathematics  definitions.  To  access  the  functions  defined  in 
the  packages,  the  packages  must  first  be  loaded.  The  file 
named  integral. ma  is  at  the  core  of  the  integral  program. 
This  file  connects  the  ten  files  together  along  with  the 
functions  defined  within  them.  When  integral. ma  is  loaded  it 
in  turns  loads  all  the  other  files  which  make  up  the  program. 

The  program  is  started  by  calling  the  subprogram  graph, 
found  in  file  integral. ma.  The  rest  of  the  program  is  built 
around  subprogram  graph  which  calls  up  the  other  subprograms. 
After  the  program  is  invoked  subprogram  graph  prompts  the  user 
to  respond  to  various  instructions.  It  is  in  this  manner  in 
which  different  options  are  selected.  Depending  on  which 
options  are  selected  the  appropriate  subprograms  are  called 
(see  Figure  2  on  next  page) . 
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Program  Architecture 


Subprogram 

graph 


To  display  a  graph  that 
is  in  the  lab  book  type 
in  the  example  number, 
otherwise  type  n. 


1  thru  23 


Type  1  for  cartesian 
coordinate  system, 
or  type  2  for  polar 
coordinate  system. 


To  just 
plot  curves 
type  y, 
otherwise 
type  n. 


Type  1  for  plot  of  area 
between  the  origin  and  a 
polar  curve,  or  type  2 
for  plot  of  area  between 
two  polar  curves. 


Type  1  for  the  area 
between  the  origin  and 
a  polar  curve,  or  type  2 
for  the  area  between 
two  polar  curves. 


If  the  user  desires  to  see  the  plot  of  the  region  for  any 
of  the  example  problems  in  Appendix  B,  subprogram  graph  calls 
subprogram  examples  (found  in  file  examples. ma)  .  If  the  user 
just  wants  to  plot  a  set  of  curves  given  in  cartesian 
coordinates  subprogram  graph  calls  up  subprogram  cartplot 
(found  in  file  cartplot. ma)  .  If  the  user  has  a  cartesian 
coordinate  problem  to  solve  subprogram  graph  calls  subprogram 
cart  (found  in  file  cart.ma) .  The  regions  in  these  types  of 
problems  must  be  given  as  functions  of  the  independent 
variable  x.  The  integration  is  then  with  respect  to  x.  For 
problems  given  in  cartesian  coordinates  subprogram  cart  calls 
up  subprogram  Inverse  (found  in  file  invar so. ma ) .  Subprogram 
inverse  solves  the  problem  by  integrating  over  the  region  with 
respect  to  the  variable  y.  Subprogram  inverse  first 
determines  if  integration  with  respect  to  the  variable  y  is 
appropriate.  If  the  functions  are  invertible  the  user  is 
then  asked  if  he  desires  to  integrate  the  same  problem  with 
respect  to  the  variable  y. 

If  the  user  wants  just  a  plot  of  the  area  between  the 
origin  and  a  polar  curve,  subprogram  graph  calls  subprogram 
polplotl  (found  in  file  polplotl.ma)  .  If  the  user  wants  just 
a  plot  of  the  area  between  two  polar  curves,  subprogram  graph 
calls  subprogram  polplot2  (found  in  file  polplot2 .ma) . 

If  the  problem  is  to  find  the  area  between  a  polar  curve  and 
the  origin,  subprogram  graph  calls  subprogram  polarl  (found 
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in  file  polarl.ma)  .  If  the  problem  is  to  find  the  area 
between  two  polar  curves,  subprogram  graph  calls  subprogram 
polar2  (found  in  file  polar2.ma). 

Many  of  the  example  problems  in  Appendix  B  have  regions 
which  are  shaded.  File  shades. ma  contains  the  Mathematics 
definitions  which  produce  the  different  shades  for  those 
example  problems. 
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▼I .  SUBPROGRAM  GRAPH 


Subprogram  graph  is  the  main  branch  of  the  program. 
Depending  upon  the  type  of  integral  problem  or  option  selected 
subprogram  graph  calls  the  appropriate  subprograms. 

After  the  program  is  started  the  first  prompt  says: 

To  display  a  graph  that  is  in  the  lab  book  type  in  the 
example  number,  otherwise  type  n. 

The  user  can  type  a  number  from  1  to  23  to  call  up  subprogram 
examples.  Subprogram  examples  displays  the  graph  which 
corresponds  to  that  number  example  problem  found  in  the  lab 
book. 

If  the  user  does  not  select  this  option  the  next  prompt 
appears : 

Type  l  for  cartesian  coordinate  system,  or  type  2  for 
polar  coordinate  system. 

If  the  user  types  1  for  cartesian  coordinate  system  the 
following  prompt  appears: 

To  just  plot  curves  type  y,  otherwise  type  n. 

If  the  user  types  y,  subprogram  graph  calls  up  subprogram 
cartplot  which  is  the  cartesian  plot-only  option.  Otherwise 
if  the  user  types  n,  subprogram  graph  calls  ups  subprogram 
cart  which  solves  those  problems  given  in  cartesian 
coordinates . 
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If  when  the  user  sees  the  prompt: 

Type  1  for  cartesian  coordinate  system,  or  type  2  for 
polar  coordinate  system. 

He  types  2  for  polar  coordinate  system,  the  following  prompt 
appears : 

To  just  plot  curves  type  y,  otherwise  type  n. 

If  the  user  types  y  to  select  the  polar  plot-only  option  he 
then  sees  the  next  prompt: 

Type  l  for  plot  of  area  between  the  origin  and  a  polar 
curve,  or  type  2  for  plot  of  area  between  two  polar 
curves . 

If  the  user  types  1,  subprogram  graph  calls  up  subprogram 
polplotl  which  allows  the  user  to  graph  the  area  between  a 
polar  curve  and  the  origin.  Otherwise  if  the  user  types  2, 
subprogram  graph  calls  up  subprogram  polplot2  which  allows  the 
user  to  graph  the  area  between  two  polar  curves. 

If  the  polar  plot-only  option  is  not  selected  the  user 
sees  the  following  prompt: 

Type  l  for  area  between  the  origin  and  a  polar  curve,  or 
type  2  for  area  between  two  polar  curves. 

If  the  user  types  l,  subprogram  graph  calls  up  subprogram 
polarl  which  finds  the  area  bounded  by  a  polar  curve  and  the 
origin.  Otherwise  if  the  user  types  2,  subprogram  graph  calls 
up  subprogram  polar2  which  finds  the  area  between  two  polar 
curves. 
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VII .  SUBPROGRAM  CART 


Subprogram  cart  solves  those  problems  given  in  cartesian 
coordinates.  Subprogram  graph  calls  up  subprogram  cart  when 
the  user  indicates  he  wants  to  solve  a  problem  given  in 
cartesian  coordinates.  Subprogram  cart  can  be  used  by  itself 
without  invoking  the  main  program.  In  order  to  use  this 
stand-alone  option  the  user  types  cart  and  then  activates  the 
cell.  The  tutorial  covers  the  many  ways  in  which  the  cell  can 
be  activated. 

Regardless  how  subprogram  cart  is  called,  the  subprogram 
attempts  to  solicit  all  known  information.  The  user  is  first 
told  to  enter  the  functions  defining  the  region.  These 
functions  must  be  functions  of  the  independent  variable  x. 
After  the  functions  have  been  entered  the  user  is  told  to 
indicate  if  any  of  the  limits  of  integration  is  known.  If  the 
user  indicates  yes  he  is  then  prompted  to  input  those  known 
limits.  The  subprogram  then  includes  the  functions  x=(left 
limit)  and  x= (right  limit)  as  appropriate,  to  those  functions 
of  x  already  entered.  If  the  left  limit  is  not  known  and  the 
right  limit  is  either  not  known  or  greater  than  zero,  the  user 
is  told  to  indicate  if  the  region  is  restricted  to  the  right 
half  plane.  If  the  right  limit  is  not  known  and  the  left 
limit  is  either  not  known  or  less  than  zero,  the  user  is  told 
to  indicate  if  the  region  is  restricted  to  the  left  half 
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plane.  In  some  problems  the  functions  will  intersect  and  form 
a  region  in  the  right  half  plane  and  one  in  the  left  half 
plane.  The  subprogram  must  know  which  region  to  consider.  An 
example  of  this  type  of  problem  is  the  intersection  of  the 
functions  y=12-x2,  y=x  and  y-0  (see  Figure  1  on  page  5)  . 

Once  all  this  information  is  obtained,  the  subprogram  is 
ready  to  solve  the  problem.  It  first  obtains  all  the  points 
of  intersection  for  each  of  the  curves  and  lines  segments. 
For  all  pairs  of  equations  it  does  this  by  simultaneously 
solving  two  equations  with  two  unknowns,  using  Mathematics' s 
NSolve  function.  The  transcendental  functions  defined  in  the 
program  are  Taylor  series  expansions.  When  equations  contain 
transcendental  functions  the  subprogram  uses  these  Taylor 
series  expansions  to  find  the  approximate  points  of 
intersection.  These  approximations  are  then  used  as  starting 
points  in  Mathematics' s  FindRoots  function  to  find  the  points 
of  intersection.  FindRoots  uses  Newton's  method  to  solve  for 
these  values. 

Subprogram  cart  uses  the  x  coordinates  of  the  points  of 
intersection,  for  each  function  entered  by  the  user,  to  define 
the  right  and  left  limits  for  the  curve.  If  these  x 
coordinates  are  the  same  value  the  subprogram  knows  that  the 
points  of  intersection,  for  the  entered  function,  form  a 
vertical  line  segment.  In  order  for  the  program  to  solve  a 
problem  each  curve  must  have  just  two  endpoints. 
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If  there  are  fewer  than  two  endpoints,  the  functions  do 
not  bound  a  closed  region.  If  there  are  more  than  two  points 
of  intersection,  the  functions  intersect  to  form  more  than  one 
region.  In  the  latter  case,  the  subprogram  is  unable  to  solve 
the  problem  correctly. 

Consider  the  example  where  the  lines  y*x  and  y*5x  are  the 
only  two  functions  entered.  These  functions  intersect  only  at 
the  point  whose  coordinates  are  (0,0).  The  program  then 
determines  only  one  endpoint  (i.e.,  the  point  (0,0))  for  each 
function.  These  two  functions  by  themselves,  do  not  bound  a 
closed  region  and  the  program  does  not  have  enough  informa¬ 
tion. 

Now  consider  the  example  where  the  functions  are  y=x,  y=0 
and  y=12-x2.  The  points  of  intersection,  and  thus  the 
endpoints  that  the  program  determines,  for  the  function  y=x 
are:  (-4,-4),  (0,0),  (3,3).  Those  for  the  function  y=0  are: 
(-3.464,0),  (0,0),  (3.464,0)  while  those  for  the  function 
y=12-x2  are:  (-4,-4),  (-3.464,0),  (3,3),  (3.464,0).  Each  of 
the  functions  has  more  than  two  endpoints  and  they  intersect 
to  form  three  separate  regions  (see  Figure  1  on  page  5) . 

If  however  the  user  indicates  that  the  region  is 
restricted  to  the  right  half  plane  the  program  discards  those 
points  of  intersection  which  are  in  the  left  half  plane. 
Similarly  if  the  user  indicates  that  the  region  is  restricted 
to  the  left  half  plane  the  program  discards  those  points  of 
intersection  which  are  in  the  right  half  plane. 
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Once  the  subprogram  knows  all  the  curves  and  their 
endpoints,  Mathematics's  Plot  function  samples  the  functions 
and  determines  the  points  which  define  each  curve.  Each  curve 
and  each  line  segment  is  then  represented  by  a  set  of  {x,y} 
points.  While  only  the  endpoints  are  needed  here,  the  list  of 
<x,y)  points  are  formed  to  be  used  later  in  subpi^gram 
inverse.  The  x  coordinates  of  the  points  of  intersection  for 
all  the  curves  are  taken  to  form  a  sorted  list  of  points. 
Each  adjacent  pair  of  points  in  this  list  are  the  limits  of 
integration  for  a  subregion  of  the  total  region  formed  by  the 
intersection  of  the  functions. 

Subprogram  cart  can  find  the  area  of  a  region  bounded  by 
two  curves  that  intersect  at  two  different  points.  In  this 
case  the  user  enters  the  two  functions  of  the  independent 
variable  x  and  then  indicates  that  no  limits  are  known.  An 
example  of  this  type  of  problem  is  the  region  bounded  by  the 
two  curves  y=x3  +  1  and  y=x2  +  x  (see  Figure  3  on  page  27)  . 
Subprogram  cart  can  find  the  area  of  the  region  bounded  by  two 
curves  that  intersect  at  one  point  and  by  a  line  segment 
perpendicular  to  the  x-axis.  In  this  case  the  user  enters  the 
two  functions  of  the  independent  variable  x,  indicates  that  a 
limit  is  known  and  then  enters  that  limit.  An  example  of  this 
type  of  problem  is  the  region  bounded  by  the  two  curves  y=x2 
and  y=0,  which  intersect  at  the  origin,  and  the  line  x=4  (see 
Figure  4  on  page  27) .  The  third  case  is  when  both  the  right 
and  left  limits  are  known.  Here  the  user  enters  the  two 
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functions  of  the  independent  variable  x,  indicates  that  the 
limits  are  known  and  then  inputs  both  limits.  An  example  of 
this  type  of  problem  is  when  the  two  curves  are  y=l  and  y=l-x  2 
with  the  given  limits  x=l  and  x=4  (see  Figure  5  on  page  28). 


y 


Figure  4 
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Figure  5 


The  maximum  number  of  functions  of  x  that  the  user  can 
input  is  three.  Then  the  maximum  number  of  subregions  which 
can  be  formed  is  two.  Hence  the  program  is  able  to  solve  for 
the  area  of  a  region  which  has  changing  boundaries.  For 
example  the  functions  y=5x,  y*x  and  y=cos  x  intersect  to  form 
a  region  which  consists  of  two  subregions  (see  Figure  6 
below) .  The  first  subregion  has  y=5x  as  its  upper  curve  and 
y=x  as  its  lower  curve.  The  second  subregion  has  y=cos  x  as 
its  upper  curve  and  y=x  as  its  lower  curve. 

y 

0  8 

0  6 

0  4 

0  2 

Figure  6 
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For  each  set  of  limit  points  and  corresponding  subregion, 
the  subprogram  determines  the  upper  and  lower  curves.  It  does 
this  by  selecting  those  curves  (now  represented  as  lists  of 
{x,y>  points)  whose  largest  and  smallest  x  coordinate  values 
bound  the  average  value  of  that  subregion's  two  limit  points. 
The  subprogram  then  uses  Mathematics 's  NIntegrate  function  to 
integrate  the  upper  curve  minus  the  lower  curve  over  the 
limits  of  that  subregion.  The  region  is  then  plotted.  The 
output  for  each  subregion  includes:  the  limits  of  integra¬ 
tion,  the  upper  and  lower  curves,  and  the  area.  The  area  of 
the  total  region  is  also  provided. 

In  order  to  solve  a  problem,  the  subprogram  cart  must  be 
able  to  identify  that  region.  It  was  mentioned  previously 
that  the  subprogram  has  limitations  if  the  functions 
intersect  to  form  more  than  one  region.  The  region  of 
interest  may  be  such  that  the  program  can  not  directly  solve 
the  problem.  An  example  of  this  type  of  problem  is  again  when 
the  functions  are  y=x,  y=0  and  y=12-x2  (see  Figure  1,  page  5)  . 
Suppose  that  the  region  of  interest  is  that  region  which  lies 
partly  in  the  left  half  plane  and  partly  in  the  right  half 
plane.  For  this  type  of  problem  the  user  would  have  to  plot 
the  functions  using  the  plot-only  option  (subprogram 
cartplot) .  When  the  user  indicates  that  his  problem  is  given 
in  cartesian  coordinates  he  is  given  the  option  to  just  plot 
the  functions.  If  this  option  is  selected  the  functions  are 
plotted  and  the  points  where  the  curves  intersect  are  provided 
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as  output.  From  the  plot  of  the  functions  and  the  points  of 
intersection  the  user  could  determine  the  limits  for  the 
region  of  interest.  The  user  could  then  use  this  information 
and  the  program  to  solve  the  problem. 

In  the  example  problem  above  the  plot  (see  Figure  1,  page 
5)  shows  that  the  upper  curve  for  the  subregion  in  the  left 
half  plane  is  y=12-x2  while  the  lower  curve  is  y=0.  For  the 
subregion  in  the  right  half  plane  the  upper  curve  is  y-12-x2 
and  the  lower  curve  is  y=x.  The  user  could  treat  each 
subregion  separately  and  use  the  program  to  solve  for  the 
areas  individually  and  then  add  the  two  areas. 

The  subprogram  has  no  difficulty  in  solving  problems  in 
which  the  functions  bound  a  single  region  which  lies,  either 
in  part  or  whole,  in  the  left  half  plane.  Thus  while 
subprogram  cart  has  limitations  the  plot-only  option  can  be 
used  in  conjunction  with  it  to  solve  certain  of  these  types  of 
problems. 
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VIII.  SUBPROGRAM  INVERSE 


Subprogram  inverse  integrates  those  problems  given  in 
cartesian  coordinates  with  respect  to  the  variable  y.  First, 
however,  the  subprogram  determines  if  it  is  appropriate  to 
integrate  a  problem  with  respect  to  y.  It  attempts  to  find 
the  inverse  of  each  function  of  x  that  was  input.  If  not  all 
the  inverse  functions  exist  subprogram  inverse  stops  execu¬ 
tion.  If  on  the  other  hand  they  all  exist  the  user  is  asked 
if  he  wants  to  integrate  the  same  problem  with  respect  to  y. 
Subprogram  inverse  uses  Mathematics's  XnverseFunction  to  find 
the  inverse  functions. 

Subprogram  inverse  works  similar  to  subprogram  cart  to 
solve  the  problem.  The  y-values  of  the  points  of  intersection 
for  the  curves  (found  in  subprogram  cart),  are  limits  of 
integration.  There  can  exist  limits,  however,  that  are  not 
any  of  the  y-values  of  the  points  of  intersection.  An  example 
of  this  type  of  problem  is  when  the  functions  are  y=12-x2  and 
y=-3+x+x2  (see  Figure  7  on  page  32).  These  functions 
intersect  at  the  points  (-3,3)  and  (2.5,5.75).  The  first 
function  is  the  upper  curve  and  takes  on  a  maximum  y-value  of 
12.  The  second  function  is  the  lower  curve  and  takes  on  a 
minimum  y-value  of  -3.25.  Hence  there  are  three  subregions. 
The  limits  of  integration  for  the  first  region  are  y=-3.25  and 
y=3.  The  limits  of  integration  for  the  second  region  are  y=3 
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and  y=5.75.  The  limits  of  integration  for  the  third  region 
are  y»5.75  and  y=12.  The  y-values  -3.25  and  12  are  limits 
while  they  are  not  any  of  the  y-values  for  the  points  of 
intersection. 

y 


Figure  7 


The  subprogram  finds  the  minimum  and  maximum  values  of  the 
functions  and  includes  them  as  limit  points.  In  subprogram 
cart  all  the  curves  are  represented  by  lists  of  points.  From 
these  lists  subprogram  inverse  selects  the  minimum  and  maximum 
values  which  y  takes  on  in  the  curves.  Once  all  the  limits  of 
integration  are  known  they  are  placed  in  a  sorted  list.  Each 
adjacent  pair  of  points  are  then  limits  of  integration  of  a 
subregion  of  the  total  region. 

For  each  subregion  the  subprogram  determines  the  right  and 
left  curves  in  the  same  manner  in  which  subprogram  cart  finds 
the  upper  and  lower  curves.  Some  inverse  functions  have  a 
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positive  branch  and  a  negative  branch.  These  functions  are 
identified  as  having  a  square  root  term  involving  the  variable 
x.  Subprogram  inverse  checks  to  see  if  either  the  right  curve 
or  left  curve  for  a  subregion  has  this  type  of  inverse 
function.  If  an  inverse  function  is  of  this  type,  the 
subprogram  must  determine  whether  it  is  the  positive  or  the 
negative  branch. 

If  it  is  the  negative  branch,  the  subprogram  calls 
routine  negbranch,  which  transforms  the  inverse  function.  The 
term  in  the  inverse  function  involving  the  square  root  is  then 
replaced  with  the  negative  of  that  term.  For  example  consider 
the  region  bounded  by  the  functions  y=0,  y^-4+x2  and  x=-4  (see 
Figure  8  on  page  34)  .  The  inverse  of  y=-4+x2  is  x=/  y+4 . 
When  the  region  is  integrated  with  respect  to  y  the  right 
curve  is  x-J  y+4  and  the  left  curve  is  x=-4.  But  the  right 
curve  is  a  negative  branch  so  we  have  x=-y  y+4  as  the  right 
curve.  In  the  calculation  of  the  integral  the  integrand  is 
the  right  curve  minus  the  left  curve,  or  (-7  y+4)  -  (-4)  = 
-J  y+4+4 . 

If  the  right  curve  for  a  subregion  is  the  same  as  the 
left  curve,  the  subprogram  calls  another  routine  bothbranches . 
This  routine  works  similar  to  the  routine  negbranch.  It 
insures  that  the  region  integrated  over  is  the  positive  branch 
minus  the  negative  branch  (or  twice  the  positive  branch) .  For 
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example  consider  the  region  bounded  by  the  functions  y=x2  and 
y=4  (see  Figure  9  on  page  35) .  The  inverse  function  of  y=x2 
is  x=y“y.  But  in  this  problem  the  right  curve  is  the  positive 
branch  and  the  left  curve  is  the  negative  branch.  Therefore 
the  integrand  is  -  (~y~x)  =  2  J~lc,  which  is  the  right 

curve  minus  the  left  curve. 

The  area  of  the  subregions  are  then  found  using 
Mathematics 'a  NIntagrata  function.  The  output  for  each 
subregion  includes:  the  limits  of  integration,  the  left  and 
right  curves,  and  the  area.  The  area  of  the  total  region  is 
also  provided. 

Subprogram  invarsa  as  mentioned  previously  works  only  when 
all  inverse  functions  exist.  Also  stated  was  that  it  only 
integrates  with  respect  to  the  variable  y  when  it  is 
appropriate  to  do  so.  There  are  functions  whose  inverse 
exists  but  which  are  complicated  functions  to  integrate 
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Figure  9 

(Mathematic®  has  no  problem  but  the  calculus  student  might) . 
For  those  problems  involving  transcendental  functions  the 
subprogram  will  not  integrate  with  respect  to  y,  even  if  such 
integration  is  possible.  Integration  with  respect  to  the 
variable  y  is  supposed  to  save  the  calculus  student  time. 
Usually  the  inverse  of  a  function  containing  transcendental 
functions  is  complicated  and  more  difficult  to  integrate. 


IX.  SUBPROGRAM  CARTPLOT 


Subprogram  eartplot  is  an  option  which  allows  the  user  to 
graph  functions  given  in  cartesian  coordinates.  After  the 
user  indicates  that  a  problem  is  given  in  cartesian 
coordinates,  subprogram  graph  asks  if  the  user  just  wants  to 
plot  the  functions.  If  the  user  selects  this  option  sub¬ 
program  graph  calls  up  subprogram  eartplot.  The  user  first 
enters  all  functions  of  the  independent  variable  x  and  then 
those  lines  perpendicular  to  the  x-axis.  The  user  enters  the 
ranges  for  both  the  independent  variable  x,  and  the  dependent 
variable  y.  The  functions  are  then  plotted  and  the  points 
where  the  curves  intersect  are  given. 

This  plot-only  option  can  be  used  to  estimate  the  coordin¬ 
ates  of  where  the  functions  intersect  and  thus  the  limits  of 
integration.  There  is  no  limit  on  the  number  of  functions 
which  can  be  plotted  on  one  graph  with  this  option.  A  single 
function  can  be  plotted  if  desired.  Subprogram  eartplot  uses 
Mathematics' a  Plot  function  to  plot  all  graphs. 

Subprogram  eartplot  can  be  used  by  itself  without  invoking 
the  main  program.  This  is  valuable  in  saving  time  if  the  user 
just  wants  to  plot  several  graphs.  In  order  to  use  this 
stand-alone  option  the  user  types  eartplot  and  then  activates 
the  cell.  The  tutorial  covers  the  many  ways  in  which  the  cell 
can  be  activated. 
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Regardless  how  subprogram  oartplot  is  called  the  first 
prompt  the  user  sees  is: 

To  just  plot  curves  type  y,  otherwise  type  n. 

Suppose  the  user  wants  to  plot  the  functions  y=x,  y=o  and 
y*12-x2.  The  user  types  y  to  select  the  plot-only  option. 
The  next  message  the  user  sees  is: 

Enter  the  curves ,  one  at  a  time.  When  finished  type  n. 

y  = 

The  user  then  types  x  and  presses  Enter  to  input  the  first 
function.  The  above  message  appears  again  and  the  user  then 
types  0  and  presses  Enter  to  input  the  second  function.  Once 
again  the  above  message  appears  and  the  user  inputs  the  last 
function  by  typing  12-z2  and  pressing  Enter.  The  above 
message  appears  again  and  this  time  the  user  types  n  and 
presses  Enter  to  indicate  that  he  is  finished. 

The  user  then  sees  the  message: 

Enter  lines  perpendicular  to  the  z-axis.  When  finished 
type  n. 

Since  there  are  none  the  user  types  n.  If  there  were  lines  to 
enter  they  would  have  been  entered  one  at  a  time  as  the 
functions  of  the  independent  variable  x  were  entered  above. 
The  user  then  sees  the  message: 

Enter  lover  limit  for  z-range  of  plot. 

The  user  must  guess  what  this  value  is  as  well  as  the  upper 
limit  for  the  x-range  and  both  the  lower  and  upper  limits  for 
the  y-range.  Suppose  the  user  knows  these  values.  For 
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example,  to  enter  the  lower  limit  for  the  x-range  he  then 
types  -5  and  presses  Enter.  The  user  then  sees  the  message: 

Enter  upper  limit  for  x-range  of  plot. 

He  types  5  and  presses  Enter  and  sees  the  message: 

Enter  lover  limit  for  y- range  of  plot. 

He  types  -5  and  presses  Enter  and  sees  the  message: 

Enter  upper  limit  for  y-range  of  plot. 

He  types  12  and  the  plot  is  displayed  along  with  the  points 
where  the  curves  intersect  (see  Figure  10  below).  After  the 
plot  is  displayed  the  user  may  want  to  change  the  x  or  y  range 
of  the  plot.  To  do  this  the  user  must  call  up  subprogram 
eartplot  again  and  enter  all  the  information. 

y 


X 
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The  points  where  the  curves  intersect  are: 

{{-4., -4.),  (-3.464,0.),  { 0 . , 0 . ) ,  (3. ,3.),  (3.464,0)) 

Figure  10 
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X.  SUBPROGRAM  POLAR1 


Subprogram  polarl  finds  the  area  bounded  by  a  polar  curve 
and  the  origin.  Subprogram  graph  calls  up  subprogram  polarl 
when  the  user  indicates  the  problem  is  to  find  the  area 
between  the  origin  and  a  polar  curve.  Subprogram  polarl  can 
be  used  by  itself  without  invoking  the  main  program.  In  order 
to  use  this  stand-alone  option  the  user  types  polarl  and  then 
activates  the  cell. 

Regardless  how  subprogram  polar  is  called,  the  user  first 
enters  the  polar  curve  r,  which  bounds  the  region.  The  polar 
curve  must  be  a  function  of  the  independent  variable  Q 
representing  the  polar  angle  theta.  The  subprogram  asks  the 
user  if  the  limits  of  integration  are  known  and  then  to  enter 
them,  if  they  are  known.  For  those  problems  in  which  the 
limits  are  known  in  advance  the  subprogram  then  graphs  the 
region  and  solves  the  problem.  The  subprogram  uses  Mathe¬ 
matics'  s  PolarPlot  function  to  produce  all  plots.  The  output 
includes  the  limits  of  integration,  the  polar  curve  and  the 
area  bounded  by  the  curve. 

A  polar  curve  may  be  negative  for  certain  ranges  of  the 
independent  variable  Q.  For  example  the  polar  curve  r=sin3Q 
takes  on  negative  values  in  the  range  Q=*/3  to  Q=2r/3.  A 
polar  curve  may  also  be  undefined  for  certain  ranges  of  the 
independent  variable  Q  (i.e.,  when  equal  to  the  square  root  of 


39 


a  negative  number) .  The  polar  curve  r=y  cos  2Q  is  defined 
only  for  Q  between  Q=-ir/4  and  Q=r/4  and  for  Q  between  Q=3»/4 
and  Q=5r/4  (see  Figure  11  below) . 
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Figure  11 

The  subprogram  determines  for  which  values  of  Q  the 
function  r  is  equal  to  zero.  It  uses  Mathematics' a  MSolve 
function  to  do  this.  For  all  regions  between  these  values  of 
Q  (where  r=0)  and  the  values  of  Q  entered  as  limits,  the 
function  r  must  be  positive,  negative,  or  undefined.  The 
subprogram  determines  the  case  for  each  such  region. 

When  the  function  r  is  defined  for  a  region  the  subprogram 
uses  Mathematics 'a  MIntegrate  function  to  find  the  area  of 
that- region.  The  limits  of  integration  are  those  values  of  Q 
which  define  the  region.  Since  the  problem  is  to  find  the 
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area  between  the  polar  curve  and  the  origin  the  absolute  value 
of  the  evaluated  integral  is  taken.  Therefore  a  region  where 
the  polar  curve  r  is  negative  contributes  positive  area  to  the 
total  area  between  the  curve  and  the  origin.  If  the  function 
r  is  not  defined  over  a  region,  no  area  is  contributed. 

If  the  user  does  not  know  the  limits  of  integration  the 
polar  curve  is  plotted  for  Q=0  to  Q*2r.  The  plot  of  the  curve 
is  displayed  along  with  those  Q  values  for  which  the  function 
r  is  equal  to  zero.  After  a  brief  delay  the  user  is  told  that 
the  polar  curve  is  plotted  from  Q-0  to  Q=2ir  and  is  asked  to 
input  the  limits  of  integration.  From  this  point  on  the 
subprogram  works  exactly  like  it  did  above  when  the  limits  of 
integration  were  known  in  advance. 

If  the  area  bounded  by  a  polar  curve  and  the  origin  is 
made  up  of  more  than  one  region,  the  limits  of  integration  for 
all  those  regions  are  provided  as  output.  For  example  the 
area  between  the  polar  curve  r=J  cos  2Q  and  the  origin  for  Q=0 
to  Q=3jt/2  consists  of  two  separate  regions  (see  Figure  12  on 
page  42) .  The  limits  of  integration  for  the  first  region  are 
al=0.  and  bl=ir/4  while  the  limits  of  integration  for  the 
second  region  are  a2=3ir/4  and  b2=5>r/4.  These  limits  are 
displayed  along  with  the  graph  of  the  region  and  the  total 
area. 
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The  subprogram  does  not  tell  how  much  area  is  contributed 
by  each  region.  If  the  user  wants  to  know  this,  he  could 
treat  each  region  as  a  separate  problem.  The  program  could 
then  be  used  to  determine  the  areas  individually,  one  at  a 

time. 
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ZZ .  SUBPROGRAM  POLPLOT1 


Subprogram  polplotl  is  an  option  which  allows  the  user  to 
graph  the  area  between  the  origin  and  a  polar  curve.  After 
the  user  indicates  that  a  problem  is  given  in  polar 
coordinates,  subprogram  graph  asks  if  the  user  just  wants  to 
plot  the  functions.  If  the  user  selects  this  plot-only  option 
the  program  then  asks  the  user  to  select  the  type  of  plot. 
The  choices  are:  (1)  the  plot  of  the  area  between  the  origin 
and  a  polar  curve,  and  (2)  the  plot  of  the  area  between  two 
polar  curves.  When  the  user  selects  (1),  subprogram  graph 
calls  up  subprogram  polplotl. 

The  user  enters  the  polar  curve  which  bounds  the  region 
and  then  indicates  if  the  limits  are  known.  If  the  limits  are 
not  known,  the  polar  curve  is  plotted  for  the  polar  angle  Q=0 
to  Q=2jt.  Those  points  where  the  polar  curve  is  equal  to  zero 
are  displayed  along  with  the  graph.  The  user  is  then  given  an 
option  to  change  the  limits  T*  the  user  selects  this  option 
the  polar  curve  is  then  replotted  for  the  new  limits  of  the 
polar  angle  Q,  which  the  user  enters.  Subprogram  polplotl 
uses  Mathematics 's  PolarPlot  function  to  plot  the  graphs. 

Subprogram  polplotl  can  be  used  by  itself  without  invoking 
the  main  program.  In  order  to  use  this  stand-alone  option  the 
user  types  polplotl  and  then  activates  the  cell. 
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Regardless  how  subprogram  polplotl  is  called  the  first 

prompt  the  user  sees  is: 

Enter  the  polar  curve  which  bounds  the  region. 

If  the  user  wanted  to  plot  the  polar  curve  r=  1  cos  Q,  he 
would  type  1  -  cos[Q]  and  press  Enter.  The  next  message  the 
user  would  see  is: 

If  the  limits  ere  known  type  y,  otherwise  type  n. 

Suppose  the  user  types  n.  The  polar  curve  is  then  displayed 
along  with  those  points  where  the  polar  curve  is  equal  to  zero 
(see  Figure  13  below) . 
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The  polar  curve  r  is  equal  to  zero  when  Q  is:  (0,  2  Pi) 


Figure  13 
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After  a  brief  delay  the  user  sees  the  message: 

The  polar  curve  is  plotted  from  Q=o  to  Q=2  pi.  if  you 
want  to  change  the  limits  type  y,  otherwise  type  n. 


Now  suppose  the  user  wants  to  change  the  lower  limit  to  Q=x/2 
and  the  upper  limit  to  Q=3v/2.  He  types  y  and  the  next 
message  that  appears  is: 

Input  new  lower  limit. 

The  user  types  pi/2  and  presses  Enter.  The  next  message 
that  appears  is: 

Input  new  upper  limit. 

The  user  types  3  pi/2  and  presses  Enter.  The  plot  of  the 
polar  curve  is  then  displayed  for  the  new  range  of  the  polar 
angle  Q  (see  Figure  14  below) . 
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XI I.  SUBPROGRAM  POLAR2 


Subprogram  polar2  finds  the  area  between  two  polar  curves. 
Subprogram  graph  calls  up  subprogram  polar2  when  the  user 
indicates  the  problem  is  to  find  the  area  bounded  between  two 
polar  curves.  Subprogram  polar2  can  be  used  by  itself  without 
invoking  the  main  program.  In  order  to  use  this  stand-alone 
option  the  user  types  polar2  and  then  activates  the  cell. 

Regardless  how  subprogram  polar2  is  called,  the  subprogram 
directs  the  user  to  enter  the  two  polar  curves.  These  polar 
curves  must  be  functions  of  the  independent  variable  Q.  After 
the  two  polar  curves  have  been  entered  the  user  is  asked  if 
the  limits  of  integration  are  known  and  to  input  them  if 
known . 

The  program  determines  those  values  of  Q  where  the  curves 
intersect.  The  functions  which  define  the  curves  are 
represented  as  Taylor  Series  expansions.  The  subprogram  uses 
Mathematics 's  NRoots  function  to  obtain  the  approximate  values 
of  Q  for  which  the  curves  intersect.  These  values  are  then 
used  as  starting  points  in  Mathematics 's  FindRoot  function. 
This  function  uses  Newton's  Method  to  find  the  points  of 
intersection. 

If  the  limits  are  not  known  the  subprogram  then  plots  the 
two  polar  curves  for  Q=0  to  Q=2*.  All  plots  are  produced 
using  Mathematics's  PolarPlot  function.  The  values  of  Q  for 
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which  the  two  polar  curves  intersect  are  also  provided. 
After  a  brief  delay  the  user  is  told  that  the  curves  are 
plotted  for  Q=0  to  Q-2r .  The  subprogram  asks  the  user  if  he 
wants  to  change  these  limits.  If  the  user  desires  he  can  then 
input  the  new  limits.  If  the  limits  of  integration  are  known 
and  inputted  iritially,  the  program  plots  the  two  polar  curves 
over  this  region. 

Once  the  limits  are  entered  the  subprogram  calls  up  a 
number  of  routines  which  check  for  certain  types  of 
conditions.  Routine  negcheck  determines  if  the  polar  curves 
are  negative.  The  user  is  told  if  one  or  more  of  the  curves 
are  negative  in  the  range  Q=lower  limit  to  Q=upper  limit.  The 
user  is  also  told  that  the  program  does  not  handle  this  type 
of  problem.  Routine  complexcheek  determines  the  ranges  of  Q 
for  which  the  curves  are  undefined  (if  they  are).  Routine 
diffcheck  checks  to  see  if  the  outer  and  inner  polar  curve 
changes  between  the  upper  and  lower  limits.  If  the  outer 
curve  and  inner  curve  changes  the  subprogram  lets  the  user 
know  this.  The  user  is  told  that  a  point  of  intersection 
cannot  lie  between  the  two  limits.  The  user  is  then  directed 
to  input  new  limits. 

Once  all  conditions  have  been  checked  and  the  program 
determines  that  the  problem  is  well  defined  it  calculates  the 
area  between  the  two  polar  curves.  The  subprogram  has  two 
routines  which  it  uses  to  find  the  area.  If  both  curves  are 
defined  over  the  entire  interval  routine  findarea2  is  used. 
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This  routine  uses  Mathenati ca'a  NIntegrate  function  to 
evaluate  the  single  integral.  If  htwever  the  inner  curve  is 
undefined  over  certain  ranges  of  Q  the  routine  findareal  is 
used  to  find  the  area. 

An  example  is  when  the  outer  curve  is  the  circle  r=l  and 
the  inner  curve  is  r=J  cos  2Q  (see  Figure  15  below) .  Here  the 
inner  curve  is  not  defined  for  the  polar  angle  Q=ir/4  to  Q=3ir/4 
and  for  Q=-3ir/4  to  -ir/4.  Routine  findareal  first  finds  the 
area  between  the  outer  polar  curve  and  the  origin  and  then 
finds  the  area  between  the  inner  polar  curve  and  the  origin. 
The  difference  between  these  two  areas  is  the  area  between  the 
two  polar  curves. 


y 


Figure  15 
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The  region  bounded  by  the  two  polar  curves  is  then 
plotted.  The  output  includes  the  limits  of  integration,  the 
outer  curve,  the  inner  curve,  and  the  total  area  of  the 
region.  The  program  does  not  provide  the  areas  of  any 
subregions.  If  the  user  wanted  to  know  this,  each  subregion 
could  be  treated  as  a  separate  problem.  The  program  could 
then  be  used  to  find  the  area  of  each  individual  subregion, 
one  at  a  time. 
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XIII.  SUBPROGRAM  POLPLOT2 


Subprogram  polplot2  is  an  option  which  allows  the  user  to 
graph  the  area  between  two  polar  curves.  After  the  user 
indicates  that  a  problem  is  given  in  polar  coordinates, 
subprogram  graph  asks  if  the  user  just  wants  to  plot  the 
functions.  If  the  user  selects  this  plot-only  option  the 
program  then  asks  the  user  to  select  the  type  of  plot.  The 
choices  are:  (1)  the  plot  of  the  area  between  the  origin  and 
a  polar  curve,  and  (2)  the  plot  of  the  area  between  two  polar 
curves.  When  the  user  selects  (2),  subprogram  graph  calls  up 
subprogram  polplot2. 

The  user  enters  the  two  polar  curves  which  bound  the 
region  and  then  indicates  if  the  limits  are  known.  If  the 
limits  are  not  known,  the  polar  curves  are  plotted  for  the 
polar  angle  Q=0  to  Q=2r.  Those  points  where  the  two  polar 
curves  intersect  are  displayed  along  with  the  graph.  The  user 
is  then  given  an  option  to  change  the  limits.  If  the  user 
selects  this  option  the  polar  curves  are  then  replotted  for 
the  new  limits  of  the  polar  angle  Q,  which  the  user  enters. 
Subprogram  polplot2  uses  Mathematics' a  PolarPlot  function  to 
plot  the  graphs. 

Subprogram  polplot2  can  be  used  by  itself  without  invoking 
the  main  program.  In  order  to  use  this  stand-alone  option  the 
user  types  polplot2  and  then  activates  the  cell. 
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Regardless  how  subprogram  polplot2  is  called  the  first 
prompt  the  user  sees  is: 

Enter  one  of  the  polar  curves  which  bounds  the  region. 
Suppose  the  user  wanted  to  plot  the  area  between  the  two  polar 
curves  r=  2  and  r=  1  +  sin  Q.  To  enter  the  first  polar  curve 
the  user  would  type  2  and  press  Enter.  The  next  message  the 
user  would  see  is: 

Enter  the  other  polar  curves  which  bounds  the  region. 

The  user  would  type  1  4-  sin[Q]  and  press  Enter  to  enter  the 
second  polar  curve.  The  next  message  the  user  would  see  is: 

If  the  limits  are  known  type  y,  otherwise  type  n. 

Suppose  the  user  types  n.  The  polar  curves  are  then  displayed 
along  with  those  points  where  the  polar  curves  intersect  (see 
Figure  16  on  page  52). 

After  a  brief  delay  the  user  sees  the  message: 

The  polar  curves  are  plotted  from  Q=C  to  Q=2  pi.  If  you 
want  to  change  the  limits  type  y,  otherwise  type  n. 

Now  suppose  the  user  wants  to  change  the  upper  limit  to 
Q=v.  He  types  y  and  the  next  message  that  appears  is: 

Input  new  lower  limit. 

The  user  types  0  and  presses  Enter.  The  next  message  that 
appears  is: 

Input  new  upper  limit. 

The  user  types  pi  and  presses  Enter.  The  plot  of  the  two 
polar  curves  is  then  displayed  for  the  new  range  of  the 
polar  angle  Q  (see  Figure  17  on  page  52) . 
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The  two  polar  curves  intersect  at  Q 


Figure  16 

y 


Figure  17 


=  Pi 
~2 
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XIV.  SUBPROGRAM  EXAMPLES 


The  file  named  examples. ma  contains  the  Mathematics  code 
for  each  of  the  graphs  found  in  the  lab  book.  All  of  the 
graphs  are  in  a  package  and  can  be  displayed  by  the  subprogram 
examples.  When  in  the  main  program  the  first  prompt  says: 

To  display  a  graph  that  is  in  the  lab  book  type  the 
example  number,  otherwise  type  n. 

The  user  can  then  type  a  number  from  1  to  23  to  call  up 
subprogram  examples.  The  subprogram  then  displays  the  graph 
which  corresponds  to  that  number  example  problem. 

There  may  be  only  limited  value  in  displaying  two- 
dimensional  graphs  which  are  already  in  the  lab  book. 
However,  when  this  integral  package  is  expanded  to  include  the 
three-dimensional  case  this  could  be  a  very  useful  option. 
Mathematics  has  a  sophisticated  three-dimensional  graphing 
capability  which  uses  a  three-color  lighting  system  to 
illuminate  three  dimensional  surfaces.  Mathematica  also  has 
the  ability  to  change  the  point  of  view  interactively  to 
explore  three-dimensional  graphics.  For  these  types  of 
problems  being  able  to  display  the  graphs  becomes  important. 
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XV.  FUNCTIONS  AND  NOTATION 


A.  FUNCTIONS 

When  inputting  information  both  the  program  and 
Mathematic*  are  case  sensitive.  Some  transcendental  functions 
are  defined  within  the  program  code  and  must  be  in  lower  case 
letters.  These  transcendental  functions  are  represented  as 
Taylor  Series  expansions.  They  were  needed  to  solve  algebrai¬ 
cally  for  where  two  curves  intersect,  using  NRoots.  These 
approximate  points  were  then  used  as  starting  points  for 
Mathematics's  FindRoot  function  which  uses  Newton's  Method. 

Mathematic*  has  its  own  built-in  transcendental  functions 
which  in  most  instances  are  spelled  exactly  the  same.  The 
only  difference  is  that  all  Mathematic*  functions  begin  with 
an  upper  case  letter. 

Below  are  those  functions  which  the  program  recognizes: 


sin[x] 

sec[x] 

cos[x] 

exp[x] 

tan[x] 

ln[x] 

cot[x] 

sqrt [x] 

csc[x] 

B.  NOTATION 

1.  Brackets  and  Parentheses 

Square  brackets  and  parentheses  are  intended  for 
different  purposes.  Square  brackets  are  used  for  specifying 


arguments  of  functions.  Parentheses  are  used  for  grouping. 
Without  parentheses,  multiplication  and  division  have  a  higher 
precedence  than  addition  and  subtraction. 

2.  Mathematical  symbols 

The  standard  mathematical  operations  are  referred  to 
with  symbols.  Those  symbols  used  are: 


Mathsmat  iggO-Saabal  operation 


+  plus,  add 

minus,  subtract 
*  times,  multiply 

/  divide 

A  power 


3.  other  Symbols 

The  symbol  Q  is  used  for  theta  in  problems  given  in 
polar  coordinates.  All  trigonometric  functions  must  be 
functions  of  the  variable  Q.  The  program  does  not  recognize 
any  other  such  variable.  The  symbol  pi  represents  the 
numerical  value  of  the  mathematical  constant  * . 
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XVI .  FUTURE  WORE 


Right  now  the  program  just  handles  two-dimensional 
problems.  The  package  could  be  extended  to  include  three- 
dimensional  problems.  The  three-dimensional  case  would 
include  problems  given  in  cartesian  as  well  as  cylindrical  and 
spherical  coordinate.  Three-dimensional  problems  could  also 
be  included  in  the  lab  book.  Subprogram  examples  then  could 
display  the  graphs  for  these  example  problems.  Mathematics 's 
representation  of  three-dimensional  plots  would  give  the 
student  a  feel  for  the  depth  of  the  region  and  insight  into 
how  the  surfaces  intersect. 

Currently  there  is  no  flexibility  in  defining  the 
dependent  and  independent  variables  for  either  cartesian  or 
polar  coordinate  problems.  In  cartesian  coordinates, 
functions  are  of  the  independent  variable  x  and  the  dependent 
variable  is  y.  In  polar  coordinates  the  polar  function  r  must 
always  be  a  function  of  the  independent  variable  Q.  The 
program  could  be  improved  by  allowing  the  user  to  determine 
the  variable  names  (i.e.,  the  user  may  want  the  dependent 
variable  to  be  t) . 

The  focus  of  all  future  work  should  support  the  intent  of 
the  integral  package,  which  is  to  teach  the  calculus  student 
how  to  find  the  limits  of  integration  for  integral  problems. 
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As  such,  all  future  work  should  be  to  improve  upon  the 
package's  ability  to  do  this. 
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APPENDIX  A 


LAB  BOOK 


I .  INTRODUCTION 

One  of  the  most  difficult  topics  in  calculus  is  deciding 
on  the  limits  of  integration.  This  lab  book  shows  the 
calculus  student  how  a  given  domain  is  plotted  and  teaches  the 
student  how  to  find  the  limits  of  integration  when  evaluating 
two-dimensional  integrals. 

A.  Integration 

Integral  calculus  is  the  mathematics  we  use  to  find 
lengths,  areas,  and  volumes  of  irregular  shapes  and  to 
calculate  the  average  values  of  functions.  The  development  of 
integral  calculus  starts  from  the  calculation  of  areas. 


y 


We  can  use  the  integral  calculus  to  find  the  areas  of 
regions  like  the  shaded  one  here. 
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B.  Lab  Book  Organisation 

This  package  uses  numerous  example  problems  for  the 
two  dimensional  case.  It  considers  cartesian  as  well  as 
polar  coordinates.  Almost  all  of  the  problems  begin  with  a 
description  of  a  given  domain  over  which  an  integral  is  to  be 
evaluated.  The  given  domain  is  plotted  and  the  way  in  which 
the  limits  of  integration  are  found  is  shown.  Finally  the 
integrals  are  set  up  and  evaluated. 


II.  CARTESIAN  COORDINATE  SYSTEM 

A.  Pinding  the  Area  under  the  Graph  of  a  Monnnegative 
Continuous  Function. 

How  to  Find  the  Area  under  the  Graph  of  a 
Non-negative  Continuous  Function  y=f(x)  from  a  to 
b. 

STEP  1:  Find  an  antiderivative  F(x)  of  f(x). 

STEP  2:  Calculate  F(b)-F(a).  This  number  will  be 
the  area  under  the  curve  from  a  to  b. 

Example  1 

Find  the  area  under  the  curve  y=x2  from  x=0  to 
x-1 . 


y 
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Solution.  In  this  example  f(x)»x*,  a*0,  and  b=l . 

We  start  with  a  graph  of  the  function  and  identify 
the  area  of  the  region  that  lies  between  the  curve 
and  the  x-axis  from  a  to  b.  We  then  find  the  area 
in  two  steps. 


STEP  1:  Find  an  antiderivative  F(x)  of  f(x)=x2. 


x3 

F(X)  -  - 
3 

STEP  2:  Calculate  F(1)-F(0). 


( 1) 3  (0)3  1  1 

F  ( 1 )  -F  ( 0 )  - — 0  =  - 

3  3  3  3 


The  area  is  1/3. 

Example  2 

Find  the  area  under  one  arch  of  the  curve  y=cos  x. 


y 


Solution.  In  this  example  f(x)=cos  x,  however  a  and 
b  are  not  explictly  given.  We  graph  the  function 
and  observe  that  it  crosses  the  x-axis  at  x=-x/2  and 
x=x/2.  Note  that  cos(-rr/2)  =0  and  cos(»r/2)  =0. 
Thus  we  have  a=-x/2  and  b=x/2.  We  then  find  the 
area  in  two  steps. 
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STEP  1:  Find  an  antiderivative  F(x)  of  f(x)*cos  x. 
F(x)=sin  x 

STEF  2:  Calculate  F(-*/2) -F(x/2) . 

F (-w/2 ) -F(ir/2)  =sin(-ff/2) -sin(w/2)=l-(-l)=2 
The  area  is  2. 

B.  The  Average  Value  of  a  Function 

The  average  value  of  the  function  f  on  [a,b]  is  the 
integral  of  f  divided  by  the  length  of  the  integral. 


1  f 

Average  value  of  f  on  [a,b]  is  -  J  f(x)dx. 

b-a  a 


If  f  is  continuous  and  nonnegative  on  [a,b],  its 
average  value  is  the  height  of  a  rectangle  whose 
area. 


f (c) (b-a ) - 


J f (x)dx, 

a 


is  the  area  under  the  graph  of  f  from  a  to  b. 

y 
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Example  3 


Find  the  average  value  of  f ( x ) =  / 4  -  x2 ,  on 
the  on  the  interval  [-2,  2]. 

Solution:  We  graph  the  integrand  f(x)  over  the 

interval  of  integration  [-2,  2]  and  see  that  the 
graph  is  a  semicircle  of  radius  2.  The  area  between 
the  semicircle  and  the  x-axis  is 


1  1 

Area  =  (-)jr  r2  =-  (-)*  (2)2=2ir. 

2  2 

Because  the  area  is  also  the  value  of  the  integral 
of  f  from  x=-2  to  x=2,  the  average  value  of  f  on 
[-2,  2)  is 


Average  value  =■ 


2  —  ( —  2 ) 


|  \/ 4  -  X2  dx=  ( - )  2x= 

-2  4  2 


y  2 

y  *  4  -  x 


Example  4 

Find  the  average  value  of  f(x)=2x  on  the  interval 
[1,4]. 

Solution.  We  graph  the  integrand  f(x)  over  the 
interval  of  integration  [1,4]  and  we  see  that  the 
graph  of  the  area  is  a  trapezoid  whose  base  is  on 
the  x-axis.  The  area  of  a  trapezoid  is 
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1  1 

Area*  -  base(hl  +  h2)  =  -(3) (8  +  2)*15. 

2  2 

Because  the  area  is  also  the  value  of  the  integral 
of  f  from  x=l  to  x=4,  the  average  value  of  f  on  [1* 
4  ]  is 


Area  Between  Curves 

This  section  shows  how  to  find  the  area  of  a  region 
in  the  coordinate  plane  by  integrating  the  functions 
that  define  the  region's  boundaries. 

Definition 

If  functions  f  and  g  are  continuous  and  if  f(x)  is 
greater  than  or  equal  to  g(x)  throughout  the 
interval  [a,b]  such  that  x  is  in  [a,bj,  then  the 
area  of  the  region  between  the  curves  y=f(x)  and 
y=g(x)  from  a  to  b  is  the  integral  of  (f-g)  from  a 
to  b: 
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Area=  J  [f (x)-g(x) ]dx. 
a 

How  to  find  tho  Aroa  botwoon  Two  Curves 

1.  Graph  the  curves  together.  This  tells  you  which  is  f, 
(upper  curve)  and  which  is  g,  (lower  curve) .  It  also 
helps  to  find  the  limits  of  integration  if  you  do  not 
already  know  them. 

2.  Find  the  limits  of  integration. 

3.  Write  a  formula  for  f(x)-g(x).  Simplify  it  if  you 
can. 

4.  Integrate  f(x)-g(x)  from  a  to  b.  The  number  you  get  is 
the  area. 


Example  5 

Find  the  area  between  the  curve  y=cos  x  and  the 
curve  y=-sin  x  from  x=0  to  x=jt/2. 

Solution. 

STEP  1:  The  graphs.  We  graph  the  curves  together. 
The  upper  curve  is  y*cos  x,  so  we  take  f(x)=cos  x 
in  the  area  formula.  The  lower  curve  is  y=-sin  x, 
so  g ( x) =-sin  x. 


y 
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STEP  2:  The  limits  of  integration.  They  are 
given: 

a=0  and  b=*/2. 

STEP  3:  The  formula  for  f(x)-g(x).  From  Step  l, 
f (x) -g(x)=cos  x  -  (-sin  x) 

*cos  x  +  sin  x 

STEP  4:  Integrate  f(x)-g(x)  from  a=0  to  b=ir/2. 

K/2  x/2 

J  (cos  x  +  sin  x)dx  =  [sin  x  -  cos  x] 

0  0 

=  (  (l-O)-(O-l) ] =2 

The  area  between  the  curves  is  2. 

Example  6 

Find  the  area  between  the  curves  y=l  and  y=l-x'2 
for  x=l  to  x»4. 

Solution. 

8TEP  1:  The  graphs.  We  graph  the  curves  together. 
The  upper  curve  is  y«l,  so  f(x)=l.  The  lower  curve 

is  y~l-x  2 ,  so  g(x)=l-x'2. 


y 
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STEP  2:  The  limits  of  integration.  The  limits  of 
integration  are  given:  a=l  and  b=4 . 

STEP  3:  The  formula  for  f(x)-g(x). 

f(x)-g(x)=l-(l-x2)=x2 
STEP  4:  Integrate. 

Area=  [  f  (x)  -g (x)  ]dx=  Jx  2  dx=  -  - 

a  lx 

The  area  of  the  region  is  3/4. 

D.  Curves  that  Cross 

When  a  region  is  determined  by  curves  that  cross, 
the  crossing  points  give  the  limits  of  integration. 

Example  7 

Find  the  area  of  the  region  enclosed  by  the 
parabola  y=2-x2  and  the  line  y=-x. 

Solution. 

STEP  1:  The  graphs.  We  graph  the  curves  together. 

The  upper  curve  is  y=2-x2,  so  f(x)=2-xz.  The  lower 
curve  is  y=-x,  so  g(x)=-x.  The  x-coordinates  of  the 
points  where  the  parabola  and  line  cross  are  the 
limits  of  integration.  We  find  them  in  Step  2. 

y 
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STEP  2:  The  limits  of  integration.  We  find  the 

limits  of  integration  by  solving  the  equations 
y=2-x2  and  y=-x  simultaneously  for  x: 


2-x2~-x  (Equate  f(x)  and  g(x)) 

x2-x-2=0  (Transpose) 

(x+1) (x-2) =0  (Factor) 

x=-l ,  x=2.  (Solve) 


The  region  runs  from  x*-l  on  the  left  to  x*2  on  the 
right.  The  limits  of  integration  are  a*-l  and  b=2. 

STEP  3:  The  formula  for  f(x)-g(x). 

f  (x)-g(x)  =  (2-x2)-(-x)=2+x-x2 

STEP  4:  Integrate. 


Area=  J  ( f  (x) -g(x)  ]dx  =  J  (2+x-x2)dx 
a  -1 


x2  x3  2 
2x  +  -  -  - 

2  3-1 


4  8  11 

=  (4  + - )  -  (-2  +  -  +  -) 

2  3  2  3 


3  9  9 

=6  +  -  -  -  =  - 

2  3  2 

The  area  of  the  region  is  9/2. 
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Exaaple  • 

Find  the  area  of  the  reaion  enclosed  by  the  curve 
y-xVl  and  the  curve  y=x2+x. 

Solution. 

STEP  1:  The  graphs.  We  graph  the  curves  together. 

The  upper  curve  is  y=x3+l,  so  f(x)*x3+l.  The  lower 
curve  is  y=x2+x,  so  g(x)=x2+x.  The  x-coordinates 
of  the  points  where  the  curves  cross  are  the  limits 
of  integration.  We  find  them  in  Step  2. 


y 


STEP  2:  The  limits  of  integration.  We  find  the 

limits  by  solving  the  equations  y=x3+l  and  y=x  +x 


simultaneously  for  x: 

x3+l  =x2+x 

(Equate  f (x)  and  g(x) ) 

x3-x2-x+l  =0 

(Transpose) 

(X+l) (X-l)2  =0 

(Factor) 

X 

II 

M 

X 

II 

• 

(Solve) 

The  region  runs  from  x=-l  on  the  left  to  x=l  on  the 
right.  The  limits  of  integration  are  a=-l  and  b=l. 
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STEP  3:  The  formula  for  f(x)-g(x). 

f  (x)  -g(x)=x3+l-(x2+x)=x3-x2--x+l 
STEP  4:  Integrate. 

Area=  f  (x) -g(x)  ]dx=  }[x3-x2-x+l]dx 
a  -1 


111 
-x4--x3--x2  +  l 
4  3  2 


1 

-1 


1  1  1 

-(l)4  -  -(l)3  -  -(1)2+  l| 

4  3  2 


111 
-(-l)4  -  -(“D3  -  -(-l)2  +1 

4  3  2 


111 

- +1 

4  3  2 


111 

-  + - +  1 

4  3  2 


4 

3 


The  area  of  the  region  is  4/3 


Example  9 

Find  the  area  of  the  region  enclosed  by  the 
parabola  y=12-x2  and  the  parabola  y=4+x2  . 

Solution. 

STEP  1:  The  graphs.  We  graph  the  curves  together. 
The  upper  curve  is  y=12-x2,  so  f(x)=12-x2.  The  lower 
curve  is  y=4+x2,  so  g(x)=4+x2.  The  x-coordin-ates  of 
the  points  where  the  two  parabolas  cross  are  the 
limits  of  integration.  We  find  them  in  Step  2. 
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STEP  2:  The  limits  of  integration.  We  find  the 

limits  by  solving  the  equations  y=12-x2  and  y=4+x2 
simultaneously  for  x: 

4+x2  =  12-x2  (Equate  f(x)  and  g(x)) 

2x2  =  8  (Transpose) 

x2  *  4  (Divide  by  2) 

x=2,  x=-2 .  (Solve) 

The  region  runs  from  x=-2  on  the  left  to  x=2  on  the 
right.  The  limits  of  integration  are  a=-2  and  b=2. 

STEP  3:  The  formula  for  f(x)-g(x). 

f  (x)  -g (x)  =  12-x2-  (4+x2)  =8-2 x2 
8TEP  4:  Integrate. 

Area*  ^ [ f (x) -g (x) ]dx  =  f  [8-2x2]dx 
a  -2 
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2  2 

**  8x  -  -  X3 

3  -2 

2  2  62 

~[8(2)  -  - ( 2 ) 3 ]  -  [8 (-2)  -  - ( -2 ) 3] =  — 
3  3  3 


The  arei  of  the  region  is  62/3. 

E.  Boundaries  with  Changing  Formulas 

If  the  formula  for  one  of  the  bounding  curves 
changes  at  some  point  across  the  region,  you  may 
have  to  add  two  or  more  integrals  to  find  the  area. 

Example  10 

Find  the  area  of  the  region  in  the  first  quadrant 
bounded  above  by  the  curve  y=/~x  and  bounded  below 
by  the  x-axis  and  the  line  y=x-2. 

Solution. 

STEP  1:  The  graphs.  We  graph  the  curves  together. 
The  entire  upper  boundary  of  the  region  consists  of 
the  curve  y=/  x,  so  f(x)«/  x.  The  lower  boundary 
consists  of  two  curves,  first  y=0  for  x=0  to  x=2  and 
then  y=x-2  for  x=2  to  x=4.  Hence  the  formula  for 
g(x)  changes  from  g(x)=0  for  x=0  to  x=2 ,  to  g(x)=x-2 
for  x=2  to  x=4 . 
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STEP  2:  The  limits  of  integration.  The  limits  of 
integration  for  the  pair  f(x)=/~x  and  g(x)=0  are  a=0 
and  b=2.  For  the  pair  f(x)=J~x  and  g(x)=x-2,  the 
left-hand  limit  is  a=2  and  the  right-  hand  limit  is 
the  x-coordinate  of  the  upper  point  where  the  line 
crosses  the  parabola.  To  find  it,  we  solve  the 
equations  y=/x  and  y=x-2  simultaneously  for  x: 


y~x=x-2  (Equate  f(x)  and  g(x)) 

x= ( x-2 ) 2 3  =x2-4x+4  (Square) 

x2-5x+4=0  (Transpose) 

(x-1) (x-4) =0  (Factor) 

x=l,  x=4.  (Solve) 

The  value  x=l  does  not  satisfy  the  equation  /"x=x- 2. 
It  is  an  extraneous  root  introduced  by  squaring. 
The  x=4  gives  our  upper  limit  of  integration. 

STEP  3:  The  formulas  for  f(x)-g(x). 

For  x=0  to  x=2 :  f (x) -g(x) =f~x-0=J~x, 

For  x=2  to  x=4 :  f (x) -g (x) =/~x- (x-2) 

=  /~x-x+2 . 

STEP  4:  Integrate.  We  have  two  integrals  to 

evaluate.  Their  sum  is  the  area. 


Area=  J [f (x)-g(x) ]dx 
a 


J  J~x  dx  +  J[,/"“x-x+2]dx 
0  2 


2  3/2  2  2  3/2  X2  4 

=  -x  +  —  X  -  -  +  2x 

3  0  3  2  2 
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2  3/2 

2  3/2  ( 4 ) 2  1 

12  3/2  (2)2  1 

«  -(2)  + 

-(4) 

-  -  +  2(4) 

- 

-(2) 

-  ~  +  2(2) 

3 

3 

2 

3 

2 

10 
=r  — — 

3 

The  area  of  the  region  is  10/3. 

Exaaple  11 

Find  the  area  of  the  region  in  the  first  quadrant 
bounded  above  by  the  curves  y=x  and  y=2-x2  and  below 
by  the  x-axis. 

Solution. 

STEP  1:  The  graphs.  We  graph  the  curves  together. 
The  upper  boundary  consists  of  two  curves,  first 

y=x  for  x=0  to  x=l  and  then  y=2-x2  for  x=l  to 
x=/  2.  Hence  the  formula  for  f(x)  changes  from 
f(x)=xfor  x=0  to  x-1 ,  to  f(x)=2-x2  for  x=l  to  x~/  2. 
The  entire  lower  boundary  of  the  region  consists  of 

the  curve  y=0 ,  so  g(x)=0. 

y 


STEP  2:  The  limits  of  integration.  For  the  pair 
f(x)=x  and  g(x)=0  the  left-hand  limit  is  a=0.  The 
right-hand  limit  is  the  x-coordinate  of  the  upper 
point  where  the  line  crosses  the  parabola.  To  find 
it,  we  solve  the  equations  y=x  and  y=2-x2  simul¬ 
taneously  for  x: 
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X=2-X2 


(Equate  f(x)  and  g(x)) 

x2+x-2=0  (Transpose) 

( x+2 ) ( x— 1 ) — 0  (Factor) 

x=-2 ,  x=l.  (Solve) 

Since  x  is  greater  than  or  equal  to  zero  in  the 
first  quadrant  we  must  pick  x=l.  The  limits  of 
integration  are  a=0  and  b*l.  For  the  pair 
f(x)=2-x2  and  g(x)=0  the  left-hand  limit  is  1  and 
the  right-hand  limit  is  obtained  by  solving  for 
where  the  curve  y=2-x2  crosses  the  x-axis: 


2-x2— 0  (Equate  f(x)  and  g(x)) 

2— x2  (Transpose) 

x»y— 2.  (Solve) 

This  region  runs  from  1  on  the  left  to  J~ 2  on  the 
right.  The  limits  of  integration  are  a=l  and  b=vr_2. 


STEP  3:  The  formulas  for  f(x)-g(x). 

For  x=0  to  x=l:  f (x) -g (x) =x-0=x, 

For  x=l  to  x=J~~2 :  f (x)  -g(x)=2-x2-0=2-x2. 

STEP  4:  Integrate.  We  have  two  integrals  to 
evaluate.  Their  sum  is  the  area. 


Area=  J [f (x) -g(x) ]dx 


=  }  xdx  +  J 


[2-x2]dx 


X2 

1 

Xs 

=  - 

+ 

2x  -  - 

2 

0 

3 

(l)2 

-  +  [2/~2 
2 


n 

i 

(</  2 ) 3 


]  -  [2(1)  - 


(l)3 

-  ] 

3 


=  .719 
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Ixaipla  12 


Find  the  area  bounded  by  the  parabola  y2  =x  and  the 
line  y»6-x. 

Solution. 

STEP  1:  The  graphs.  We  graph  the  curves  together. 
The  upper  boundary  consists  of  two  curves,  first 
y=y  x  for  x=0  to  x=4  and  then  y=6-x  for  x=4  to 
x=*9 .  Hence  the  formula  for  f(x)  changes  from 
for  x=0  to  x=4  to  f(x)=6-x  for  x=4  to 
x=9.  The  entire  lower  boundary  of  the  region 
consists  of  y=-J~x,  so  g ( x) =-/x . 


y 


STEP  2:  The  limits  of  integration.  For  the  pair 
f (x)=/^c  and  g(x)=-/"lc  the  left-hand  limit  is  0. 
The  right-hand  limit  is  the  x-coordinant  of 
the  upper  point  where  the  line  crosses  the 
parabola.  To  find  it,  we  solve  the  equations 
y -y~x  and  y=6-x  simultaneously  for  x: 

7~x=6-x  (Equate  f(x)  and  g(x)) 

x=x2-12x+36  (Square) 

x2-13x+36=0  (Transpose) 

( x-4 ) (x-9)=0  (Factor) 

x=4 ,  x=9.  (Solve) 
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The  upper  point  where  the  parabola  and  line  cross 
is  at  x=4.  The  limits  of  integration  are  a=0  and 
b=4.  For  the  pair  f (x) =6-x  and  q (x) =-J~x  the 
left-hand  limit  is  x=4  and  the  right-hand  limit  is 
the  x-coordinate  of  the  lower  point  where  the 
parabola  and  line  cross.  We  found  this  point  when 
we  solved  the  two  equations  above  simultaneously 
for  x.  The  right-hand  limit  is  x=9.  The  limits  of 
integration  are  a=4  and  b=9. 

STEP  3:  The  formulas  for  f(x)-g(x). 

For  x=0  to  x=4:  f (x) -q(x)=J~x- (-/~x) 

-2j~x , 

For  x=4  to  x=9:  f (x) -g (x) =6-x- (-J~x) 

=6 -x+/~x. 

STEP  4:  Integrate.  We  have  two  integrals  to 
evaluate.  Their  sum  is  the  area. 

Area=^[ f (x) -g(x) ]dx=  ^2/“x  dx  +  J[6-x+/“x]dx 
a  0  4 

4  3/2  4  X2  2  3/2  9 

=  -  x  +  6X--  +  -X 

3  0  2  3  4 

4  3/2  (9) 2  2  3/2  ( 4 ) 2  2  3/2 

=-(4)  +[6(9)  -  -  +  -(9)  ]  -  [6(4)  -  -  +  -(9)  ) 

3  2  3  2  3 

=20.833 

The  area  of  the  region  is  20.833. 

F.  Integrating  with  Respect  to  y 

When  a  region's  bounding  curves  are  described  by 
giving  x  as  a  function  of  y,  the  basic  formula 
changes , 

For  regions  like  these. 
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y 


use  the  formula: 


Area= 


J [f (y)-g(y) ]<*y- 

a 


The  only  difference  is  that  we  are  now  integrating  with 
respect  to  y  instead  of  x.  We  can  sometimes  save  time  by 
doing  so.  The  basic  steps  are  the  same  as  before. 

Example  13 

Find  the  area  of  the  region  between  the  curves  x=y2  and 
x=y+2  in  the  first  quadrant. 

Solution. 

STEP  1:  The  graphs.  We  graph  the  curves  together.  The 
right-hand  curve  is  x=y+2,  so  f(y)=y+2.  The  left-hand 
curve  is  x=y2,  so  g(y)=y2. 
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STEP  2:  The  limits  of  integration.  The  lower  limit  of 
integration  is  y=0.  The  upper  limit  is  the  y-coordinate 
of  the  upper  point  where  the  line  crosses  the  parabola. 
We  find  it  by  solving  the  equations  x=y+2  and  x=y2 
simultaneously  for  y: 

y+2=y2  (Equate  f(y)  and  g(y)) 
y2+y+2=0  (Transpose) 


(y+1) (y-2)=0  (Factor) 


y=-l ,  y=2.  (Solve) 


The  upper  limit  of  integration  is  2.  (The  value  y=-l 
gives  the  point  of  intersection  below  the  x-axis.) 


step  3:  The  formula  for  f(y)-g(y). 
f (y)-g(y)*y+2~y2 

STEP  4:  Integrate. 

Area=  ? ( f (y) -g(y) ]dy=  f[2+y-y2]dy 
a  0 


y2  y3 

2y  +  -  -  - 
2  3 


2  4  8 

=  4  +  -  -  - 

0  2  3 


10 

3 


The  area  of  the  region  is  10/3. 

Example  14 

Find  the  area  of  the  region  between  the  curves  x=y 
and  x-J  2 -y  in  the  first  quadrant. 

Solution. 


STEP  1:  The  graphs.  We  graph  the  curves  together.  The 
right-hand  curve  is  x=7  2-y ,  so  f(y)=/  5-y.  The  left-hand 
curve  is  x=y ,  so  g(y)=y. 
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y 


X 

STEP  2:  The  limits  of  integration.  The  lower  limit  of 
integration  is  y=0.  The  upper  limit  is  the  y-coordinate 
of  the  upper  point  where  the  line  crosses  the  parabola. 
We  find  it  by  solving  the  equations  x=y  and  x~J  2-y 
simultaneously  for  y: 

y =7 2-y  (Equate  f(y)-g(y)) 

y2 3=2-y  (Square) 

y2+y-2=0  (Transpose) 

(y+2)(y-l)=0  (Factor) 

y=-2,  y=l.  (Solve) 

The  upper  limit  of  integration  is  1.  (The  value  y=-l  gives 
the  point  of  intersection  below  the  x-axis.)  The  limits 
of  integration  are  a=0  and  b=l. 

STEP  3:  The  formula  for  f(y)-g(y). 

f(y)-g(y)=7  2-y  -  y 

STEP  4:  Integrate. 


Area=  J[f  (y) -g(y)  Jdy  =  J[S~T= y  -  y)dy 
a  0 


2  3/2  y2  1 

=  -  -(2-y)  -  - 

3  2  0 
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2  3/2  (l)2  212  3/2 

=[-  -(2-1)  -  -]=-  -  +  -(2)  -  0= .719 

3  2  3  2  3 

The  area  of  the  region  is  .719. 

Exaaple  15 

Find  the  area  bounded  by  the  parabola  y2=x  and  the  line 
y=6-x . 

Solution. 

STEP  1:  The  graphs.  We  graph  the  curves  together.  The 
right-hand  curve  is  x=6-y,  so  f(y)=6-y.  The  left-hand 
curve  is  x=y2,  so  g(y)=y2. 


y 


STEP  2:  The  limits  of  integration.  The  lower  limit  of 
integration  is  the  y-coordinate  of  the  lower  point  where 
the  parabola  crosses  the  line.  The  upper  limit  is  the 
y-coordinate  of  the  upper  point  where  the  parabola  and  the 
line  cross.  We  find  these  points  by  solving  the  equations 
x=y2  and  x=6-y  simultaneously  for  y: 

y2=6-y  (Equate  f(y)  and  g(y)) 

y2+y-6=0  (Transpose) 

(y+3) (y-2)=0  (Factor) 

y=-3,  y=2.  (Solve) 

The  upper  limit  of  integration  is  2  while  the  lower  limit 
is  -3.  The  limits  of  integration  are  a=-3  and  b=2. 

STEP  3:  The  formula  for  f(y)-g(y). 

f (y)-g(y)=6-y-y2 
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STEP  4:  Integrate. 


Area= 


J [f (y)-g(y) ]dy 

a 


[6_y_y2]dy 


6y 


y 2  y3 

2  3 


=[6(2) 


(2) 2  ( 2 ) 3  (~3) 2  (-3)3 

-  -  -  ]  -  [  6  ( -3 )  - —  ] 

2  3  2  3 


8  9 

=[12  -  2  -  -  [-18  -  +  9]  =  20.833 

3  2 

The  area  bounded  by  the  curves  is  20.833. 
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III.  POLAR  COORDINATE  SYSTEM 


This  section  shows  how  to  calculate  areas  of  plane 
regions  and  lengths  of  curves.  The  general  methods  for 
setting  up  the  integrals  are  the  same  as  for  cartesian 
coordinates,  although  the  resulting  formulas  are  somewhat 
different. 

A.  Ares  in  the  Plane 


We  calculate  areas  of  plane  regions  like  the  shaded  one 
above  using  polar  coordinates. 

Area  Between  the  Origin  and  r=f(Q),  Q  between  a  and  b: 


Area  =  j  -  r2  dQ. 
a  2 

Example  16 

Find  the  area  of  the  region  enclosed  by  the  cardioid 
r=2 ( 1+cosQ) . 

Solution. 

STEP  1:  The  graph.  We  graph  the  cardioid  and  determine 
that  the  radius  OP  sweeps  out  the  region  exactly  once  as 
Q  runs  from  0  to  2w . 
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STEP  2:  The  limits  of  integration.  The  limits  of 

integration  are  a=0  and  b=2x. 

STEP  3:  Determine  the  integrand. 

1  1 

-  r2  =  -(4) (1+cosQ)2  =  2 ( l+2cosQ+cos?Q) 

2  2 

1+COS2Q 

=  2+4cosQ+2( - J  =  3  +  4cosQ+cos2Q 

2 

STEP  4:  Integrate. 


Jjr  sin2Q  2if 

[ 3+4cosQ+cos2Q]dQ  =  3Q+4sinQ+[ - ] 

0  2  0 

=  6w+0+0=6ir 

The  area  of  the  region  is  6ir. 

Example  17 

Find  the  area  of  the  region  bounded  by  the  rays  Q=-*r/4  and 
Q=»r/4,  and  the  graph  of  the  cardioid  r=l+sinQ. 

Solution. 

STEP  1:  The  graph.  The  region  is  the  portion  of  the 
cardioid  swept  out  by  the  radius  of  length  r=l+sinQ  as  Q 
increases  from  Q=-ir/4  to  Q=jt/4. 
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y 


8TEP  2:  The  limits  of  integration.  The  limits  are  given: 

a=-w/4  and  to=w/4. 

STEP  3:  Determine  the  integrand. 

11  1 
-  r2  =  -[1+sinQ]2  =  - [ l+2sinQ+sin2Q] 

2  2  2 

1  11 

=  -[l+2sinQ+(-  -  -cos2Q) ] 

2  2  2 

13  1 

=  - [ -  +  2sinQ  -  -cos2Q] 

2  2  2 

3  1 

=  -  +  sinQ  -  -cos2Q 

4  4 

STEP  4:  Integrate. 


1  w/4  3  1 

-  rz  dQ  =  J  [-  +  sinQ  -  -cos2Q]dQ 
a  2  -w/4  4  4 

3  1 
=  -Q  -  cosQ  -  -sin2Q 

4  8 

The  area  of  the  region  is  .928. 

Example  18 

Find  the  area  inside  the  smaller  loop  of  the  limacon 
r=2cosQ+l . 


w/4 

=  .928 

-w/4 
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Solution. 


8TBP  1:  The  graph.  After  sketching  the  curve  we  see  that 
the  smaller  loop  is  traced  out  by  the  point  (r,Q)  as  Q 
increases  from  Q=2ir/3  to  Q=4x/3. 

y 


STEP  2:  The  limits  of  integration.  Since  the  curve  is 
symmetric  about  the  x-axis,  we  may  calculate  the  area  of 
the  shaded  half  of  the  inner  loop  by  integrating  from 
Q=2*/3  to  Q=x.  The  limits  are  then  a=2x/3  and  b=ir.  The 
area  we  seek  will  be  twice  the  value  of  the  resulting 
integral : 


Area=  2 


J 

2»/3 


-  rz  dQ  =  J  r2  dQ. 
2  2ff/3 


STEP  3:  Determine  the  integrand. 


r2  =  (2cosQ  +  l)2  =  4cos2Q  +  1 

1+COS2Q 

*  4  [  -  — - ]  +  4cosQ  +  1 

2 

=  2  +  2cos2Q  +  4cosQ  +  1 
=  3  +  2cos2Q  +  4cosQ 


STEP  4:  Integrate. 

n 

Area=  J  [ 3+2cos2Q+4cosQ]dQ 
27T/3 
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w 


3Q+sin2Q+4sinQ 


2ir/3 


r*  4  n 

■  3»  -  [2 IT  -  -  +  - ] 

2  2 


3/“T 


=  w  - 


The  area  of  the  smaller  loop  is  ir  -  3 J  3/2. 

Example  19 

Find  the  area  of  the  region  enclosed  by  the  graph  of  the 
polar  equation  r=sin3Q. 

Solution. 


STEP  1:  The  graph.  When  we  graph  the  region  we  see  that 
it  is  a  three-leaved  rose. 
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STEP  2:  Limits  of  integration.  Since  the  three  leaves 
are  congruent,  and  since  the  leaf  determined  by  Q  from  Q=0 
to  Q=ir/3  is  symmetric  about  the  ray  Q=*/6,  we  may  obtain 
the  area  of  the  figure  by  calculating  the  area  of  the  half 
leaf  determined  by  Q=0  to  Q=ir/6  and  multiplying  the  result 
by  6.  The  limits  of  integration  are  a=0  and  b=x/6. 


w/6  1 

Area=  6  J  -  r2  dQ 
0  2 

STEP  3:  Determine  the  integrand. 

1 

6  [-  r2]  =  3r2  =  3  (sin3Q) 2 
2 

11  3  3 

=  3 [ -  -  -cos6Q]  =  -  -  -cos6Q 

2  2  2  2 

STEP  4:  Integrate. 


n/6  3  3  3  1  jt/6  k 

Area-  J  [ -  -  -cos6Q]dQ  =  -Q  -  -sin6Q  =  - 
0  2  2  2  4  0  4 

The  area  of  the  region  is  <r/4. 

Example  20 

Find  the  area  of  the  region  enclosed  by  the  lemniscate  r2 
=cos2Q. 

Solution. 

STEP  1:  The  graph.  By  extracting  square  roots  on  both 
sides  of  the  equation  we  obtain  r =/  cos2Q  and  we  graph  the 
figure.  v 
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STEP  2:  The  limits  of  integration.  We  note  that  the 

function  r  is  defined  only  for  Q  between  Q=-»/4  and  Q=ir/4 
and  for  Q  between  Q*3ir/4  and  Q=-*/4.  As  Q  ranges  through 
these  two  intervals  the  two  lobes  of  the  lemniscate  are 
traced  out.  The  rays  Q=-jr/4  and  Q-w/4  determine  half  the 
area  of  the  lemniscate.  We  may  therefore  integrate  from 
Q=-ir/4  to  Q=r/4  to  obtain  the  area  of  one  lobe  and  double 
the  result.  The  limits  are  then  a=-ir/4  and  b=»/4. 

STEP  3:  Determine  the  integrand. 

r2  =  cos2Q 

8TBP  4:  Integerate. 


Area= 


1 

cos2QdQ  =  -sin2Q 


w/4 


=  1 


-*V  4 

The  area  of 

B.  Area  of 


2 


-*/ 4 


the  region  is  1. 

Region  which  lies  Between  two  Polar  Curves 


y 


To  find  the  area  of  a  region  which  lies  between  two  polar 
curves  from  Q=a  to  Q=b,  we  subtract  the  integral  of 
(1/2) rl2  dQ  from  the  integral  of  (1/2) r22  dQ.  This  leads 
to  the  following  formula. 

Area  of  the  Region:  rl(Q)  less  than  or  equal  to  r  and  r 
less  than  or  equal  to  r2(Q),  Q  from  Q=a  to  Q=b 
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b  i  b  1  b  1 

Area=  J  -  r22  dQ  -  J  -  rl2  dQ  =  J  -(r22  -rl2)dQ. 
a  2  a  2  a  2 

Example  21 

Find  the  area  of  the  region  that  lies  inside  the  circle 
r=l  and  outside  the  cardiod  r=l-cosQ. 

Solution. 

STEP  1:  The  graphs.  We  sketch  the  region  to  determine 
its  boundaries  and  find  the  limits  of  integration.  The 
outer  curve  is  r2=l,  and  the  inner  curve  is  rl=l-cosQ. 

y 


STEP  2:  Limits  of  integration.  We  see  that  Q  runs  from 
-k/2  to  k/2.  Because  of  symmetry  we  find  the  area  for  Q=0 
to  Q=ir/2.  The  limits  of  integration  are  a=-*/2  and  b=ir/2. 


k/2  1  k/2  1 

Area=  J  -(r22  -rl2)dQ  =  2  J  -(r22  -rl2)dQ 
-K/2  2  02 


k/2 

=  J  (r22  -  rl2)  dQ 
0 

8TEP  3:  Determine  the  integrand. 


r22-rl2  =  1- ( l-2cosQ+cos2Q)  =  2cosQ-cos2Q 
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COS2Q+1 


=  2cosQ  - 


2 


STEP  4:  Integrate. 


ir/2  cos2Q+l 

Area=  J[2cosQ  -  - ]dQ 

0  2 


2sinQ  - 


sin2Q 


jr/2 

0 


* 

-  2  -  - 
4 


The  area  of  the  region  is  2  -  ir/4. 

Example  22 

Find  the  area  of  the  region  lying  inside  the  circle 
r=3cosQ  and  outside  the  cardioid  r=l+cosQ. 

Solution. 

STEP  lj  The  graphs.  We  sketch  the  graphs  of  the  two 
equations  r2=3cosQ  and  rl=l+cosQ. 

y 


STEP  2:  The  limits  of  integration.  By  solving  the  two 
equations  simultaneously,  we  find  the  point  of  inter¬ 
section: 
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3cosQ=l+cosQ 


(Equate  rl  and  r2) 
(Transpose) 


2 cosQ* 1 
1 

cosQ*-  (Divide) 

2 

Q=ir / 3  ,  Q=-*/3  .  (Solve) 

The  limits  of  integration  are  a=-jr/3  and  b=ir/3. 
STEP  3:  Determine  the  integrand. 


- (r22  -  rl2)  =  - [  ( 3cosQ) 2  -  (1  +  cosQ)2] 

2  2 

1 

=  -[9cos2Q  -(1  +  2cosQ  +  cos2Q)  ] 
2 

1 

=  -[8cos2Q  -  1  -  2COSQ] 

2 

1 

=  4cos2Q  -  -  -  cosQ 
2 

11  1 

=  4 [ -  +  -cos2Q]  -  -  -  cosQ 
2  2  2 

1 

=  2  +  2cos2Q  -  -  -  cosQ 
2 

3 

=  -  +  2cos2Q  ~  cosQ 
2 

STEP  4:  Integrate. 
tt/3  3 

Area=  J  [-  +  2cos2Q  -  cosQ]dQ 
-ir/ 3  2 

3  ir/3 

=  -Q  +  sin2Q  -  sinQ 

2  -n/3 


=  Pi 


The  area  of  the  region  is  w. 

C.  The  Length  of  a  Curve 

We  calculate  the  length  of  a  curve  r=f(Q),  Q  from  a  to  b 
by  expressing  the  differential  ds=y~(dx2  +dy2)  in  terms  of 
Q  and  integrating  from  a  to  b. 

If  r=f(Q)  has  a  continuous  first  derivative  for  Q  from  a 
to  b  and  if  the  point  P(r,Q)  traces  the  curve  r=f(Q) 
exactly  once  as  Q  runs  from  a  to  b,  then  the  length  of  the 
curve  is  given  by: 


b  /  dr2 

Length=  J  V  r2  +  ( — )  dQ. 
a  dQ 

Example  23 

Find  the  length  of  the  cardioid  r=l-cosQ. 

Solution. 

STEP  l:  The  graph.  We  sketch  the  cardioid  to  determine 
the  limits  of  integration.  The  point  P(r,Q)  starts  at  the 
origin  and  traces  the  curve  once,  counterclockwise. 

y 

P(r.Q) 


STEP  2:  Limits  of  integration.  We  see  that  Q  runs  from 
0  to  2 w,  so  the  limits  of  integration  are  a=0  and  b=2n . 

STEP  3:  Determine  the  integrand. 

dr 

r=l-cosQ,  —  =  sinQ, 
dQ 


dr2 

r2  +  (— )  ®(l-cosQ)2  +(sinQ)2  =l-2cosQ+cos2Q+sin2Q 
dQ 

=1-2cosQ+1=2-2cosQ 
STEP  4:  Integrate. 


Length®  Jx/r2  +  (  — )  dQ  -  J 
a  v  dQ  0 


2-2cosQ  dQ 


The  length  of  the  curve  is  8. 
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APPENDIX  B 


TUTORIAL 


A.  INTRODUCTION 

This  tutorial  instructs  the  user  on  how  to  start  Mathe- 
matica  and  how  to  load  and  use  the  program.  The  tutorial 
includes  three  example  problems  that  are  intended  to  show  the 
user  how  the  program  works.  The  first  example  problem  is  an 
integral  problem  given  in  cartesian  coordinates.  In  the  first 
problem,  the  program  is  used  to  find  the  area  of  a  region  by 
first  integrating  with  respect  to  the  variable  x  and  then  with 
respect  to  the  variable  y.  The  second  problem  is  an  integral 
problem  given  in  polar  coordinates.  In  this  problem  the 
program  finds  the  area  between  a  polar  curve  and  the  origin. 
The  final  example  problem  is  also  given  in  polar  coordinates. 
In  this  problem  the  program  is  used  to  find  the  area  of  the 
region  bounded  by  two  polar  curves.  These  three  example 
problems  are  representative  problems  that  the  program  can 
solve.  They  show  some  of  the  program's  many  capabilities. 
The  tutorial,  while  not  comprehensive,  gives  the  user  some 
idea  of  the  types  of  problems  the  program  can  solve.  It 
demonstrates  how  the  program  prompts  the  user  for  input  and 
how  that  input  is  entered. 
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B.  STARTING  MATHEMATICA  AMD  LOADING  PROGRAM 

To  start  Mathematica,  double-click  the  Matheaatica  icon. 
Mathematica  documents  are  called  Notebooks.  When  you  start 
Mathematica,  an  empty  Notebook  window  appears  on  your  screen. 
As  soon  as  you  start  typing,  a  cell  is  created,  indicated  by 
a  cell  bracket  along  the  right  edge  of  the  Notebook  window. 
Each  piece  of  information  in  a  Notebook  is  contained  in  a 
separate  cell.  To  load  the  program  and  start  the  kernel  first 
type:  <<integral.ma. 

You  now  must  activate  the  cell.  There  are  several  ways  to 
do  this.  You  can  use  the  Mouse  and  select  Evaluate  Selection, 
in  the  Action  menu  which  appears  at  the  top  of  the  screen  (see 
Figure  18  on  page  96) .  The  three  keyboard  equivalents  are: 

1.  Shift+Enter 

2 .  Insert 

3.  5  on  the  numeric  keypad. 

Once  you  have  activated  the  cell  you  are  ready  to  use  the 
program. 
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Mathematics  -  (Newnb-1 1 


Figure  18 
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C.  EXAMPLE  1 


Find  the  area  of  the  region  bounded  by  the  functions  y  = 
5x,  y  =  x  and  y  =  4  -  x2  in  the  first  quadrant. 

*  Press  the  down  arrow  to  create  a  new  cell  and  then  type 
graph . 

*  For  this  problem  press  the  Insert  key  to  activate  the 
cell. 

*  You  will  see  appear  the  message: 

To  display  a  graph  that  is  in  the  lab  book  type  the 
example  number,  otherwise  type  n. 

Because  you  are  going  to  enter  our  own  example  type  n  and 
press  Enter. 

*  You  will  see  the  message: 

Type  1  for  cartesian  coordinate  system,  or  type  2  for 
polar  coordinate  system. 

Since  your  first  example  is  a  problem  given  in  cartesian 
coordinates  type  1  and  press  Enter. 

*  You  will  see  the  message: 

To  just  plot  curves  type  y,  otherwise  type  n. 

Since  you  do  not  want  to  just  plot  the  functions  type  n. 

*  You  will  see  the  message: 

Enter  the  functions  which  define  the  region's 
boundaries.  When  finished  type  n. 

The  functions  must  be  functions  of  the  independent 
variable  x  and  are  entered  one  at  a  time.  They  can  be 
inputted  in  any  order.  You  are  now  ready  to  enter  the 
functions.  Type  5x  then  press  Enter,  type  x  then  press 
Enter,  type  4  -  x2  then  press  Enter. 

*  You  will  see  the  message: 

If  any  of  the  limits  are  known  type  y,  otherwise  type  n. 

Since  none  of  the  limits  are  explicitly  given  you  type  n. 

*  You  will  see  the  message: 
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If  the  region  is  restricted  to  ths  right  half  plans 
typs  y ,  otherwise  typs  n. 

Because  we  are  told  in  the  problem  statement  that  the 
region  is  restricted  to  the  first  quadrant  type  y. 

*  The  region  is  then  plotted  and  the  program  solves  the 
problem.  You  will  see  the  following  output: 


y 


The  limits  of  integration  for  the  first  region  are  x  =  0 
and  x  =  0.702 

The  upper  curve  for  the  first  region  is  y  =  5  x 

The  lower  curve  for  the  first  region  is  y  =  x 

The  area  of  the  first  region  is  0.985608 

The  limits  of  integration  for  the  second  region  are  x  = 

0.702  and  x  =  1.562 

The  upper  curve  for  the  second  region  is  y  =  4  -  x2 

The  lower  curve  for  the  second  region  is  y  =  x 

The  area  of  the  second  region  is  1.31145 
The  area  of  the  total  region  is  2.29706 

All  of  the  functions  are  invertible  and  and  it  is 
appropriate  to  integrate  the  problem  with  respect  to  the 
variable  y.  As  a  result,  after  a  breif  delay  you  will 
see  the  message: 

Type  y  if  you  wish  to  integrate  the  same  problem  with 
respect  to  y,  otherwise  type  n. 
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Type  y  and  you  will  see  the  following  output: 

The  limits  of  integration  for  the  first  region  are  y  =  0 
and  y  =  1.56 

The  right  curve  for  the  first  region  is  x  =  y 

y 

The  left  curve  for  the  first  region  is  x  = 

5 

The  area  of  the  first  region  is  0.973443 

The  limits  of  integration  for  the  second  region  are  y  = 
1.56  and  y  =  3.51 

The  right  curve  for  the  second  region  is  x=Sqrt[4-y] 

y 

The  left  curve  for  the  second  region  is  x  = 

5 

The  area  of  the  second  region  is  1.32362 
The  area  of  the  total  region  is  2.29706 


D.  EXAMPLE  2 

Find  the  area  of  one  leaf  of  the  three-leaf  rose  formed  by 
the  polar  curve  r  -J  sin3Q. 

*  Press  the  down  arrow  to  create  a  new  cell  and  then  type 

graph . 

*  Activate  the  cell  by  using  the  Mouse  to  select  Evaluate 
Selection  from  the  Action  menu. 

*  You  will  see  appear  the  message: 

To  display  a  graph  that  is  in  the  lab  book  type  the 
example  number,  otherwise  type  n. 

Because  you  are  going  to  enter  our  own  example  type  n  and 
press  Enter. 

*  You  will  see  the  message: 

Type  l  for  cartesian  coordinate  system,  or  type  2  for 
polar  coordinate  system. 

Since  your  second  example  is  a  problem  given  in  polar 
coordinants  type  2  and  press  Enter. 

*  You  will  see  the  message: 

To  just  plot  curves  type  y,  otherwise  type  n. 

Since  you  do  not  want  to  just  plot  the  polar  curve  type  n. 
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*  You  will  see  the  message: 

Type  1  for  ares  between  the  origin  and  a  polar  curve,  or 
type  2  for  area  between  two  polar  curves. 

Type  1  since  your  problem  is  to  fii.d  the  area  between  the 
origin  and  the  polar  curve  r  =y  sin  3q  for  one  leaf. 

*  You  will  see  the  message: 

Enter  the  polar  curve  which  bounds  the  region. 

Type  eqrt[sin[3Q] ] . 

*  You  will  see  the  message: 

If  the  limits  are  known  type  y,  otherwise  type  n. 

Since  the  limits  are  not  known  type  n. 

*  The  polar  curve  is  then  plotted.  You  will  see  the 
following  output: 

y 


The  polar  curve  intersects  the  origin  at  Q: 

Pi  2  Pi  4  Pi  5  Pi 

(0,  ,  i  Pi,  ,  ,  2Pi ) 

3  3  3  3 

After  a  brief  delay  you  will  see  the  message: 
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The  polar  curva  is  plotted  from  Q=o  to  Q=2  Pi.  if  you  want 
to  change  the  limits  type  y ,  otherwise  type  n. 

Since  you  are  only  interested  in  the  area  of  one  leaf 
type  y. 

*  You  will  see  the  message: 

Input  new  lover  limit. 

Type  0. 

*  You  will  see  the  message: 

Input  new  upper  limit. 

Type  pi/ 3. 

*  The  region  in  which  you  are  interested  is  then  plotted. 
You  will  see  the  following  output: 


The  limits  of  integration  for  the  regior  which  is  bounded 
by  the  polar  curve  r  =  Sqrt[sin[3Q] ]  and  the  origin  are: 


Pi 

a=0  and  b= — 

3 

The  area  bounded  by  the  polar  curve  and  the  origin  is 
.333 
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B.  EXAMPLE  3 

Find  the  area  of  the  region  bounded  by  the  circle  r  *  1 
and  the  cardioid  r  =  1  -  cos  Q  if  the  cuter  curve  is  the 
cardioid  and  the  inner  curve  is  the  circle. 

*  Press  the  down  arrow  to  create  a  new  cell  and  then  type 

graph . 

*  Activate  the  cell  by  pressing  Shift+Enter. 

*  You  will  see  appear  the  message: 

To  display  a  graph  that  is  in  the  lab  book  type  the 
example  number,  otherwise  type  n. 

Because  you  are  going  to  enter  our  own  example  type  n  and 
press  Enter. 

*  You  will  see  the  message: 

Type  1  for  cartesian  coordinate  system,  or  type  2  for 
polar  coordinate  system. 

Since  your  third  example  is  a  problem  given  in  polar 
coordinates  type  2  and  press  Enter. 

*  You  will  see  the  message: 

To  just  plot  curves  type  y,  otherwise  type  n. 

Since  you  do  not  want  to  just  plot  the  polar  curves  type 

n. 

*  You  will  see  the  message: 

Type  l  for  area  between  the  origin  and  a  polar  curve,  or 
type  2  for  area  between  two  polar  curves. 

Type  2  since  your  problem  is  to  find  the  area  between  two 
polar  curves. 

*  You  will  see  the  message: 

Enter  one  of  the  polar  curves  which  bounds  the  region. 
Type  1  -  cos [Q] . 
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*  You  will  see  the  message: 

Enter  the  other  polar  curve  which  bounds  the  region. 

Type  1 . 

*  You  will  see  the  message: 

If  the  Units  are  known  type  y,  otherwise  type  n. 

Since  we  do  not  know  where  the  two  curves  intersect  and 
what  the  limits  of  integration  are  type  n. 

*  The  two  polar  curves  will  then  be  plotted.  You  will  see 
the  following  output: 


y 


The  two  polar  curves  intersect  at  Q: 

Pi  3  Pi 

<  — ,  - ) 

2  2 

*  After  a  breif  delay  you  will  see  the  message: 

Input  lower  limit. 

Since  the  two  curves  intersect  and  bound  two  different 
regions  you  must  be  careful  in  selecting  the  correct 
limits.  Type  pi/2. 
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You  will  see  the  message: 
Input  upper  limit. 

Type  3/2  pi. 


The  region  is  then  plotted.  You  will  see  the  following 
output : 


y 


The  outer  polar  curve  is  r  =  1  -  cos[Q] 

The  inner  polar  curve  is  r  =  1 

Pi  3  Pi 

The  limits  of  integration  are  a= —  and  b= 

2 

The  area  of  the  region  bounded  by  the  two 
polar  curves  is  2.7854 
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APPENDIX  C  “  PROGRAM  CODE 
I.  PILE  CART.MA 


(*  Copyright  Dennis  A.  Polaski,  May  4,1993  *) 

cart : = 

Module [ { } , 

cltranscend;  In [ x_] restrict ionR=n;  restrictionL=n; 
endpoints* { } ;  intersect* { } ; 

myfunc={};  rightlimit*n; leftlimit=n;  type2D*l; 


While [type2D==l,  { 

c2func=Input ["Enter  the  functions  which 

define  the  region's  boundaries.  When  finished  type  n. 


y  ="]  ; 

If[c2func  =!=n,  { 

functotal=ToExpression[StringJoin["y==" , "c2func" ] ] ; 
AppendTo[myfunc, functotal] } ,  type2D=0] ; 

I f [ Length [ my f unc ]==3 , type2  D*0 ] 

}]; 

myfunc=myfunc  /.  sqrt->Sqrt; 
myfunc=myfunc  /.  pi->Pi; 

If [Length [myfunc]==2 ,  { 

limitsknown=Input["If  any  of  the  limits  is  known  type 
y,  otherwise  type  n." ]? 

While [ limitsknown  ===y, 

leftlimit=Input[ "Input  left  limit  if  known,  otherwise 
type  n. 

x  ="  ] ; 

leftlimit=leftlimit  /.  sqrt->Sqrt; 
leftlimit=leftlimit  /.  pi->Pi; 

If[leftlimit  =!=  n,  AppendTo[myfunc, 
ToExpression[StringJoin["x==", "leftlimit"] ] ] ] ; 
rightliroit=Input[" Input  right  limit  if  known,  otherwise 
type  n. 
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X  ="]? 

rightlimit=rightlimit  /.  sqrt->Sqrt; 
rightlimit*rightlimit  /.  pi->Pi; 

If [right limit*!*  n,  AppendTo[myfunc, 
ToExpressionCStringJoinfx**", "rightlimit"]  ]  ]  ] ; 
limitsknown  =  n]  ; 

>]? 

If[leftlimit  ***  n  &&  (rightlimit  >0  |j  rightlimit===n) , 
restrictionR=Input["If  the  region  is  restricted  to  the 
right  half  plane  type  y,  otherwise  type  n."] 

]; 

If [rightlimit  ===  n  &&  restrict ionR==*n  && 

(leftlimit  <0  J  j  leftlimit===n) , 

restrictionL=Input[MIf  the  region  is  restricted  to  the 
left  half  plane  type  y,  otherwise  type  n."] 

]; 

last=Length[myfunc] ; 

For[ii=l,  ii<*  last, 

{ival*ii;  ++ii; 

For  C  j  j —1 ,  jj<*last, 

{ jval=j j ;  ++j j ; 

If [ jval!*ival  && 

(my f unc [ [ ival , 1 ] ] ***y  \\  myfunc[ [ jval, 1] ]*=*y) , 

{ 

ln[x_] :=Log[x] ; 

If [myfunc[ [ival, 1] ]==*y  &&  myfunc[ [ jval , 1] ]*==y, 

{ 

If [ {myfunc[ [ival] ] ,  myfunc[ [ jval] ] }=!= 

({myfunc[ [ival] ] ,  myfunc[ [ jval] ]  /.  Log[x_]->x) ) , 
(theroots*{nSolve[ {myfunc[ [ival] ] ,myfunc[ [jval] ] } ] } 

), 

{ 

ranscend ; 

theroots=NSolve[ {myfunc[ [ival] ] ,  myfunc[ [jval] ] } , (x,y) ] ; 
theroots=N[ {x,y}  /.  theroots]; 
ther oots=DeleteCases [ DeleteCases [ DeleteCases [ 
theroots,  (x_Complex,y_) ] , {x_,y_Complex} ] , 
{x_Complex,y_Complex} ] ; 
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I f [ Length [ theroots ] ==0 , { 

theroots=Solve [ { my f unc [ [ ival ] ] ,  my f unc [ [ j val ] ] } , { x , y } ] ; 
theroots=N [ { x , y )  / .  theroots ] ; 
theroots-DeleteCases [ DeleteCases [ DeleteCases [ 
theroots ,  { x_Complex , y_ } ] , { x_, y_Complex } ] , 

{ x_Complex , y_Complex ) ]  }  ] ; 

cltranscend ; 

If (theroots=!={x,y)  &&  theroots  !  =  {}  && 

(myfuncf [ival] ] ,myfunc[ [ jval] ] }=== 

(( (myfunc[ [ival] ] ,myfunc[ [jval]] })  /.  Sqrt[x_]->x) , 

{ 

rmtranscend ; 

theroots=Table [ FindRoot [my f unc [ [ ival ,2]]- 
myfunc[ [ jval, 2] ]==0,  {x,theroots[ [m, 133)3* 

{m, Length [theroots] } ] ; 
theroots=x  / .  theroots ; 

theroots=Table[ { theroots [ [m] ] ,  (myfunc[ [ival, 2] ]  /. 
x->  theroots[ [m] ] ) } , {m, Length [theroots] } ] ; 
temproots= 

Table [ (TrueQ[ (Abs[myfunc[ [ival, 2] ]-myfunc[ [jval, 2] ] ]  /. 
x  ->  theroots [ [ s , 1 ] ] )  <=.5] }, {s,Length[ theroots] }] ; 
theroots=Transpose[ { theroots , temproots } ] ; 

theroots=Transpose[  DeleteCases  [theroots,  {_,  (False) } ,  1]  ]  [  [1]  ] 
I f [ theroots==  {}[[!]], theroots* { }  ] ; 


x=. 

>]; 

}]; 

}]; 

rmtranscend ; 

If [myfunc[ [ival,l] ]===y  &&  myfunc[ [ jval , 1] ]===x, 
theroots=N[ { (myfunc[ [jval, 2] ] , (myfunc[ [ ival, 2] ]/. 
x->myfunc[ [jval, 2] ] ) }}]] ; 

If [myfunc[ [ival, 1] ]===x  &&  myfunc[ [ jval , 1] ]===y , 
theroots=N[ { (myfunc[ r -val, 2] ] , (myfunc[ [jval, 2] ]/. 
x->my f unc [ [ ival , 2 ]])})]] ; 

theroots=DeleteCases [ DeleteCases [ DeleteCases [ 
theroots,  (x_Complex,y_) ] , (x_,y_Complex) ] , 

{ x_Complex , y_Complex } ]  ; 

theroots=ReplaceAll [theroots,  0.  ->  0]; 

If [ restrict ionR===y , 

theroots=Select [theroots, (#[ [1] ]>=0) &] ] ; 

If [ restrictionL===y , 

theroots=Select [theroots, (#[ [1] ]<=0) &] ] ; 
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If[theroots  (x,y)  &&  theroots  =!=  (},{ 
endpoints=Join[ endpoints,  theroots]  }]; 

}]; 

If [leftlimit  =  !=n,  endpoints=Select [ endpoints, (#[ [1] ]>= 
leftlimit) &]  ] ; 

If [rightlimit  =  !=n,  endpoints=Select [ endpoints, (#[ [1] ]<= 
r ight limit) &] ] ; 

}]  ; 

endpo i nt s=N [ Round [endpoints*1000]/1000]; 
endpoints=Union[ endpoints, { } ] ; 
inter sect*Jo in [ intersect , endpoints ] ; 
intersect=Union[ intersect , { } ] ; 

ycoord=Sort[ Union [Transpose [Union [intersect, {}]][ [2] ],{}]] ; 
ycoord=N[Union[Round[ycoord*100] ]/100] ; 
mm=Transpose [ { endpoints [ [ 1 ] ] ,  endpoints [ [ 2 ] ] } ] ; 
mnew=Sort [mm[  [  1  ]  ]  ] ;  x= . ; 

If [ival==l,  { 

If  [xnnew[  [1]  ]  1=  mnew[[2]], 

{linel=Plot[myfunc[ [ival,2] ] , (x,mnew[ [1] ] ,mnew[ [2] ] } , 
Display  Functions  Identity];  linell=linel[  [1, 1,  l,  1]  ] } , 

{ 1 inel=Show[ Graphics [ (Line [endpoints] } ] , 

DisplayFunction->  Identity];  linell=endpoints} ] 

>]; 

If [ival==2,  { 

If [mnew[ [1] ]  !=  mnew[[2]], 

{line2=Plot[myfunc[ [ival, 2] ] , (x,mnew[ [1] ] ,  mnew[ [2] ] } , 
DisplayFunction->  Identity];  Iine22=line2 [ [1, 1, 1, 1] ] } , 

{ line2=Show[ Graphics [ {Line [endpoints] > ] , 

DisplayFunction->  Identity];  line22=endpoints) ] 

>]; 

If [ival==3 ,  { 

If[xnnew[[l]]  !=  mnew[[2]], 

{ linea^Plottmyfuncf [ival, 2] ] , {x,mnew[ [1] ] ,  mnew[ [2] ] } , 
DisplayFunction->  Identity];  Iine33=line3 [ [1, 1, 1, 1] ] } , 

{ 1 ine3*Show[ Graphics [ (Line [endpoints] } ] , 

DisplayFunction->  Identity];  line33=endpoints} ] 

}]; 

If [ival==4,  { 

I f [ mnew [ [ 1 ] ]  ! =  mnew [ [ 2 ] ] , 

{line4=Plot[myfunc[  [ival, 2]  ] ,  {x,mnew[  [1]  ] ,  xnnew[  [2]  ] } , 
DisplayFunction->  Identity];  Iine44=line4 [ [1, 1, 1, 1] ] } , 
{line4=Show[ Graphics [ {Line [endpoints] } ] , 

DisplayFunction->  Identity];  line44=endpoints) ] 
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>]; 

If  [ival=*=5,  { 

If [mnew[ [1] ]  !*  mnew[[2]], 

{line5=Plot[myfunc[ [ival,2] ] , {x,mnew[ [1] ] ,  mnew[ [2] ] } , 
DisplayFunction->  Identity];  line55=line[ [1, 1, 1, 1] ] } , 

{ line5=Show[ Graphics [ {Line [endpoints] } ] , 

DisplayFunction->  Identity];  line55=endpoints) ] 

)  ]  ; 

If [ival==6,  { 

If [mnew[ [1] ]  !=mnew[[2]], 

(line5=Plot[myfunc[  [ival,2]  ] ,  {x,mnew[ [1]  ] ,  mnew[  [2]  ] } , 
DisplayFunction->  Identity];  Iine66=line6[ [1, 1, 1, 1] ] } , 

{ line6=Show[ Graphics [ {Line [endpoints] > ] , 

DisplayFunction->  Identity];  line66=endpoints) ] 

}  ]  ; 

endpoints={ } 

} 

]  ? 

I f [ ival==2 ,  { Show [ 1 inel , 1 ine2 , DisplayFunct ion-> 

$Display Function] ; 
thelines={linell, line22 } } ] ; 

If [ival==3 ,  { Show [1 inel, line2 , line3 ,DisplayFunction-> 

$DisplayFunction] ; 

thelines={ linell , line22 , line33 } } ] ; 

I f [ ival==4 ,  { Show [ 1 inel , 1 ine2 , 1 ine3 , 1 ine4 , DisplayFunct ion-> 
$DisplayFunction] ; 

thelines*{ linell, Iine22,line33,line44) }] ; 

I f [ ival==5 ,  { Show [ 1 inel , 1 ine2 , 1 ine3 , 1 ine4 , 1 ine5 , 

DisplayFunct ion->  $DisplayFunction] ; 
thelines={ linell, line2 2 , line33 , line44 , lines 5 } > ] ; 

If [ ival==6 ,  { Show [ 1 inel , line2 , line3 , line4 , lines , line6 , 
Display Function->  $DisplayFunction] ; 

thelines={linell,line22,line33,line44,line55,line66) } ] ; 

intersect=Transpose [ intersect ] ; 
intersect=Sort[ intersect [ [1] ] ] ; 
intersect=ReplaceAll [ intersect , 0->0 . ] ; 

xcoord=Union [ intersect , { } ] ; 
xcoord=Union[xcoord/l. ] ; 
rmtranscend ; 

For [kk=l ,  kk<=Length[xcoord]-l, 

{ kval=kk ;  ++kk; 

xx=(xcoord[ [kval] ]+xcoord[ [kval+1] ] )/2 ; 

integrfunc=Select[thelines, (xx>  #[[1,1]]  &&  xx<  # [ [—1 ,!]])&] ; 
placel=Position[thelines,  integrfunc[ [ 1] ] ] [ [1] ] [ [1] ] ; 
place2=Position[thelines,  integrfunc[ [2] ] ] [ [1] ] [ [1] ] ; 
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x= .  ;  rmtranscend ; 

area=NIntegrate[ Evaluate [my func[ [placel, 2] ]-myfunc[ [place2,2 
]]» 

{x,  xcoord [[ kval ]] ,  xcoord [ [kval+1] ]}]] ; 
cltranscend;  ln[x_]=.; 

If[area  <0, 

{upper*myfunc[ [place2] ] ;  lower=myfunc[ [placel] ] } , 
{upper=myfunc[ [placel] ] ;  lower=myfunc[ [place2] ] } ] ; 

If [kval==l  &&  Length [ xcoord ] ==2 , 

{Print ["The  limits  of  integration  are  x  =  ", 
xcoord[ [kval] ] , "  and  x  =  ",  xcoord [ [kval+1] ]] ; 

Print [MThe  upper  curve  is  y  -  ",  upper[[2]]]; 

Print ["The  lower  curve  is  y  *  ",  lower [[2]]]; 

Print["The  area  of  the  region  is  ",  Abs[area]]  }]; 

If [kval==l  &&  Length [xcoord] ==3, 

{Print["The  limits  of  integration  for  the  first  region 
are  x  =  ",  xcoord [ [kval] ] , "  and  x  =" ,  xcoord [ [kval+1] ]] ; 
Print ["The  upper  curve  for  the  first  region  is  y  =  ", 
upper [ [ 2 ] ] ] ; 

Print ["The  lower  curve  for  the  first  region  is  y  =  ", 
lower [ [2] ] ] ; 

Print ["The  area  of  the  first  region  is  ",  Abs[area]]; 
area l=Abs[ area] } ] ; 

If  [kval==2  &&  Length [ xcoord ]~3, 

{Print ["The  limits  of  integration  for  the  second  region 
are  x  =  ",  xcoord [ [kval] ], "  and  x  *  ",  xcoord[ [kval+1] ]] ; 
Print ["The  upper  curve  for  the  second  region  is  y  =  ", 
upper [ [ 2 ] ] ] ; 

Print ["The  lower  curve  for  the  second  region  is  y  -  ", 
lower [ [2] ] ] ; 

Print ["The  area  of  the  second  region  is  ",  Abs[area]]; 

Print ["The  area  of  the  total  region  is  ",  areal+ 

Abs[area] ] } ] ; 

}]; 

inverse ; 

]  f 


cltranscend: “Module [ {}, sin[x_]=. ;  cos[x_]=.;  tan[x_]=. ; 
cot[x_]=.;  csc[x_]=.;  sec[x_]=. ;  exp[x_]=.] 

transcend : =Module [ { } , 

sin[X_]:=  X  -  XA3/6  +  XA5/120  -  XA7/5040  +  XA9/362880  - 
xA 11/3 99 168 00 ; 

COS[X_] :=  1  -  xA2/2  +  XA4/24  -  xA 6/720  +  XA8/40320  - 
XA 10/3 628800  +  XA 12/479001600  - 
XA 14/87 178 2 9 12 00  +  XA 16/20922789888000 ; 
tan[x_] :=  X  +  XA3/3  +  2xA5/15  +  17xA7/315  +  62xA9/2835  + 
1382xA 11/155925; 

cot [ X_] :=  1/x  -X/3  -xA3/45  -  2xA5/945  -  XA7/4725  - 


110 


2xA9/93555; 

CSC[X_] :=  1/X  +  X/6  +  7XA 3/360  +  31XA5/15120  + 

127xA7/604800  +  73xA9/3421440 ; 
sec[X_] :=  1  +  XA2/2  +  5xA4/24  +  61XA6/720  +  277xA8/8064  + 
5052 lxA 10/ 3 62 8 8 00 ; 

exp[x_] :=  1  +  X  +  XA2/2  +  XA3/6  +  XA4/24  +  XA5/120  +  xA6/720 
+ 

XA7/5040  +  XA8/40320  +  XA9/362880] 

rmtranscend : =Module [  { } ,  sin[x_] :=Sin[x] ;  cos[x_] :=Cos[x] ; 
tan[x_] :=Tan[x] ;  cot[x_] :»Cot[x] ?  csc[x_] :=Csc[x] ; 
sec[x_] :=Sec[x] ;  exp[x_] :=Exp[x]  ] 

nSolve [  (x==a_,y==d_:  0+c_:  l*Log[b_:  l*x__] }  ]  :  = 

{a,N[d+c*Log[b*a] ] } 

nSolve [ { y==d_ : 0+c_ : l*Log [b_: l*x_] , x==a_ ) ] : * 

{a,N[d+c*Log[b*a] ] } 

nSolve[ {y==dl_: 0+cl_: l*Log[bl_: l*x] , 
y==d2  :0+c2_:l*Log[b2_:l*x] }]:= 

N[ { ( (blAcl) / (b2Ac2) Exp[dl-d2 ] ) A (1/ (c2-cl) ) ,dl+cl*Log[bl*x]/ . 
x->( (blAcl)/(b2Ac2)Exp[dl-d2]) A (l/(c2-cl) ) >] 

nSolve [ F_List ] :  =  { x ,  y }  / .NSolve[F, {x,y} ] 


II.  FILE  CART PLOT. KA 


(*  Copyright  Dennis  A.  Polaski,  May  4,1993  *) 

cartplot:=Module[ { } ,allyfunc={ } ;allxfunc={ } ;xlines*{ } ; 

myfunc*{ } ; endpoints* { } ; intersect* { } ; 

cltranscend; 

f[{x _ }]:=x; 

type2D*l; 

Wh i 1 e [ type2 D=* 1 , 

yfunc=Input[ "Enter  the  curves,  one  at  a  time.  When 
finished  type  n. 

Y  =  "3  • 

yfunc=yfunc  /.  sqrt->Sqrt? 
yfunc=yfunc  /.  pi->Pi; 

If [yfunc=!=n,  { 

AppendTo[allyfunc,  yfunc], 

AppendTofmyfunc,  ToExpression[StringJoin[ Hy==" , "yfunc"] ] ] } , 
{ type2D=0) ] 

]; 

While [type2D==0, 

xfunc*Input ["Enter  lines  perpendicular  to  the  x-axis. 

When  finished  type  n. 


x  =  "  ]  ; 

xfunc=xfunc  /.  sqrt->Sqrt; 
xfunc=xfunc  /.  pi->Pi; 

If [xfunc=!=n,  { 

AppendTofallxfunc,  xfunc] , 

AppendTofmyfunc,  ToExpression[StringJoin[ "x=*" , "xfunc"] ] ] } , 
{type2D=-l } ] 

3 ; 

xlower=Input ["Enter  lower  limit  for  x-range  of  plot."]; 
xlower=xlower  /.  sqrt->Sqrt; 
xlower=xlower  / .  pi->Pi ; 

xupper=Input[ "Enter  upper  limit  for  x-range  of  plot."]; 
xupper=xupper  / .  sqrt->Sqrt ; 
xupper=xupper  /.  pi->Pi; 

ylower=Input[ "Enter  lower  limit  for  y-range  of  plot."]; 
ylower=ylower  /.  sqrt->Sqrt; 
ylower=ylower  /.  pi->Pi; 

yupper=Input[ "Enter  upper  limit  for  y-range  of  plot."]; 
yupper=yupper  /.  sqrt->Sqrt; 
yupper=yupper  /.  pi->Pi; 
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I f [ Length [ al lxf unc ] >*1 , { 

For[i=l,  i<=Length[allxfunc] ,  {ival=i;  ++i; 

1 ine=Show[ Graphics [ { 

Line[ { {allxfunc[ [ival] ] ,  ylower) , {ailxf unc [ [ival] ] / 

yupper}}]  }  ] ,  AxesLabel-> { "xw , "y" } , DisplayFunction 
-> Identity]; 

AppendTo [ xl ines , 1 ine ]  } ]  ; 

}]  ; 

rmtranscend ; 

theplot=Plot[Evaluate£allyfunc] , (x,  xlower,  xupper), 
AxesLabel— > { Mxw , "yM ) ,  PlotRange-> { { xlower , xupper } , 

{ylower, yupper) ) ,DisplayFunction->Identity] ; 

Show [xl ines,  theplot, DisplayFunction— >$DisplayFunction] ; 

• f y= • i 

last=Length[royfunc] ; 

For [ ii=l,  ii<=  last, 

{ ival=ii ;  ++ii; 

For [ j j=l ,  j j<=last , 

{ jval=j j ;  ++j j ; 


cl transcend; 

If [jval!=ival  && 

(my f unc [  [ ival ,  1  ]  ]  ===y  jj  myfuncf  [  jval ,  1)  ]~=y) , 
ln[x_] :=Log[x] ; 

If[myfunc[[ival,i]]=«=y  &&  myfunc[ [ jval , 1] ]===y. 

If [ (myfunc[ [ival] ] ,  myfunc[ [ jval) ] }=!= 
({myfunc[[ival]],  myfunc[ [ jval] ]  /.  Log[x_]->x)) , 
(theroots-{nSolve[ (myfunc[ [ival] ] ,myfunc[ [ jval] ] } ] } 
) » 


transcend ; 

theroots=NSolve [ {myfunc[ [ival] ] ,  myfunc[ [jval] ] ) , (x,y)  1 ; 
theroots=N [ { x , y )  /.  theroots] ; 
theroots=DeleteCases [ DeleteCases [ DeleteCases [ 
theroots,  {x_Complex,y_> ] , {x_,y_Complex) ] , 

{ x__Complex ,  y_Complex }  ] ; 

I f [ Length [ theroots ] ==o , { 

theroots=Solve[ {myfunc[ [ival] ] ,  myfunc[ [ jval] ] } , (x,y ) ] ; 
theroots=N [ { x , y }  /.  theroots]; 
theroots=DeleteCases [ DeleteCases [ DeleteCases [ 
theroots,  {x_Complex,y_} ] , {x_,y_Complex} ] , 

{x_Complex,y  Complex}] 

}]; 


cltranscend ; 


If [ theroots®! =  {x,y }  66  theroots  !={}  66 
{myfunc[ [ival] ] ,myfunc[ [jval] ] }=== 

(( {myfunc[ [ival] ] ,myfunc[ [jval] ]} )  /.  Sqrt[x_]->x) , 

{ 

rmtranscend ; 

theroots=Table[FindRoot[myfunc[ [ival, 2] ]- 
myfunc[ [jval, 2] ]==0,  {x,theroots[ [m, 1] J } ] , 

{m, Length [theroots] } ] ; 
theroot s=x  / .  theroots ; 

theroots=Table [ { theroots [ [ m ] ] , ( my f unc [ [ ival , 2 ] ]  / . 
x->  theroots [ [m] ] ) ) , { m , Length [ theroots ] } ] ; 
temproots= 

Tablet (TrueQf (Abs[myfunc[ [ ival, 2] ]-myfunc[ [jval, 2] ] ]  /. 
x  ->  theroots[ [s, 1] ] )  <=. 5] } , (s, Length [theroots] } ] ; 
theroots=Transpose[ {theroots, temproots} ] ; 

theroots=Transpose[DeleteCases[ theroots,  {False} } ,  1]  ]  [  [1]  ] ; 

I f [ theroots==  { } [ [ 1 ] ] , theroots= { }  ] ; 

x=. 

)]; 

)]; 

)]? 

rmtranscend ; 

If [myfunc[ [ival, 1] ]===y  66  myfunc[ [ jval, 1] ]===x, 
theroots=N[ { {myfunc[ [jval, 2] ] , (myfuncf [ival, 2] ]/. 
x->myfunc  [  [jval ,  2  ]  ] )  }}]],* 

If [myfunc[ [ival, 1] ]=»=x  66  myfunc[ [ jval, 1] ]===y, 
theroots=N[ { {myfunc[ [ival, 2] ] , (myfuncf [jval, 2] ]/. 
x->my f unc [ [ ival , 2 ]])}}]] ; 

theroots=DeleteCases [ DeleteCases [ DeleteCases [ 
theroots ,  { x_Complex , y_ } ] , { x_, y_Complex } ] , 

{ x_Complex , y_Complex } ]  ; 

theroots=ReplaceAll[ theroots,  0.  ->  0]; 

If [theroots  =!=  {x,y}  66  theroots  =!=  {},{ 
endpoints=Join[ endpoints,  theroots]  }]; 

}]; 

I f [ Length [ endpoints ] >0 , { 

endpoints=S elect [ endpoints , ( # [ [ 1 ] ] >=xlower ) 6 ] ; 
endpoints=Select [ endpoints , ( # [ [ 1 ] ] <=xupper ) 6 ] ; 
endpoints=Select [ endpoints , ( # [ [ 2 ] ] >=y lower ) 6 ] ; 
endpoints=Select [ endpoints , ( # [ [ 2 ] ] <=yupper) 6 ] ; 

}]; 
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)]  ; 

}  ] : 

intersect=Union[N[Round[endpoints*1000]/1000] , { } ] ; 


I f [ Length [ intersect ] == 1 , 

Print ["The  curves  intersect  at  the  point:  ", f [intersect] ]] ; 
I f [ Length [ intersect ] >=2 ,  { 

Print ["The  curves  intersect  at  the  points:"]; 

Print [ intersect ]  } ] ; 
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III.  FILE  EXAMPLES . MX 


(*  Copyright  Dennis  A.  Polaski,  May  4,1993  *) 


gr [ 1 ] : =Shov [ shades , 

Plot[  .4  Sin[x]  +  .2  , (x, .2,1.25} , 

DisplayFunction  ->  Identity, 

AxesLabel  ->  ("x","y"),  Ticks->  None, AspectRatio  ->  1, 
PlotRange  ->  { {0,1.35} , (-.07, .6} } ] , 

Graphics [ ( 

Text ["a",  {.2, -.05)], 

Text["b" ,  (1.25, -.05)] , 

Text [ "y  =  f(x)",  { .75, .6)] 

>], 

DisplayFunction  ->  $DisplayFunction] ; 

example [ 1 ] : =Show  [ 

Plot [xA2 ,  (x, 0, 1} , 

DisplayFunction  ->  Identity, 

AxesLabel  ->  {"xw/"yM)( 

Ticks  ->  {(1},{1)}, 

AspectRatio  ->  1, 

PlotRange  ->  { {-. 06, 1. 1 ) , {-. 1, 1 } } ] , 

Graphics [ { 

Line  [{(1,0),  (1,1))], 

Text["0", { .04, -.08)] , 

Text  ["y  =  x«,  (.5, .8)], 

Text  ["2", { .615, .85}] 

>], 

DisplayFunction  ->  $DisplayFunction] ; 

example [ 2 ] : =Show [ 

Plot[  Cos[x],  (x, -Pi/2,  Pi/2), 

DisplayFunction  ->  Identity, 

AxesLabel  ->  { MxH , "y" } , 

AspectRatio-> . 5 , 

Ticks  ->  None, 

PlotRange  ->  {  (-1.65,1.65),  (-.3,  1.2}}], 

Graphics [ { 

Text[  ”  y  =  cos  x" ,  (1.4,1.05)], 

Text [  "-Pi/2",  {-Pi/2,-. 2} ] , 

Text [  "Pi/2",  { Pi/2 , - . 2 } ] 

)], 

DisplayFunction  ->  $DisplayFunction] ; 
gr[2 ] :=Show[ 

Plot [ . 5* (x-1) A3  +.2  ,{x,.5,1.5),  DisplayFunction  ->  Identity, 
AxesLabel  ->  ("x","y"),  Ticks->  None,  AspectRatio  ->1, 
PlotRange  ->  { { 0, 2 } , { -. 085, . 265} } ] , 

Graphics [ ( 
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Line [ {  (.5,  0},  (.5,  ( . 5-1) A3+. 325}  }], 

Line[{  (1.5,  0),  {1.5,  .27}}], 

Line [ {  {1.68,0},  { 1. 68, .088} } ] , 

Line [ {  {1.68, .13},  {1.68, .2}}], 

Line[{  {.5, .2}, {1.5, .2}}], 

Line[{  {1.58, .2},  {1.78, .2)}], 

Line [ {  {.5,  -.08}, {.5, -.0525}}], 

Line [ {  {1.5, -.08},  { 1 . 5, -. 0525 } } ] , 

PointSize[ .01] ,  Table [ Point [{ l,n  .005}], (n,  0,  39}], 

Text ["a",  {.5,  -.025}], 

Text [ "b" ,  {1.5,  -.025}], 

Textfc",  {1,  -.025}], 

Text [ M  f(c)",  {1.68,  .11}], 

Text f 0" ,  {-.075, -.025} ] , 

Text["< -  (b  -  a)  - >",  {1.005,  -.0675}], 

Textfy  -  f  (x)",  {1,  .25}] , 

Text [ "  V  ",  {1.68, .01)3, 

Text [ "  /\  ",  {1.68, .19}] 

}]# 

DisplayFunction  ->  $DisplayFunction] ; 

example[3] : “Show [ Plot [ {Sqrt[4-xA2] , -Sqrt[4-xA2] } , {x,-2,2} , 
DisplayFunction  ->Identity, 

Ticks  ->  { { —2 ,0,2},{}}, 

Axes Label  ->  { "x" , "y" ) , 

AspectRatio  ->  Automatic, 

PlotRange  ->  { {-2 . 15, 2 . 2 } , { 0, 2 . 1} } ] , 

Graphics [{ 

Textfy  «  4  -  x  ",  {2,2.35}], 

Text["2", {2.8,2.65}] 

>], 

DisplayFunction  ->  $Display Function] ; 

example [4 ] :=Show[ Plot [2x, {x,l,4}, 

DisplayFunction  ->Identity, 

Ticks  ->  { { l, 4 } , { } } , 

AxesLabel  ->  {"x","y"}, 

AspectRatio  ->  l, 

PlotRange  ->  { { 0, 4 . 6} , { 0, 8 } } ] , 

Graphics [ { 

Textfy  =  2x" ,  {2,7}], 

Text ["(1,2) ",{.5, 2}], 

Text ["(4,8) ", {4.6,8}] , 

Line [{{1,0}, {1,2}}], 

Line [{{4,0}, {4, 8}}] 

)] , 

DisplayFunction  ->  $DisplayFunction] ; 

example [5] : =Show [ Plot [{Cos [x] ,  -Sin[x] }, {x,-. 1,1.6} , 
DisplayFunction  ->  Identity,  AxesLabel  ->  fx","y"), 

Ticks  ->  { {0} , {-1,0, 1} } ,  AspectRatio  ->  Automatic, 
PlotRange  ->  { {-. 15, 1 . 8 } , {-1. 2 , 1. 1 } } ] , 
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Graphics [ { 

Line[{ {Pi/2,0}, {Pi/2,-1}}], 

Text["Pi/2",  {1.75, .2}], 

Text["y  =  cos  x",  {1,1.1}], 

Text["y  «  -sin  x",  {.75,-1. 1}] 

)). 

DisplayFunction  ->  $DisplayFunction] ; 

example [6] : =Show[ shade 1, Plot [ l-xA (-2) , {x, .5,5}, 
DisplayFunction  ->  Identity, 

AxesLabel  ->  { "x" , "y" } ,  Ticks  ->  {{1,4},{1}}, 

AspectRatio  ->  Automatic, 

PlotRange  ->  { {-. 1, 5} , {-1.75, 1.5) } ] , 

Graphics [ { 

Line [{{0,1}, {5. 5,1}}], 

Line [{{1,0}, {1,1}}], 

Line[{{4, 15/16), {4,1}}], 

Text [ "y  =  1", {2,1.3}], 

Text [My  =  1  -  x", {1.6,-1}] , 

Text ["-2",  {2. 6, -.75}] 

}], 

DisplayFunction  ->  $DisplayFunction] ; 

example[7] :=Show[Plot[ {-x,  2-xA2) , {x,-l,2} , 
DisplayFunction  ->  Identity,  AxesLabel  ->  { "x" , "y" } , 
Ticks  ->  { {-1, 0, 1, 2 } , {-2, 0,2 } } ,  AspectRatio  ->  Automatic 
PlotRange  ->  { {-1. 95, 2 .2 } , {-2 . 4, 2 . 1} } ] , 

Graphics [ { 

Textf" (-1,1)",  {-1.65,1}], 

Text [ M (2,-2)",  {2, -2. 2}], 

Text [ "2" ,  {2.35,1.9}], 

Text [ "y  =  2  -  x  ",  {1.7, 1.7}], 

Textf'y  «  -  X",  {.9, -1.8}] 

>], 

DisplayFunction  ->  $DisplayFunction] ; 

example[8] :=Show[Plot[ {xA3+l,xA2+x} , {x, -1. 5, 1. 2 } , 
DisplayFunction  ->Identity,  AxesLabel  ->  {"x","y"}, 

Ticks  ->  {{1},{2}},  AspectRatio  ->  Automatic, 

PlotRange  ->  { { -1. 75, 1 . 5} , { -. 5, 2 . 3 } } ] , 

Graphics [ { 

Text [ "y  =  x  +  1", {-.7,1.2}] , 

Text ["3", {-.625,1.35}] , 

Text["2", {1.175, .65}] , 

Text [ "y  =  x  +  x", {1.1, .5}] , 

Text {"(-1,0)", {-1.5,-.2}] , 

Text ["(1,2) ",{1.4, 2}] 

)]# 

DisplayFunction  ->  $DisplayFunction] ; 

example[9] :=Show[Plot[ {4+xA2 , 12-xA2 } , {x,-2 .2, 2. 2} , 
DisplayFunction  ->  Identity,  Ticks  ->  None, 


Axes  Label  ->  ("x^^y"),  AspectRatio  ->  Automatic, 
PlotRange  ->  {{-4. 5, 4. 5),  {2,13}}], 

Graphics [ { 

Text["(-2,8)",{-3.2,8}], 

Text ["(2, 8)", (3.1,8)], 

Text [ "y  =  12  -  x", {-3,12.5}] , 

Text ["2", {-1.125,12.95}] , 

Textfy  -  4  +  x",  {3,3.5}], 

Text["2", (4.725,3.95}] 

>], 

DisplayFunction  ->  $DisplayFunction] ; 

example [10] :=Show[ Plot [Sqrt [x] , {x,0,4} , 

DisplayFunction  ->  Identity, 

Axes Label  ->  (wx",wyM),  Ticks  ->  { { 2 , 4 } , ( 1, 2 } } , 
AspectRatio  ->  .75, 

PlotRange  ->  { { -. 2 , 4 . 5) , { -. 2 , 2 . 2 } } ] , 

Graphics [ { 

Line [ { { 2 , 0 } , { 2 , Sqrt [ 2 ] } } ] , 

Line [{{2,0}, {4, 2}}], 

Text ["2 ",{1.7, 2. 2}], 

Text["y  =  x",  {1.85,2.1}], 

Text [" (4,2)" ,{4.4,2}], 

Text["y  =  0",  { 1 , - . 15 } ] , 

Textf’y  -  x  -  2",  {4,  .9}] 

>]» 

DisplayFunction  ->  $Display Function) ; 

example [11] :=Show[shade2 ,  Plot[2-xA2, (x,0,Sqrt[2] } , 
DisplayFunction  ->  Identity, 

AxesLabel  ->  { "x" , "y" } , 

Ticks  ->  { { 1 } , { 1 , 2 } } ,  AspectRatio  ->  Automatic , 
PlotRange  ->  { { 0, 2 ) , { -. 07 , 2 . 2 } ) ] , 

Graphics [ { 

Text ["y  =  2  -  x",  {.65,2.1}], 

Text["2", {1,2.2}] , 

Textf'y  =  x",  {1.65,1.25}], 

Text["y  =  0", { .55,-.l}] , 

Line [{{0,0}, {1.7, 1.7}}], 

Line [{{2. 1,2.1}, {2.3,1.55}}], 

Line [{{1,0}, {1,1}}] 

>]» 

DisplayFunction  ->  $DisplayFunction] ; 

example [12] :=Show[shade3,Plot[ {Sqrt[x] , -Sqrt [x] , 6-x} , 
{x,-.l,ll} , 

DisplayFunction  ->  Identity,  AspectRatio  ->  Automatic, 
AxesLabel  ->  { Mxn , "y" } ,  Ticks  ->  None, 

PlotRange  ->  { {-. 1, 11} , {-4 . 2 ,4 . 2 } } ] , 

Graphics [ { 

Line [{{0,0}, {9,0}}), 

Text[M (9,-3)", {10,-2. 2}] , 
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Text["(4,2)",{4.2,3.2}], 

Text [ "y  *  x** ,  {2. 2, -3}] 

Text ["2", {1.85,-2.67}] , 

Text [ "y  *  6  -  x", {7. 2, 1.1}] 

>3, 

DisplayFunction  ->  $DisplayFunction] ; 

gr [ 3 ] : =Show [ Plot [4(x+.l)A3,{x,.2,.5}, 

DisplayFunction  ->  Identity, 

AspectRatio  ->  .4,  PlotRange  ->  {{-.1,1. 2}, {0,1}}, 
AxesOrigin  ->  {-.1,0}, 

AxesLabel  ->  { "x",  "y"} ,  Ticks  ->  None], 

Plot [4 (x— .4)3, {x, .7,1},  DisplayFunction  ->  Identity, 
AspectRatio  ->.5,  PlotRange  ->  {{0,1}, {-.1,1. 3}}], 

Graphics [ { 

Line[{ { .485, .81} , { .99, .81} }] , 

Line[{ { .248, .175} , { .756, .175} } ] , 

Text["a", {-.15, .175}] , 

Text["b", {-.15, .81}] , 

Text [ "x  =  g(y) ", { .2, .95} ] , 

Text [ "x  =  f (y}«, {l.l, .2}] 

} ] ,  DisplayFunction  ->  $DisplayFunction] ; 

gr[4] :=Show[Plot[2.1-(x+.3) A2, {x, .59,1.1} , 

DisplayFunction  ->  Identity, 

AspectRatio  ->.4,AxesLabel->  {"x", "y" } , Ticks  ->  None, 
AxesOrigin  ->  {-.1,0}, 

PlotRange  ->{{-. 1,1. 3}, {0,1. 4}}], 

Plot[ (x+.5) A3-.2, {x, .2, .66} ,  DisplayFunction  ->  Identity, 
Ticks  ->  None, 

AspectRatio  ->.4,  PlotRange  ->  {{0,1. 3), {-.1,1. 3}}], 
Graphics [ { Line [{{.25, .2}, {1.08, .2}}], 

Text["a", {-.15, .2}], 

Text ["b", {-.15,1.235}] , 

Text["x  =  f (y) ", {1.2, .9} ] ,  Text [ "x  =  g(y) ", { .2, 1} ] } ] , 
DisplayFunction  ->  $DisplayFunction] ? 

gr[5] :=Show[Plot[ {-Sqrt[xA2-.5]+1.5,Sqrt[xA2-.5]+1.5} , {x, . 
711,1.5}, 

DisplayFunction  ->  Identity,  AspectRatio  ->.5, 

AxesLabel  ->  {"x","y"},  Ticks  ->  None, AxesOrigin  ->  {-.1,0}, 
PlotRange  ->  {{-.1,2. 55), {0,3}}], 

Plot[  {-Sqrt[3-xA2]  <-1.5,Sqrt[3-xA2]+1.5}  ,  {x,  1 . 2 , 1 . 732051 } , 
DisplayFunction  Identity,  AspectRatio  ->  .5, 

Ticks  ->  None] , Graphics [ {  Text["x  =  f (y) " , {2.2, .5}] , 

Line[ { {1.73205, 1.45} , {1.73205,1.55}}  ] , 

Text ["a", {-.25, .38} ] ,  Text [ "b" , {-. 25, 2 . 62 } ] , 

Text [ "x  =  g(y)", { .55,2.7}] }] , 

DisplayFunction  ->  $DisplayFunction] ; 
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example[13] :=Show[Plot[Sqrt[x] , {x,0,4} , 

DisplayFunction  ->  Identity, 

Axes Label  ->  { "x" , "y" } ,  Ticks  ->  { {2,4} , {1,2} ) , 

AspectRatio  ->  Automatic, 

PlotRange  ->  { {-.2,4.5} , {-.45, 2.2 )}] , 

Graphics [ { 

Line [{{2,0}, {4,2}}], 

Text [ w 2 " ,  {2.35,  2.3}], 

Text["y  -  x",  {1.8, 2.1}], 

Text ["(4,2)", {4.4,2}], 

Text [ "y  =  0",  { 1, -.3} ] , 

Text["x  »  y  +  2",  {4,. 9}] 

>]# 

DisplayFunction  ->  $DisplayFunction] ; 

example [14] :=Show[shade2 ,  Plot[2-xA2, {x,0,Sqrt[2] } , 
DisplayFunction  ->  Identity, 

AxesLabel  ->  { Mx" , "y" } , 

Ticks  ->  {{1},{1}},  AspectRatio  ->  Automatic, 

PlotRange  ->  { {0,2.25} , {-.22,2.05} }] , 

Graphics [ { 

Text ["(1,1) ",{1.3,1}], 

Text["2", {1.5, .67}] , 

Text["x  =  2  -  y",  {1.7, .6}], 

Text["y  =  x" ,  {.4, .8}], 

Text [ "y  =  0",  {.55, -.15}], 

Line [{{0,0}, {1.5, 1.5}}] 

>], 

DisplayFunction  ->  $DisplayFunction] ; 

example[l5] :=Show[shade3 , Plot [ {Sqrt[x] ,-Sqrt[x] ,6-x} , {x,-.l, 
11}, 

DisplayFunction  ->  Identity,  AspectRatio  ->  Automatic, 
AxesLabel  ->  {Mx","y"),  Ticks  ->  None, 

PlotRange  ->  { {-. 1, 11} , {-4 . 2 , 4 . 2 } } ] , 

Graphics [ { 

Line [{{0,0}, {9,0}}], 

Text ["(9,-3)”, {10,-2.2}] , 

Text ["(4,2) ",{4.2, 3. 2}], 

Text [ "x  =  y", {2. 8,-3. 2}] , 

Text ["2",  {3.875,-2.8}], 

Text["x  =  6  -  y",{7.5,.9}] 

)]» 

DisplayFunction  ->  $DisplayFunction] ; 

gr[6] :=Show[shadell,Plot[Sqrt[l-xA2] , {x, .59, .975} , 
DisplayFunction  ->  Identity, 

AxesLabel->  { "x" , "y " } , AspectRatio  ->  Automatic , 

Ticks  ->  None, PlotRange  ->  {{0,1.2}, {0, .8}}], 

Plot [Sqrt [ .5-xA2] , {x, .4375, .7071} , DisplayFunction  ->  Identity, 
Ticks  ->  None, AspectRatio  ->  1], 

Plot[Sqrt[ ,25-xA2] , {x, .45, .5} , DisplayFunction  ->  Identity, 
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Ticks  ->  None, AspectRatio  ->1] , 

Graphics [ { 

Text["r-f (Q)", (1.15, .255)3, 

Text [ "b" , { . 575 , . 5 } ]  , 

Text [ "a" , { . 55, . 125 ) ] , 

Text [ "A" , { . 9 , . 55 } ] , 

Text[MB",  { .625, .865}] , 

Line [{{0,0),{1,.5}}], 

Line[ { {0,0}, {.7, .88}}] 

}]# 

DisplayFunction  ->  $DisplayFunction] ; 
example [ 16 ] : =Show [ PolarPlot 

[2 (l+Cos[Q] ) , {Q, 0, 2Pi) ,  DisplayFunction  ->  Identity, 
AspectRatio  ->  Automatic,  AxesLabel  ->  {Hxw,Hy"}, 

Ticks  ->  None,PlotRange  ->  { { -1, 6} , {-2 . 65, 3 } } ] , 

Graphics [ { 

Line [{{0,0}, (2. 4 14,  2.414}}], 

Linet { {4,0} , {4.25, .225} }] , 

Text["2", {-.2,2.1}], 

Text [ "-2", {-.25, -2.2}] , 

Text [ "Q=0 , 2Piw ,{5.1, .35}], 

Text [ "4" ,  {4. 2, -.25}], 

Text["P(r,Q) w,  {2.6,  2.75}], 

Textf'r",  {1.2,  1.45}], 

Text["r=  2(l+cosQ)", {2.5,3.35}] 

}], 

DisplayFunction  ->  $Display Function] ; 

example [17] :=Show[ shade 12 , PolarPlot [ { l+Sin[Q] } , {Q, 0,2Pi} , 
DisplayFunction  ->  Identity,  Ticks  ->  None, 

AspectRatio  ->  Automatic, 

AxesLabel  ->  {"x","y"},  PlotRange  ->  { { -1. 8 , 2 . 5} , { -1. 1, 2 . 3 } } ] , 
Graphics [ { 

Text["r  ■  1  +  sinQ", {-.95,2.2}] , 

Text [ "Q=Pi/4 " ,  {1.65,1.65}], 

Text [ "Q=-Pi/4 " ,  {.85, -.85}], 

Line [{{0,0}, {1.4, 1.4)}], 

Line [{{0,0}, {.7, -.7}}] 

}], 

DisplayFunction  ->  $DisplayFunction] ; 

example [18] : =Show [ shade7 , PolarPlot [ { 2  Cos[Q]-H) , {Q,0,2Pi} , 
DisplayFunction  ->  Identity,  Ticks  ->None, 

AspectRatio  ->  Automatic,  AxesLabel  ->  { "x" , ”y" } , 

PlotRange->  { { -1 . 9 , 4 } , { -2 , 2 . 3 } } ] , 

Graphics [ { 

Line [{{-.75, 1.3}, {.5, -.866} }], 

Line[ { {-.75, -1.3} , { .5, .866} } ] , 

Line [{{3,0}, {3. 2, .3}}], 

Line [{{1,0}, {1.2, .3}}], 

Line [{{0,0}, {3,0}}], 
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Text["r  *  2  cosQ  +  1",  {1.75,  2.15)], 

Text["Q*2Pi/3", (-1.2,1.45)], 

Text ["Q=4 Pi/3", {-1.2, -1.45)] , 

Text["Q=Pi", {1.65, .4}] , 

Text I "Q*0", {3.55, .4)] 

)], 

DisplayFunction  ->  $DisplayFunction] ; 

example[19] :=Show[shadelO, PolarPlot [ {Sin[3Q] ), (Q,0,2Pi), 
DisplayFunction  ->  Identity,  Ticks  ->  None, 

PlotRange  ->  { { -. 9, 1. 2 ) , { -1. 05, 1. 1) ) , 

AspectRatio  ->  Automatic,  AxesLabel  ->  {"x","y")], 
Plot[Sqrt[l. a-xA2] , {x, .908,1.05) , DisplayFunction  ->  Identity, 
Ticks  ->  None,  AspectRatio  ->  Automatic], 

Graphics [ ( 

Line [{{0,0), {1.039, .6))], 

Line [{{0,0), {.6, 1.039))], 

Text["r  =sin3Q" ,  {— .6,.9)], 

Text["Q=Pi/3", { .6,1.15)] , 

Text ["Q=Pi/6",{ 1.15,. 725}], 

Text [ "Q=0" , { 1 , - . 1 ) ] , 

}]» 

DisplayFunction  ->  $DisplayFunction] ; 

example ( 20 ] :=Show[ PolarPlot [  {Sqrt[Cos[2Q] ] ) ,  {Q,0,2Pi), 
DisplayFunction  ->  Identity,  AxesLabel  ->  {"x","y"). 

Ticks  ->  None,  AspectRatio  ->  Automatic, 

PlotRange  ->  { { -1. 1, 1. 1) , {-. 8, 1. 1) ) ] , 

Graphics [ { 

Line [{{-.6, -.6), {.6, .6))], 

Line  C({  —  .6, .6), {.6,  —  .6) )], 

Text [ "r  =  cos2Q", {-.4, .95)] , 

Text ["2",  {-.6,1.05}], 

Text [ "Pi/4" , { .7, .7}] , 

Text ["-Pi/4", {.7,— .7}], 

Text["3Pi/4", {-.75, .7}] , 

Text["5Pi/4", {-.75, -.7}], 

Text ["1" , { 1. 05, -. 1) ] 

)]» 

DisplayFunction  ->  $DisplayFunction] ; 

gr[7] :=Show[shade6, Plot[Sqrt [l-xA2 ] , {x, .555, .896} , 
DisplayFunction  ->  Identity,  AxesLabel  ->  {"x","y" }, 

Ticks  ->  None,  AspectRatio  ->  Automatic], 

Plot [. 5x, {x, 0, 1. 275 }, DisplayFunction  ->  Identity, 

Ticks  ->  None], 

Plot [1. 5x, {x, 0, . 665} , DisplayFunction  ->  Identity, 

Ticks  ->None], 

Plot [Sqrt [ . 75-xA2+X] , {x, .658,1.27), 

DisplayFunction  ->  Identity,  Ticks  ->  None] , 

Graphics [ { 

Text["rl",{.65,.65}], 
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Text["r2" , { 1, .95} ] , 

Text [ "Q^a", {1.2,. 475 } ] , 

Text["Q»b",  { .45,  .95}] 

)]. 

DisplayFunction  ->  $DisplayFunction] ; 

example [21]: “Show [ shades , shade4 , 

PolarPlot[l, {Q,0,2Pi} ,AxesLabel  ->  { "x" , "y" } , 

Ticks  ->  None,  AspectRatio  ->  Automatic, 

DisplayFunction  ->  Identity], 

PolarPlotf .85, {Q, 0,1.3), 

Ticks  ->  None,  AspectRatio  ->  Automatic, 

DisplayFunction  ->  Identity] , 

PolarPlot [{ 1-Cos [Q] }, {Q, 0,2 Pi},  DisplayFunction  ->  Identity, 
Ticks  ->  None,  AspectRatio  ->  Automatic] , 

Graphics [ { 

Line [{{0,-1}, {.2, -1.3}}], 

Line [{{0,1}, {.2, 1.3}}], 

Line [{{0,0}, (.2, .98}}], 

Line [{{0,0}, {.3, .954}}], 

Line [{{0,0}, {1,0}}], 

Text [ "Q=-Pi/2" , { . 5 , -1 . 4 } ] , 

Text["Q=Pi/2", { .4,1.4}] , 

Text [ MQM , { . 6 , . 4 } ] , 

Text [ "r2  =  1", { .55,1.15}] , 

Text[''rl  *  1  -  cosQ", {-1.1, 1.5} ] 

>], 

DisplayFunction  ->  $Display Function] ; 

example  [  22  ]  :=Show[shade9,PolarPlot[{3Cos[Q]  , 
l+Cos[Q] } , {Q, 0, 2Pi } , 

DisplayFunction  ->  Identity,  AspectRatio  ->  Automatic, 
AxesLabel  ->  { "x" , "y" } ,  Ticks  ->  None, 

PlotRange  ->  { {-.3 , 3.75} , {-2 . 5, 2 .5} } ] , 

Graphics [ { 

Line[{ {0,0} , {1.25,2.16} }] , 

Line [{{0,0}, {3,0}}], 

Line [{{0,0}, {1.25,-2.16}}], 

Line [{{1. 6,-1}, {2. 2, -1.9}}], 

Line [{{1.9, 1.8}, {1.7, 1.5}}], 

Text ["Q=Pi/3", {1.25,2.3}] , 

Text ["Q=-Pi/3", {1.25, -2.3}] , 

Text[”r2  =  3COSQ", {2.7,1.95}] , 

Text["rl  «  1  +  cosQ", {3.05,-2.05}] 

>], 

DisplayFunction  ->  $DisplayFunction] ; 

example[23] :=Show[ PolarPlot [ { 1-Cos [Q] } , {Q, 0, 2Pi} , 
DisplayFunction->Identity, 

Ticks->None,  AxesLabel->{ "x" , "y" } , AspectRatio->Automatic, 
PlotRange-> { { -2 . 75 , 1 } , { -1 . 3 , 1 . 3 } } ] , 

Plot[Sqrt[ ,l-xA2] , {x,-.161, .318} , DisplayFunction->Identity , 
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Ticks->  None,AspectRatio->Automatic] , 
Graphics [ { 

Line[{ (0,0), (-.75, 1.3}}], 

Text [ "QH , { . 35, .2 } ] , 

Text["r", (-.55,  .55}] , 

Text[Hr  *  1-cosQ", {-2.2, 1.2}] , 

Text t"P(r,Q)", {-.85,1.45}], 

Text f "-2",  {-2. 15, -.15}] , 

Text["l", {-.1,1.2}] 

>], 

DisplayFunction->$DisplayFunction] ; 
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IV.  VILE  INTEGRAL. HA 


(*  Copyright  Dennis  A.  Polaski,  May  4,1993  *) 


BeginPackage [ " integral ' "  ] ; 


graph:: usage  =  "Type  graph,  to  invoke  a  program  which 
plots  a  region,  finds  the  limits  of  integration,  and 
evaluates  the  integral." 


SetOptions [ Plot ,  AxesLabel  ->  { "x" , »y» } ] ; 
SetOptions [ParametricPlot,  PlotPoints  ->  600]; 
SetOptions  [Plot,,  PlotPoints->150]  ; 

Of  f [ Drop : : drop ] ; 

Off [EndPackage: :noctx] 

Of f [FindRoot : :cvnwt] ; 

Off [General: : spell] ; 

Off [General: : spell 1] ; 

Off [General: :stop] ; 

Off [Get: : noopen] ; 

Off [Graphics: :gprim] ; 

Off [Graphics: :gptn] ; 

Off [Greater: :nord] ; 

Of  f [ Infinity : : indet ] ; 

Off [NIntegrate: :ncvb] ; 

Off [NIntegrate: :nlim] ; 

Off [NIntegrate: :slwcon] ; 

Off [ParametricPlot: :pptr] ; 

Off [Part: :partd] ; 

Off [Part: :partw] ; 

Off [Plot : :plln] ; 

Off [Plot: :plnr] ; 

Off [Power: :infy] ; 

Off [Remove: :remal] ; 

Off [ReplaceAll: :rmix] ; 

Of f [ReplaceAll : :reps] ; 

Off [Show: :gtype] ; 

Off [Solve: : ifun] ; 

Off [Solve: :tdep] ; 

Off [Unset: : norep] ; 


Needs["Graphics'Master'M] ; 
«c :  polar  l.ma; 

«c:polar2  .ma; 

«c :  cartplot .  ma ; 

«c :  inverse .  ma ; 
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«c :  examples .  ma ; 

«c:  shades,  ma; 

«c:cart.ma; 

«c :  polplotl .  ma ; 

«c :  polpl  ot  2 .  ma ; 

Begin [ " ' Private ' " ] ; 

graph : * 

Module [ { } , 

labexample= Input [ "To  display  a  graph  that  is  in  the  lab 
book  type  in  the  example  number,  otherwise  type  n.N]; 

If [labexample  =1=  n,  example [ labexample ] ,  { 

type2D=Input["Type  1  for  cartesian  coordinate  system,  or 
type  2  for  polar  coordinate  system."]; 

If [type2D==l,  { 

plotonly=Input["To  just  plot  curves  type  y,  otherwise 
type  n . " ] ; 

If [plotonly===y,  cartplot,  cart] 

>]? 


If [type2D==2,  { 

plotonly=Input [ "To  just  plot  curves  type  y,  otherwise 
type  n . " ] ; 

If [plotonly===y,  { 

typepolplot=Input["Type  1  for  plot  of  area  between  the 
origin  and  a  polar  curve,  or  type  2  for  plot  of  area 
between  two  polar  curves."]; 

If [typepolplot==l, polplotl] ; 

If [typepolplot==2,polplot2] } ,  { 

typepolar=Input [ "Type  1  for  area  between  the  origin  and 
a  polar  curve,  or  type  2  for  area  between  two  polar  curves."]; 
If [typepolar==l,  polarl]; 

If [typepolar==2 ,  polar2]; 

}]; 

>]; 

)]? 

] 

End [ ]  ; 

EndPackage [ ] ; 
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V.  FILE  INVERSE. ME 


(*  Copyright  Dennis  A.  Polaski,  May  4,1993  *) 

inverse : ^Module [ { } , 
ln[x_] :«Log[x] ; 

If [myfunc===(myfunc  /.  Log[x_]->x) , { 

transcend ; 
f[{x _ )]:-  x; 

allycoord=Transpose [Join[ f [thelines] ] ] [ [2] ] ; 
ycoord»Union[ycoord,  {Max[all vcoord] ) ,  [Minfallycoord] } ] ; 
ycoord=N [ Union [ Round [ycoord* 100] ]/100 j ; 
ycoord=Union[ReplaceAll[ycoord,  0.->0] , { )] ; 
xfuncs-Transpose [ Partition J Join [my func,  thelines] , 

Length [ my func] ] ] ; 

yfunc=Complement[xfuncs, Cases [xfuncs,  { y==b_Integer , c_ }  j 
{ y==b_Real , c_ } ] ] ; 

For [ww=l ,  ww<*Length[yfunc] , {wval=ww;  ++ww; 

If [yfunc[ [wval , 1, 1] ]===y, 

yfunc[ [wval, 1] ]=ToExpression[StringJoin[ "x==" , 

"(x/.  Solve[yfunc[ [wval, 1, 2] ]==y,  x, 

InverseFunctions->  True ] ) [ [ 1 ] ] " ] ]  ] ; 

If [yfunc[ [wval , 1] ]===True  JJ 

(yfunc[ [wval, 1,2] ]  /.  (y _ ) * (1/3) ->y)«!*yfunc[ [wval, 1, 2] ] , 

ww=99 ] 

}]; 

If [WW<99,  { 

revlines=Transpose[yfunc] [ [2] ] ; 

For [w=l ,  vv<=Length [ ycoord ] - 1 ,  (wal=w;  ++w; 
yy=(ycoord[  [wal]  ]+ycoord[  [wal+1]  ]  )/2 ; 

integrfunc=Select[revlines, ( (yy>  Min[Transpose[#]  [[2]]])  && 
(yy<  Max[Transpose[#]  [[2]]]))S]; 

If  [Length  [integr  func]  >2,  [w=99;  ww=99}  ]  ? 

>]; 

)]; 

If [ww<99,  { 

Pause [ 8 ] ; 

yintegr=Input  [  ’’Type  y  if  you  wish  to  integrate  the  same 
problem  with  respect  to  y,  otherwise  type  n."]? 

If [yintegr===y, { 

Print [ "  M ] ; 

Priirt  [m**********************************************wj ; 
invfunc=Transpose[yfunc] [ [1] ] ; 
tempinvfunc=invfunc ; 
revlines=Transpose[yfunc] [ [2] ] ; 
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For[zz=l,  zz<=Length[ycoord]-l,  {zval=zz;  ++zz; 
yy=(ycoord[ [zval] ]+ycoord[ [zval+1] ] )/2 ; 

integrfunc=Select(revlines, ( (yy>  Min [ Transpose [ #]  [[2]]])  && 
(yy<  Max[Transpose[ #]  [[2]]]))S]; 

invf unc=tempinvf unc ; 

I f [ Length [ integrf unc ] ==2 , { 
cl=Drop [Select [ integrf unc[ [1] ]  , 

(#[ [2] ]>=ycoord[ [zval] ]  &&  #[[2]]<=ycoord[[zval+l]])i],l]; 
c2=Drop [ Select [ integrfunc[ [2] ]  , 

( # [ 1 2 ] ] >=ycoord [ [ zval ] ]  &&  #[[2]]<=ycoord[[zval+l]])i],-l]; 
avgcurvel= (Transpose [cl] ( [1, 1] ]+Transpose[cl] [ [1,-1] ] )/2 ; 
avgcurve2=( Transpose [c2 ] [ [1,1] ]+Transpose[c2] [ [1,-1] ])/2; 

adjcl=Transpose [Reverse [Transpose [cl] ] ] ; 
maxcl=Max [Union [ (adjcl [ [-1, 1] ] } , { adjcl [[1,1] ]}]]; 
mincl=Min [Union [ {adjcl[ [-1,1]] },{adjcl[[l,l]])]]? 
curvel=NIntegrate[ Fit [adjcl, [l,y,yA2},y], 

{y,mincl,maxcl} ]/ (maxcl-mincl) ; 
adjc2=Transpose [Reverse! Transpose [c2] ] ] ; 
maxc2=Max [Union [ (adjc2[ [-1, 1] ] } , (adjc2[ [1,1]]}]]; 
minc2=Min [Union [ {adjc2[ [-1, 1] ] } , (adjc2[ [1,1]]}]]; 
curve2=NIntegrate[Fit[adjc2, {l,y,yA2),y], 

{y ,minc2 ,maxc2 ) ]/ (maxc2-minc2) ; 

>]; 

I f [ Length [ integrf unc ] ==2 , { 

If [avgcurve2>avgcurvel, { integrf unc=Reverse[ integrf unc] ; 
tempcurve=curve2 ; curve2=curvel ; curve l=tempcurve ; 
tempcurve=avgcurve2 ;  avgcurve2=avgcurvel ;  avgcurvel=tempcurve }  ] 
>] ; 

I f [ Length [ integrf unc ] ==2 , { 

placel=Position[revlines,  integrf ^nc[ [1] ] ] [ [1] ] [ [1] ] ; 
plare2=Position[revlines,  integrfunc[ [2]]][ri]][[l]]}, 
{theplace=Position[revlines,  integrfunc[ [1]]][[1]][[1]]}]; 

I f [ Length [ integrf unc ] ==2 , { 

If [Replace [invfunc[ [placel,2] ] , Sqrt [x_] ->0]  =1= 
invf unc [ [placel, 2] ] , 

I  f  [  curveKavgcurvel , 

invfunc[ [placel, 2] ]=negbranch [invf unc [ [placel, 2] ]  ]  ] 

] 

>]; 

If [Length [ integr f unc] ==2, { 

If [Replace [ invf unc [ [place2,2] ] , Sqrt [x_] ->0]  =!= 
invf unc [ [ place2 , 2 ] ] , 

I f [ curve2<avgcurve2 , 

invfunc[ [place2 , 2] ]=negbranch[ invfunc[ [place2 , 2] ]  ]  ] 

] 

)]? 
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y*. ; 

If [ Length [integrfunc] ==2,  { 

area=NIntegrate[invfunc[ [placel, 2] ]-invfunc[ [place2, 2] ] , 

(y,  ycoord[ [zval] ] ,ycoord[ [zval+1] ]}]},{ 

area=NIntegrate [ Evaluate [2bothbranches[invfunc[ [theplace,2] ] 
}}' 

ty ,ycoord[ [zval] ] ,  ycoord[ [zval+l] ]}]}]; 

If [Length [ integrfunc] ==1,  { 

right=invfunc[ [theplace,2] ] ; 

lef t=negbranch [ invfunc [ [ theplace , 2 ] ]  ] } , 

{right=invfunc[ [placel, 2] ] ; left=invfunc[ [place2 , 2] ] } ] ; 

Print ["  "]; 

If [ zval==l  &&  Length [ ycoord ] ==2 ,  { 

Print ["The  limits  of  integration  are  y  =  ", 
ycoord[ [zval] ] , "  and  y  =  ",  ycocrd[ [ zval+1] ]  ]; 

Printf'The  right  curve  is  x  =  ",  right]; 

Print["The  left  curve  is  x  =  ",  left]; 

Print["The  area  of  the  region  is  ",  Abs[area]  ]  }]; 

If [zval==l  &&  Length [ycoord ] >=3 ,  { 

Print["The  limits  of  integration  for  the  first  region"]; 
Print ["are  y  =  " , ycoord [ [zval] ] , 

"  and  y  =  ", ycoord [ [zval+1] ]] ; 

Print ["The  right  curve  for  the  first  region  is  x  *  ",  right]; 
Print["The  left  curve  for  the  first  region  is  x  *  ",  left]; 
Print ["The  area  of  the  first  region  is  ",  Abs[area]  J; 
areal=Abs [ area ]  } ] ; 

If [ zval==2  &&  Length [ ycoord ]>=3,  { 

Print["The  limits  of  integration  for  the  second  region"]; 
Print ["are  y  =  ", ycoord [ [zval] ] , 

"  and  y  =  ", ycoord [ [zval+1] ]] ; 

Print ["The  right  curve  for  the  second  region  is  x  =  ",  right]; 
Print ["The  left  curve  for  the  second  region  is  x  =  ",  left]; 
Print ["The  area  of  the  second  region  is  ",  Abs[area]  ]; 

I f [ Length [ ycoord ]== 3 , 

Print ["The  area  of  the  total  region  is  ",  areal+Abs[area]  ], 
area2=Abs[area]  ] 

}]? 

If[zval==3  &&  Length [ ycoord ] >=4 ,  { 

Print["The  limits  of  integration  for  the  third  region"]; 
Print ["are  y  =  ", ycoord [ [zval] ] , 

"  and  y  =  " .ycoord [ [zval+1] ]] ; 

Print["The  right  curve  for  the  third  region  is  x  =  ",  right]; 
Print["The  left  curve  for  the  third  region  is  x  =  ",  left]; 
Print["The  area  of  the  third  region  is  ",  Absfarea]  ]; 

I f [ Length [ ycoord ] ==4 , 

Print ["The  area  of  the  total  region  is  ",  areal+area2+ 

Abs [ area ] ] ,  area3=Abs [ area ]  ] 
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If[zval==4  &&  Length [ycoord]*=5,  { 

Print ["The  limits  of  integration  for  the  fourth  region”]; 
Print ["are  y  -  ”,  ycoord[ [zval] J , 

"  and  y  **  " ,ycoord[ [zval+1] ] ] ; 

Print ["The  right  curve  for  the  fourth  region  is  x  =  right]; 
Print ["The  left  curve  for  the  fourth  region  is  x  =  ",  left]; 
Print ("The  area  of  the  fourth  region  is  ",  Abs[area]  ]; 

Print ("The  area  of  the  total  region  is  ",  areal+area2+area3+ 
Abs[area] ] 

}]; 

}]; 

}]; 

}]; 

)]; 

] 

gl[(b_  +  Sqrt[a_] ) *c_] :=(b  -  Sqrt[a])*c; 
g2[(b_  +  d_*Sqrt[a_] ) *c_] :=(b  -d  Sqrt[a])*c; 
g3[  b_  +  d_*Sqrt[a_]  ]:=  b  -d  Sqrt [a]; 
g4[  Sqrt[a_]  *c_]:=  -Sqrt[a]  *c; 

g5[  Sqrt[a_]  ]:=  -Sqrt[a] ; 

g6[  b_  +  Sqrt[a_]  ]:=  b  -Sqrt[a]; 

hl[(b_+  Sqrt[a_] ) *c_] :=  Sqrt[a]*c; 

h2((b_  +  d_*Sqrt(a_] ) *c_] ;»d*Sqrt[a]*c; 
h3[  b_  +  d_*Sqrt[a_]  ] ;«d*Sqrt [a] ; 
h4[  Sqrt[a_]  *c_]:=  Sqrt[a]*c; 

h5[  Sqrt [a_]  ]:=  Sqrtfa] ; 

h6[  b_  +  Sqrt[a_]  ];=  Sqrt [a]; 

negbranch[x_] := 

If [Cases[ {x} , (b_  +  Sqrt[a_] ) *c_]=!={ } ,gl[x] , 

If [Cases[ {x} , (b_  +  d_*Sqrt[a_] ) *c_]=!={ } ,g2 [x] , 

If [Cases[ {x} ,  b_  +  d_*Sqrt[a_]  ]=!={} ,g3 [x] , 

If [Casesf (x) ,  Sqrt[a_]  *c_]=!={ } ,g4 [x] , 

If [Cases[ {x} ,  Sqrt [a_]  ]=!={} ,g5[x] , 

If [Cases [ (x) ,  b_  +  Sqrt[a_]  ]=!=(},g6[x]]]]]]]; 

bothbranches[x_] ;= 

If [Cases [ {x} , (b_  +  Sqrt[a_] ) *c_]=!={ } ,hl[x] , 

If [Cases[ {x} , (b_  +  d_*Sqrt [a_] ) *c_]=!={ } ,h2 [x] , 

If [Cases[ (x> ,  b_  +  d_*Sqrt[a_]  ]=!={} ,h3 [x] , 

If [Cases[ (x) ,  Sqrt[a_]  *c_]=!={ } ,h4 [x] , 

If [Cases[ {x} ,  Sqrt[a_]  ]=!={} ,h5[x] , 

If [Cases [ {x} ,  b_  +  Sqrt[a_]  ]=!={} ,h6[x] ]]]]]] ; 


▼I.  FILE  PO LARI. HA 


(*  Copyright  Dennis  A.  Polaski,  May  4,1993  *) 

polarl:= 

Module [ ( 

func,  limitsknown,  theroots,  lowerlimit, 
upperlimit,  newcurve,  area,  totalarea={ } ,  regions={), 
curvesknown,  fund,  func2,  temp 

>, 

Q=.  ; 

cltranscend; 

func=Input[ "Enter  the  polar  curve  which  bounds  the  region. 


r  =  "  ]  ; 

func=func  /.  sqrt->Sqrt; 
func=func  /.  pi->Pi; 
func=func  /.  q->Q; 

revfunc=Replace[ func,  Sqrt[x_]  ->  x] ; 
limitsknown=Input ["If  the  limits  are  known  type  y, 
otherwise  type  n."]; 

I f [ 1 imitsknown===y , 

{ lower limit=Input[" Input  lower  limit."]; 
lowerlimit=lowerlimit  /.  sqrt->Sqrt; 
lowerlimit=lowerlimit  /.  pi->Pi; 

While [lowerlimit  <-N[Pi],  {lowerlimit=Input["The  lowerlimit 
can  not  be  less  than  -Pi,  input  lower  limit  again."]; 
lowerlimit=lower limit  /.  sqrt->Sqrt; 
lowerlimit=lower limit  /.  pi->Pi; 

)]; 

upperlimit=Input[" Input  upper  limit."]; 
upperlimit=upperlimit  /.  sqrt->Sqrt; 
upperlimit=upperlimit  /.  pi->Pi; 

While [upperlimit>2  N[Pi],  (upperlimit=Input["The  upperlimit 
can  not  be  greater  than  2Pi,  input  upper  limit  again."]; 
upperlimit=upper limit  /.  sqrt->Sqrt; 
upperlimit=upper limit  /.  pi->Pi; 

>]; 

If [N [upperlimit ]<N[ lowerlimit] , lowerlimit=N[lowerlimit-2Pi] ] 
{lowerlimit=0. ;  upperlimit*  2Pi}]; 
rmtranscend ; 

max=Max[DeleteCases[ Table [func  /.  Q  ->  values, 

{values,  0,  2Pi,  . 25} ] ,x_Complex] ] ; 

rayl=N[(func  /.  Q->lowerlimit) ] ; 
ray2=N[(func  /.  Q->upperlimit) ] ; 
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llequivalent*N (Mod [lower limit,  2Pi]] ; 
ulequi valent =N[ Mod [upper limit,  2Pi] ] ; 

If [lowerlimit  <  0,  llequivalent*  2  N[Pi]+lowerlimit, 
llequivalent*  lowerlimit] ; 

If [upperlimit  <  0,  ulequivalent=  2  N [ Pi ]+upper limit, 
ulequivalent*  upperlimit]; 

Q=.  ; 

If [ (rayl*=ray2)  &&  (llequivalent==ulequivalent) j | 

(rayl==0  &&  ray2**0) , { 

PolarPlot [ f unc ,  { Q ,  topi [lowerlimit] ,  topi [upperlimit] } , 
AxesLabel  ->  {"x","yM},  AspectRatio  ->  Automatic, 

PlotRange->  {{-1.02  max,  1.02  max), {-1.02  max,  1.02  max}}]}, 
{ 

Show[ PolarPlot [ func, {Q,  topi [lowerlimit] ,  topi [upperlimit] } , 
DisplayFunction  ->  Identity,  AxesLabel  ->  { "xM , My" } , 
AspectRatio  ->  Automatic,  PlotRange->  {{-1.02  max, 1.02  max), 
{-1.02  max,  1.02  max}}], 

Graphics [ { 

Line[ {{0,0}, {Cos [lowerlimit] *rayl.  Sin [lowerlimit] *rayl) } ] , 
Line[ { {0,0} , {Cos [upper limit]* ray 2,  Sin [upperlimit] *ray2 } } ] , 

)3» 

DisplayFunction->  $DisplayFunction] } ] ; 
ptranscend ; 

theroots*NSolve [ { r==0 ,  r== rev func } , { r , Q } ] ; 
theroots=N[ {r,Q}  /.  theroots] ; 

If [theroots=!={r,Q}  &&  theroots  !={},{ 
theroots=DeleteCases [theroots , { r_, Q_Complex  / ; 

Abs[Im[Q] ]>. 1} ] ; 

theroots=ReplaceAl 1 [ theroots , { r_ , Q_ }  ->  { r , Re [ Q ] } ] ; 
rmtranscend ; 

Q=. ; 

theroots=Table [ FindRoot [ revf unc==0 , 

{Q,theroots[ [m, 2] ] } ] , {m, Length [theroots] } ] ; 
theroots* { 0 , Q }  /.  theroots; 
theroots=N[ Union [theroots, { } ] ] ; 

theroots=Union [theroots,  theroots+  N [Table [{ 0, 2Pi } , 

{ i , Length [ theroots ]}]]]; 

theroots=Select [theroots, (#[ [2] ]<=N[ upperlimit]  && 

#  [  [ 2 ] ]  >=  N[ lowerlimit] )&] ; 

theroots=Union[N[ (Round [100*theroots] )/100] , { }] ; 

)]? 

rmtranscend ; 

If [limitsknown=!=y,  { 

If [theroots* I ={r,Q}  &&  theroots!={}  && 

Length [ theroots ]>=1, { 

Print [ "  " ]  ; 
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Print ["The  polar  curve  r  is  equal  to  zero  when  Q  is: 

* , Hap [ topi , Transpose [ theroots ] [ [ 2 ] ] ]  ] 

>].* 

Pause [ 8 ] ; 

newcurve=Input["The  polar  curve  is  plotted  from  Q=o 

to  Q—2  Pi.  If  you  want  to  change  the  limits  type  y,  otherwise 

type  n . " ] ; 

I f [ newcurve===y , { 

lowerl imit=Input[" Input  new  lower  limit."]; 
lower limit=lower limit  /.  sqrt->Sqrt; 
lowerl imit=lower limit  /.  pi->Pi; 

While [lower limit  <-N[Pi],  {lowerlimit=Input["The  lower  limit 
can  not  be  less  than  -Pi,  input  lower  limit  again."]; 
lowerl imit=lowerl imit  /.  sqrt->Sqrt; 
lowerlimit=lowerlimit  /.  pi->Pi; 

>]; 

upperlimit=Input[" Input  new  upper  limit."]; 
upperlimit=upperlimit  /.  sqrt->Sqrt; 
upperlimit=upperlimit  /.  pi->Pi; 

While [upperlimit>2  N[Pi],  {upperlimit=Input ["The  upperlimit 
can  not  be  greater  than  2Pi,  input  upper  limit  again."]; 
upper limit=upper limit  /.  sqrt->Sqrt; 
upperlimit=upperlimit  /.  pi->Pi; 

}]; 

If  [N [upperlimit ]<N[lowerlimit] ,  lowerlimit=N[lowerlimit-2Pi]  ] ; 

rayl=N[ (func  /.  Q->lowerlimit) ] ; 
ray2=N[(func  /.  Q->upperlimit) ] ; 

llequivalent=N [Mod [lower limit,  2Pi]] ; 
ulequivalent=N[ Mod [upperlimit,  2Pi] ] ; 

If [lowerlimit  <  0,  llequivalent=  2  N[Pi]+lowerlimit, 
llequivalent=  lowerlimit] ; 

If [upperlimit  <  0,  ulequivalent=  2  N [ Pi ]+upper limit, 
ulequivalent=  upperlimit]; 

Q=  * ; 

If [ (rayl==ray2)  &&  (llequivalent==ulequivalent) J { 

(rayl==0  &&  ray2==0) , { 

PolarPlot [ func , { Q ,  top i[ lower 1 imit ] ,  topi [upperlimit] } , 
AxesLabel  ->  {"x","y"},  AspectRatio  ->  Automatic, 

PlotRange->  {{-1.02  max,  1.02  max}, {-1.02  max,  1.02  max}}]}, 
{ 

Show [PolarPlot [func, {Q,  topi [lowerl imit ] ,  topi [upperlimit] } , 
DisplayFunction  ->  Identity,  AxesLabel  ->  {"x","y"}, 
AspectRatio  ->  Automatic,  PlotRange->  {{-1.02  max, 1.02  max), 
{-1.02  max,  1.02  max}}], 

Graphics [ { 

Line[ { {0,0} , { Cos [lowerl imit ] *rayl,  Sin [lowerl imit ]*rayl) } ] , 
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Line[ { {0,0} , (Cos [upperl imit ]*ray2,  Sin [upper limit] *ray2 } } ] , 


H. 

DisplayFunction->  $Display Function] } ] 

)  3 
>]? 

ptranscend ; 

theroots=NSolve [ (r==0,  r==revfunc), {r,Q}] ; 
theroots=*N[  (r,Q)  /.  theroots] ; 

If [theroots*! ={r,Q}  &&  theroots  !={},{ 
theroots=DeleteCases [ theroots, (r_,Q_Complex  /; 

Abs[Im[Q] ]>. 1 } ] ; 

theroots=ReplaceAll [theroots, (r_,Q_)  ->  (r,Re[Q]}]; 
rmtranscend ; 

Q-  •  • 

theroots=Table [ FindRoot [ revfunc==0 , 

{Q,theroots[ [m,2] ] } ] , {m, Length [theroots] } ] ; 
theroots={0,Q)  /.  theroots; 
theroots«N [Union [theroots, { } ] ] ; 

theroots=Union[ theroots,  theroots+  N[Table[ { 0, 2Pi } , 

{ i , Length [ theroots ])]]]; 

theroots=Select [theroots, (#[ [2] ]<=N[upperlimit]  && 

#  [  [2 ]  ]  >=  N [ lower 1 imit] )&] ; 

theroots=Union[N[ (Round [100*theroots] )/100] , { } ] ; 

}]? 

I f [ theroots  === ( r , Q  >  J  j  Length [ theroots ]<=1,  ( 
area=NIntegrate[ .5*(func) A2,  {Q, topi [lower limit] , 
topi [upperl imit] } ] ; 
cltranscend; 

Print [ "  " ] ; 

Print ["The  limits  of  integration  for  the  region  bounded  by 
theH] ; 

Print [ "polar  curve  r  =  " ,func,  ",  and  the  origin  are:"]; 
Print [ "  " ] ; 

Print ["  a=", topi [Chop [lower limit, .01] ] ,"  and  b=", 

topi [Chop [upperl imit, .01] ],"."]; 

Print ["  "]; 

Print ["The  area  bounded  by  the  polar  curve  and  the  origin  is 

t« 

N[area, 3] ] } ] ; 

If [theroots  =!={r,Q)  &&  theroots!={)  &&  Length [ theroots ]>1,  { 

theroots=Transpose[ theroots] ; 
theroots=theroots [ [ 2 ] J ; 
theroots=Sort [Union [ReplaceAll [Join[ 

{ N [ (Round[100*lowerlimit] )/100] } , 

theroots, {N[ (Round [ 100*upper limit ] )/100] } ] ,0.->0] ,{}]]; 
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rmtranscend  ; 

For[i=l,  i<=Length [ theroots ] -l,  {ival=i;  ++i;  area=. ; 

Q= (theroots [ [ival] ]  +  theroots[ [ival+lj ] )/2 ; 

If [func=!=revfunc,  { 

If [N[func]  >0,  {  Q». ; 

area=NIntegrate[ . 5* (func) A2 , { Q , topi [ theroots [ [ival] ]  ] , 
topi[theroots[ [ival+1] ]  ]}]; 

AppendTo[ total area,  Abs[area] ] ; 

AppendTo [ regions ,  ival ]  } ] 

>]; 

I f [ f unc===revf unc ,  { 

If [Abs[N[func] ]  >0,  {  Q=.; 

area=NIntegrate[ .5* (func) A2, {Q, topi [theroots [ [ival] ]  ] , 
topi [ theroots [ [ival+1] ]  ]}]; 

AppendTo [ totalarea ,  Abs [ area ] ] ; 

AppendTo [regions,  ival]  )] 

)] 

}]; 

area=Apply[Plus,  totalarea]; 
munregions=Length[ regions] ; 
cltranscend; 

If [numregions==l,  {Q=. ; 

Print [ "  " ] ; 

Print ["The  limits  of  integration  for  the  region  bounded  by 
the"] ; 

Print ["polar  curve  r  *  ",  func,"  and  the  origin  are;"]; 

Print ["  "]; 

Print ["  a=", topi [Chop [theroots [ [  regions[[l]] 

]],.01]], 

"  and  b=", topi [ Chop [ theroots [ [regions [ [1] ]  +  1] ] , . 01] ] , " . "  ]; 
Print [ "  " ] ; 

Print ["The  area  bounded  by  the  polar  curve  and  the  origin 
is  ",N[area,3]  ]}]; 

If [numregions  >1,  { 

Print [ "  " ] ; 

Print ["The  limits  of  integration  for  the  regions  which 
are  bounded"];  Q=. ;  cltranscend; 

Print ["by  the  polar  curve  r  =  " ,  func,"  and  the 
origin  are;"];  Print["  "]; 

If [Abs[N[Mod[theroots[ [regions[ [1] ]  ]]+2Pi,2Pi]]  - 

N[ Mod [ theroots [ [regions! [numregions] ]+l] ]+2Pi, 2Pi] ] ]<. 1, { 

If [theroots [ [regions [ [l] ]  ]]<0  && 
theroots [ [ regions [ [ numregions ] ]  ] ] >=N [ Pi ] , 

Print ["  al=", 

topi [Chop[N[theroots[  [regions!  [numregions]  ] 
] ]-2Pi] , .01] ] , "bl=", 

topi [ Chop [ theroots [ [regions! [1] ]+l] ] , .01] ]  ] ] ; 
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If [theroots[ [regions! [1] ]  ]]<0  && 
theroots [ [ regions [ [ numregions ] ]  ] ] <-N [ Pi ] , 

Print ["  al-", 

topi [ Chop [ theroots [ [regions [ [numregions] ]  ] ] , .01] ] , 

"  and  bl-", topi [Chop [N[ theroots [ [regions! [1] ]+l] ]+2Pi] , .01] ] 

]]; 

If [theroots! [regions [ [1] ]  ]]>*0  && 
theroots [ [ regions [ [ numregions ] ]  ]  ]  >0 , 

Print ["  al=", 

topi [ Chop [N[ theroots [ [ regions [ [numregions] ]  ] ]-2Pi] , .01] ] , 

"  and  bl®1*,  topi  [  Chop  [  theroots  [  [regions!  [1]  ]+l]  ] ,  .01]  ]  ]]; 

For[ j=3 , j<=numregions,  {jval=j;  ++ j ; 

Print ["  a", jval-1, 

topi [ Chop [ theroots [ [regions[ [ jval-l] ]  ] ] , .01] ] , 

M  and  b", jval-l, M=w, 

topi [ Chop [ theroots [ [ regions [ [jval-l] ]+l] ] , .01] ]  ] } ] 

}/ 

{For[j=l,  j<=numregions, { jval=j ;  ++ j ; 

Print ["  a",  jval,M=", 

topi [ Chop [ theroots [ [  regions [ [jval] ]  ]],.01]], 

"  and  b", jval, , 

top i[ Chop [ theroots [ [  regions [ [jval] ]+l  ] ] , . 01] ]  ]  }]; 

>]; 

Print ["  "]?  Q=. ; 

Print [MThe  total  area  bounded  by  the  polar  curve  and  the 
origin  isM] ; 

Print [ "  " , N [ area , 3 ] ] ; 

>]; 

)]? 

] 


ptranscend:=Module[ { } , 

sin[X_]:=  X  -  XA3/6  +  XA5/120  -  XA7/5040  +  XA9/362880  - 

XAll/39916800-t-XA13/6227020800-XA  15/1307674368000+ 
XA17/355687428096000-XA 19/121645100408832000 ; 

XA 2 1/5109 094 2 17 17 094 4 00 00  - 
XA 2 3/2 5852 016738884 97664 0000  + 

XA 2 5/ 1551 12 1004 3 33 09 8 59 8 40 000 00 ; 
cos[x_]:=  1  -  xA2/2  +  XA4/24  -  xA6/720  +  xA8/40320  - 
XA10/3628800  +  XA 12/479001600; 
tan[X_] :=  x  +  XA3/3  +  2xA5/15  +  17xA7/315  +  62xA9/2835  + 
1382 XA1 1/155925; 

COt [X_] :=  1/X  -x/3  -XA3/45  -  2xA5/945  -  XA7/4725  - 
2xA9/93555; 

csc[x_]:=  1/x  +  x/6  +  7xA  3/360  +  31XA5/15120  + 

127xA7/604800  +  73xA9/3421440 ; 
sec[X_] :=  1  +  xA2/2  +  5xA4/24  +  61XA6/720  +  277xA8/8064  + 
5052 lXA 10/3 62 88 00 ; 
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exp[x_]:=  1  +  X  +XA2/2  +  XA3/6  +  XA4/24  +  XA5/120  +  XA6/720  + 
XA7/5040  +  XA8/40320  +  XA9/362880] 

cltranscend:=Module[ { } , sin[x_]=. ;  cos[x_]=.;  tan[x_]=.? 
cot[x_]*.;  csc[x_]*.;  sec[x_]=. ;  exp[x_]*.] 

rmtranscend:=Module[ { } ,  sin[x_] :=Sin[x] ?  cos[x_] :=Cos[x] ; 
tan[x_] :=Tan[x] ;  cot[x_] :=Cot[x] ;  csc[x_] :=Csc[x] ; 
sec[x_] :=Sec[x] ;  exp[x_] :=Exp[x]  ] 

topi[-.26] :=-Pi/12;  topi[-. 28] :=-Pi/li;  topi[-. 29] :*-Pi/ll; 
topi[-. 31]  :=-Pi/10;  topi[-.34] Pi/9;  topi[-.35]:=-Pi/9; 
topic-. 39]  :=-Pi/8;  topi [-. 44 ];— Pi/7?  topi [-. 45]  :=-Pi/7 ; 
topi[-. 52] :*-Pi/6;  topi[-.57] :=-2Pi/ll;topi[-.62] :=-Pi/5; 
topic-. 63] :®-Pi/5;  topi [-.69] :=-2Pi/9;  topi [-.70] :=-2Pi/9; 
topi [-.78] :=-Pi/4;  topi [-.79] :=-Pi/4 ;  topi[-. 85] :=-3Pi/ll ; 
topic-. 86]  :=-3Pi/ll?  topi[-. 94]  :=-3Pi/10;  topi[-l.04] Pi/3? 
topi [-1.05] :«-Pi/3;  topi[-l. 14] :=-4Pi/ll;  topi[-l.l7]:=-3Pi/8; 
topi[-l. 18] :=-3Pi/8;  topi [-1.25] :=-2Pi/5;  topi[-l. 26] :=-2Pi/5; 
topi [-1.30]  :=-5Pi/12?  topi  [-1.  31]  :=-5Pi/12  ? 
topi [-1.34] :®-3Pi/7 ;  topi[-1.35] :=-3Pi/7;  topi[-1.39] :=-4Pi/9; 
topi [ - 1 . 4  o ]  s ■ - 4  P i / 9  ?  topi[-1.42]  :»-5Pi/ll? 

topi [-1.43]  s  »-5Pi/ll  ?  topi [ -1  .  57  ]  :=-Pi/2; 
topi [ -1 .71]  : *-6Pi/ll ;  topi [ -1  .  74  ]  : =-5Pi/9 ; 
topi[-1.75] :=-5Pi/9;topi[-1.79] :=-4Pi/7 ;topi[-1.80] :=-4Pi/7; 
topi [-1.83] :=-7Pi/12; 

topi  [-1.88]  :  ■  -  3  P  i  /  5  ;  topi  [-1  .  96]  : -  5Pi/8  ? 
topi [-1.99] s*-7Pi/ll?  topi [ -2 . 00 ] :*-7Pi/ll? 
topi [-2. 09]  : =-2  Pi/ 3  ;  topi [-2.  19]  :=  -7  Pi/ 10? 

topi  [-2.20]  :  *  -  7  P  i  /  1  0  ;  topi  [-2  .  24  ]  :-5Pi/7  ? 

topi  [-2.28]  :=-8Pi/ll;  topi  [  -2  .  3 5 ]  :  =-3Pi/4  ; 

topi [-2.36]: — -3Pi/4 ;  topi [-2.44]: =-7 Pi/ 9 ;  topi [-2.51]: =-4Pi/ 5  J 
top i [ -  2 . 5  7 ]  :=-9Pi/ll?  topi[-2.61]  :  =  -  5  P  i  /  6  ; 

topi [-2.62]: =-5Pi/6 ;  topi [ -2.69]: =-6Pi/7 ;  topi [-2.74]: =-7Pi/8 ; 
topi  [-2.75]  : * -  7  P i / 8  ;  topi [ -2  .  79  ]  : =-8Pi/9  ; 
topi  [-2.82]  :=-9Pi/10;  topi  [-2.  83  ]  :-9Pi/10? 
topi [-2 .85] :=-10Pi/ll;  topi[-2.86]:— 10Pi/ll; 

topi  [-2.87]  -llPi/12;  topi  [-2. 88]  :=-HPi/12; 

topi [-3.14] : =-Pi ;  topi [ . 26] :=Pi/12 ;  topi[ .28] :=Pi/ll? 
topi [ . 29] :=Pi/ll ;  topi[ . 31] :=Pi/10;  topi [ . 34 ] :=Pi/9 ; 
topi [ . 35] : =Pi/9 ;  topi [ . 39 ] : =Pi/8 ;  topi [ . 44 ] : =Pi/7 ; 
topi[  *45]  :=Pi/7;topi[ .  52]  :=Pi/6;topi[  .57]  :=2Pi/ll;topi[  .62]  :  = 
Pi/5;  topi[ .63] :=Pi/5;  topi[ .69] :=2Pi/9;  topi[ . 70] :=2Pi/9 ; 
topi[ .78] :=Pi/4;  topi[ .79] :=Pi/4 ;  topi[ .85] :=3Pi/ll; 
topi[ .86] :=3Pi/ll;  topi[ .94] :=3Pi/10;  topi[1.04] :=Pi/3 ; 
topi [1.05] :=Pi/3 ;  topi[l. 14] :=4Pi/ll;  topi[l. 17] :=3Pi/8; 
topi [l. 18] :=3Pi/8 ;  topi [ 1 . 25] :=2Pi/5 ;  topi[1.26] :=2Pi/5; 
topi[1.30] :=5Pi/12;  topi[1.31] :=5Pi/12;  topi[1.34] :=3Pi/7; 
topi[1.35] :=3Pi/7;  topi[1.39] :=4Pi/9;  topi[1.40] :=4Pi/9; 
topi[1.42] :=5Pi/ll;  topi[1.43] :=5Pi/ll;  topi[l. 57] :=Pi/2 ; 
topi [1.71] :=6Pi/ll;  topi [ 1 . 74 ] :=5Pi/9 ;  topi[1.75]:=5Pi/9; 
topi [1.79] :=4Pi/7;  topi [1.80] :=4Pi/7;  topi [1.83] :=7Pi/12; 
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topi [1.88] :*3Pi/5;  topi [1.96] :«5Pi/8;  topi[1.99]  :=7Pi/U; 
topi [2. 00] :«7Pi/ll;  topi[2.09] :-2Pi/3;  topi[2.19]:»  7Pi/lO; 
topi [2. 20] :=7Pi/10;  topi[2.24] :»5Pi/7;  topi[2.28] :-8Pi/ll; 
topi [2 .35] :=3Pi/4 ;  topi[2 . 36] :*3Pi/4 ;  topi[2.44] :=7Pi/9; 
topi [2 . 51] :*4Pi/5;  topi[2.57] :«9Pi/ll;  topi[2.61] :=5Pi/6; 
topi [2. 62] :»5Pi/6;  topi [2 . 69] : »6Pi/7 ;  topi [2. 74] :=7Pi/8; 
topi [2 . 75] :=7Pi/8 ;  topi[2.79] :=8Pi/9;  topi[2.82] :*9Pi/10; 
topi [2. 83] :=9Pi/10;  topi[2.85]  :*10Pi/ll,*topi[2.86] :=10Pi/ll; 
topi [2. 87] :*llPi/12;  topi [2. 88] :®llPi/12?  topi[3. 14] :=Pi; 
topi [3.40] :*13Pi/12;topi[3.42]:=12Pi/ll;topi[3.43]:«12Pi/ll; 
topi [3. 45]  :*llPi/10;  topi [3. 46]  :«HPi/10?topi[3.49]  :«10Pi/9; 
topi [3. 59] :-8Pi/7;  topi [3. 53] :*9Pi/8?  topi [3. 66] :»7Pi/6; 
topi [3.67] : =7Pi/6 #  topi[3.7l] :«l3Pi/ll;  topi[3.76]:-6Pi/5; 
topi [3. 77]  :=6Pi/5;  topi [3. 83]  :*HPi/9;  topi[3.84]:*HPi/9; 
topi [3 .92] :=5Pi/4 ;  topi [3. 93] :=5Pi/4 ;  topi [3. 99] :=14Pi/ll; 
topi [4. 00] :=14Pi/ll;  topi [4. 03] :=9Pi/7;  topi [4. 04] :«9Pi/7; 
topi [4.08] :=13Pi/10 ;  topi[4. 18] :=4Pi/3;  topi[4 . 19] :=4Pi/3 ; 
topi [4.28] :  =15Pi/ll ;  topi [4. 31]  :=HPi/8?  topi[4.32]:=HPi/8; 
topi [4.39] :=7Pi/5 ;  topi[4 . 40] :=7Pi/5 ;  topi[4.45]:=17Pi/i2; 
topi [4 .48] :=10Pi/7 ;  topi[4.49] :=10Pi/7;  topi[4.53]:=13Pi/9; 
topi [4.54]: =13Pi/9 ;  topi[4.56] :=16Pi/ll;topi[4.57]:=16Pi/ll? 
topi[4.71] :=3Pi/2;  topi[4.85] :=17Pi/ll;  topi[4.86]:*17Pi/ll; 
topi [4 .88] :*14Pi/9 ;  topi [4. 89] :=14Pi/9;  topi [4. 93] :=llPi/7; 
topi [4.94]: =llPi/7 ;  topi [4. 97] :=19Pi/12;  topi [5. 02] :=8Pi/5? 
topi[5.03] :=8Pi/5;  topi[5.10] :=13Pi/8;  topi[5.11] :=13Pi/8; 
topi[5.14] :=18Pi/ll;  topi[5.23] :=5Pi/3;  topi[5.24]:=5Pi/3; 
topi [5.34]: =17Pi/10 ?  topi [5. 38] :=12Pi/7;  topi [5. 39] :=12Pi/7 ; 
topi [5. 42] :=19Pi/ll;  topi[5.43] :»19Pi/ll;  topi[5.49] :=7Pi/4; 
topi [5.50] : =7Pi/4 ;  topi [5. 58] :=16Pi/9 ;  topi[5.59]:=16Pi/9; 
topi [5.65] : =9Pi/5 ;  topi[5.71] :=20Pi/ll;  topi[5.75]:=HPi/6; 
topi [5. 76] :=HPi/6;  topi [5. 83] :=13Pi/7;  topi[5.89]:=15Pi/8; 
topi [5.93] : =17Pi/9 ;  topi[5.96] :=19Pi/10;topi[5.97] :=19Pi/10; 
topi [5. 99] :=21Pi/ll;topi[6. 00] :=21Pi/ll;topi[6.02] :=23Pi/12? 
topi [6. 28] :=2Pi;  topi [6. 54] :=25Pi/12;  topi [6. 56] :«23Pi/ll; 
topi [6.57] : =23Pi/ll ; topi [6.59] :=21Pi/10;topi[6. 60] :=21Pi/10; 
topi [6. 63] :=19Pi/9;  topi [6. 67] :=17Pi/8 ;  topi[6.68]:=17Pi/8; 
topi [6.73] : =15Pi/7 ;  topi[6.80] :=13Pi/6;  topi[6.81] :=13Pi/6; 
topi[6.85] :=24Pi/ll;  topi[6.91] :=HPi/5;  topi[6.98]:=20Pi/9; 
topi [7.06] : =9Pi/ 4 ;  topi[7.07] :=9Pi/4;  topi[7. 13] :=25Pi/ll; 
topi[7.14] :=25Pi/ll;  topi[7.18] :=16Pi/7;topi[7.22] :=23Pi/10; 
topi [7 . 23 ] : =2 3 Pi/ 10 ;  topi[7.33] :=7Pi/3;  topi[7.42]:=26Pi/ll; 
topi [7. 43] :=26Pi/ll;  topi [7. 46] :=19Pi/8;  topi[7.53]:=12Pi/5; 
topi [7.54]: =12Pi/5 ;  topi [7. 59] :=29Pi/12;  topi [7. 62] :=17Pi/7; 
topi [7 .63] :=17Pi/7 ;  topi[7.67] :=22Pi/9;  topi [7 . 68 ] :=22Pi/9 ; 
topi[7.7l] :=27Pi/ll;  topi[7.85] :=5Pi/2y  topi[7.99] :=28Pi/ll; 
topi[8.00] :=28Pi/ll;  topi[8.02] :=23Pi/9;  topi[8.03]:=23Pi/9; 
topi [8.07] : =l8Pi/7 ;  topi [8. 08] :=18Pi/7;  topi[8 . 11] :=31Pi/12 ; 
topi[8.12] :=31Pi/12;  topi[8.16] :=13Pi/5;  topi[8.17] :=13Pi/5; 
topi [8.24] :=2 lPi/8 ;  topi[8.25] :=21Pi/8;  topi[8.28] :=29Pi/ll; 
topi [8.37]: =8Pi/3 ;  topi[8.38] :=8Pi/3 ;  topi[8.48]:=27Pi/lO; 
topi [8.52] : =l9Pi/7 ;  topi [8 . 53 ] :=19Pi/7 ;  topi[8.56]:=30Pi/ll; 
topi [8.57]  :=30Pi/ll ;  topi [8. 63]  :=HPi/4 ;  topi [8 . 64]  :=HPi/4 ; 
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topi[8.72] :«25Pi/9;  topi[8.73] 
topi[8.80] :*14Pi/5;  topi[8.85] 
topi [8. 97] :=20Pi/7;  topi [8. 98] 
topi [9. 07] :=26Pi/9;  topi [9. 08] 
topi [9. 13] :*32Pi/ll;topi[9.14] 
topi[9. 42] :=3Pi;  topi[x_]:=x; 


25Pi/9;  topi[8.79] :=14Pi/5; 
31Pi/ll;  topi[8.90] :=17Pi/6; 
20Pi/7 ;  topi [9.03]: *23Pi/8 ; 
26Pi/9 ;  topi[9.11] :=29Pi/10; 
3 2 Pi/11; topi [9. 16] :=35Pi/12 ; 
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VII.  FILE  POLAR1.MA 


(*  Copyright  Dennis  A.  Polaski,  May  4,1993  *) 
polar 2 : “Module [ { 

>, 

cltranscend ? 

funcl“Input[ "Enter  one  of  the  polar  curves  which  bounds  the 
region. 

r  *  "  ] ; 

funcl“funcl  /.  sqrt->Sqrt; 
funcl“funcl  /.  pi->Pi; 
funcl=funcl  /.  q->Q; 

func2=Input[ "Enter  the  other  polar  curve  which  bounds  the 
region. 

r  =  "  ] ; 

func2=func2  /.  sqrt->Sqrt ; 
func2=func2  /.  pi->Pi; 
func2=func2  /.  q->Q; 

limitsknown=Input["If  the  limits  are  known  type  y, 
otherwise  type  n . " ] ? 

If [limit sknown===y,  { 

lowerlimit=Input [ "Input  lower  limit . " ] ; 
lowerlimit=lowerlimit  /.  sqrt->Sqrt; 
lowerlimit*=lowerlimit  /.  pi->Pi; 
upperlimit=Input[" Input  upper  limit."]? 
upperlimit=upperlimit  /.  sqrt->Sqrt; 
upper limit*upper limit  /.  pi->Pi; 

If [N[upperlimit]<N[lowerlimit] , lowerlimit=N[lowerlimit-2Pi] ] 

{ lowerl imit=0 ;  upper 1 imit=2  N [ Pi ] } ] ; 
bothf unc=  { fund ,  f unc2 } ; 
ptranscend  ? 

If  [{ fund,  func2 }  =  !  =  ({ fund,  func2 }  /.  Sqrt[x _ ]->x), 

theroots=NRoots [ f unci A  2==f unc2  A  2 , Q ] , 
theroots=NRoots[funcl==func2,  Q] ; 

]? 


If [theroots=! “False,  { 

theroots=  Table [theroots[ [i, 2] ] , { i , Length [ theroots ] } ] ? 
theroots=DeleteCases[theroots,  Q_Complex  /;  Abs[Im[Q] ]>.l] ; 
theroots=ReplaceAll [theroots,  Q_  ->  Re[Q]  ]; 
rmtranscend ; 

Q=. ; 
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theroots*Table [ FindRoot [ funcl**func2 ,  (Q, theroots [ [m] ] )  ] , 

{ m ,  Length [ theroots ] }  ] ; 

theroots*  Q  /.  theroots; 

theroots*N [Union[theroots ,  {>33? 

theroots*Union[ theroots,  theroots  +  N[Table[2Pi, 

{ j ,  Length [ theroots ]}]]]; 

testroots*Union [ Select [ N [ Round [ theroots* 1 00 ] / 100 ] , 
(#>=N[-Pi]  &&  #<*N[(5/2)  Pi  ])&],{}] ; 
theroots=Select [ theroots, (#<=  N[upperlimit]  && 

#>=  Nflowerlimit] ) 4] ; 

theroots*Union[N[ (Round[ 100*theroots] )/100] , { ) ] ; 
theroots=Map [ topi ,  theroots ] ; 
testroots=Map[topi,  testroots] ; 

}]; 

theroots*Transpose [ Select [Tablet { theroots [ [i] ] , 
(N[Abs[funri-func2]  /.  Q->  theroots [ [i] ]  ])<.05), 

{ i , 1 , Length [ theroots ]}],(#[[ 2 ] ] ===True )*]][[ 1 )] ? 

I f [ theroots=== {}[[!]]/ theroots* { } ] ; 

max={0,0) ; 

For[rr=l,  rr<=2,  { rval=rr ?  ++rr; 

max[ [rval] ]=Max[DeleteCases[ Table [bothfunc[ [rval] ] 

/.  Q->values,  (values,  0,  2Pi,  .25}],  x_Complex]]  }]; 
max=Max [ max ] ; 

plotroutine; 

ptranscend ; 

For [hh=l ,  hh<=2,  (hval=hh;  ++hh; 

otherroots=NSolve ( { r==0 ,  r==bothf unc [ [ hval ] ] } ,  { r , Q } ] ; 
otherroots=N [ { r , Q }  / .  otherroot s ] ; 


If [otherroots=!=(r,Q)  &&  other roots  !={},( 
otherroots=DeleteCases [ otherroots , ( r_ , Q_Complex  / ; 
Abs[Im[Q] ] > . 1 >  3 ; 

otherroots=ReplaceAl 1 [ otherroots , { r_ , Q_ }  ->  { r , Re [ Q ] } ] ; 
rmtranscend ; 

Q=.  ; 


otherroots=Table [ FindRoot (bothf unc [ [hval] ]*=0, 

{ Q , otherroots [[m,2]]}],{m, Length [ otherroots ] } ] ; 
otherroots=(0,Re[Q] }  /.  otherroots; 
otherroots=N[ Union [otherroots, { } ] ] ; 

otherroots=Union [ otherroots ,  otherroots+  N[Table[ (0,2Pi) , 
( i , Length [ otherroots ]}]]]; 

otherroots=Select[ otherroots, (#[ [2] ]<=N[upperlimit]  && 

# [ [2] ]  >=  N[lowerlimit] ) &] ; 

otherroots=Union[N[ (Round [100*otherroots] )/100] , { } ] ; 

>]; 
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otherroots=Union[ReplaceAll[otherroots,  {0,6.28}->{0,0)],{}]; 
otherroots“Union[ReplaceAll [otherroots,  {Q, r}->{ {Q, r }}],{}] ; 
If[hval*=l,  funclroots=otharroots,  func2roots*otherroots] ; 
}]; 


allroots={ funclroots,  func2roots} ; 
endprograic=0 ; 
runoption*0 ; 
diff check; 

If [ (theroots===False  jj  testdif f===yes)  &&  runoption==0  && 
Length [ funclroots ]<=1  &&  Length[ func2roots]<=l, { 

If [limitsknown===n,  { 

I f [ Length [ ther oot s ] == 1 , 

Print [ "  M , 

"The  two  polar  curves  intersect  at  Q  =  ", 
topi[theroots[ [1] ]  ]  ]]; 

If [Length [ther oot s]  >1, 

Print ["The  two  polar  curves  intersect  at  Q:  "," 

",  Map [topi,  theroots]  ]]; 

Pause [ 8 ] ; 

newcurve=Input [ "The  polar  curves  are  plotted  from  Q=0 
to  Q=2Pi.  If  you  want  to  change  the  limits  type  y,  otherwise 
type  n . " ] ; 

If [newcurve*==y, { 

lowerlimit=Input[ "Input  new  lower  limit."]; 
lowarlimit=lowerlimit  /.  sqrt->Sgrt; 
lower limit=lower limit  /.  pi->Pi; 

While [lowerlimit  <-N[Pi],  { lowerlimit=Input [ "The  lower  limit 
can  not  be  less  than  -Pi,  input  lower  limit  again."]; 
lowerlimit=lowerlimit  /.  sqrt->Sqrt; 
lowerlimit=lowerlimit  /.  pi->Pi; 

)]; 

upperlimit=Input[ "Input  new  upper  limit."]; 
upper limit=upper limit  /.  sqrt->Sqrt; 
upperlimit=upperlimit  /.  pi->Pi; 

While [upperlimit>2  N[Pi],  {upperlimit=Input["The  upperlimit 
can  not  be  greater  than  2Pi,  input  upper  limit  again."]; 
upperlimit=upper limit  /.  sqrt->Sqrt; 
upper limit=upper limit  /.  pi->Pi; 

}]; 

If  [N [upperlimit ]<N[ lower l  i  t] ,  lowerlimit=N[lowerlimit-2Pi]  ] ; 
rmtranscend;  plotroutine, 

)3  ; 

}]; 

complexcheck;  negcheck; 

If [Length [complexf unci] <=2  &&  Length [complexfunc2]<=2  && 
endprogram==0 ,  findarea2]; 

If  [  (Length  [comp  lex  fund]  >2  !!  Length [ complexf unc2]  >2)  && 
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endprogram==0 ,  f indarea 1 ] ;  runopt ion=l 

)]; 


If [testdiff==*yes  &&  (Length [ fv nclroots]  >1  || 

Length [func2roots]>l)  &&  runopt ion==0,  ( 

I f [ 1 imitsknown===n ,  { 

I f [ Length [ theroots ] ==1 , 

Print [ "  " , 

"The  two  polar  curves  intersect  at  Q  =  ", 
topi [theroots [  [  1  ]  ]  ]  ] 3 ; 

If [Length [theroots]  >1, 

Print["The  two  polar  curves  intersect  at  Q: 

",  Map [topi,  theroots]  ]]; 

Pause [ 8 ] ; 

newcurve=Input["The  polar  curves  are  plotted  from  Q=0 
to  Q=2Pi.  If  you  want  to  change  the  limits  type  y,  otherwise 
type  n . " ] ; 

I f [ newcurve===y , { 

lowerlimit=Input[w Input  new  lower  limit."]; 
lowerlimit=lowerlimit  /.  sqrt->Sqrt; 
lower limit=lower limit  /.  pi->Pi; 

While [lowerlimit  <-N[Pi],  { lower limit=Input[ "The  lower  limit 
can  not  be  less  than  -Pi,  input  lower  limit  again."]; 
lowerlimit=lowerlimit  /.  sqrt->Sqrt; 
lowerlimit=lowerlimit  /.  pi->Pi; 

}]; 

upperlimit=Input[" Input  new  upper  limit."]; 
upper limit=upper limit  /.  sqrt->Sqrt; 
upperlimit=upperlimit  /.  pi->Pi; 

While [ upper limit>2  N[Pi],  [upperlimit=Input["The  upperlimit 
can  not  be  greater  than  2Pi,  input  upper  limit  again."]; 
upperlimit=upperlimit  /.  sqrt->Sqrt; 
upper limit=upper limit  /.  pi->Pi; 

>]; 

If  [N  [upperlimit  ]<N[  lowerlimit] ,  lower  limit=N[  lower  limit-2  Pi]  ] ; 
rmtranscend;  plotroutine; 

)]; 

)]; 

complexcheck;  negcheck; 

If [ Length [complexfuncl]<=2  &&  Length [ complexfunc2 ] <=2  && 
endprogram==0 ,  findarea2]; 

If [ (Length [complexfuncl]  >2  |{  Length [complexfunc2 ]  >2)  && 
endprogram==0 ,  findareal];  runoption=l 
}]  ; 

If [ (theroots===False  Jj  Length [ theroots ]  <=1)  && 
runoption==0  && 

(  Length [funclroots]  >1  Length [ func2roots ] >1 ) ,  { 

If [limitsknown===n,  { 
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I f [ Length [ ther oots ] *=1 , 

Print  ["  "," 

The  two  polar  curves  intersect  at  Q  *  ", 
topi [ theroots [ [ 1 ] ]  ]  ] ] ; 

Pause [  8  ] ; 

newcurve=Input["The  polar  curves  are  plotted  from  Q=0 
to  Q*2Pi.  If  you  want  to  change  the  limits  type  y,  otherwise 
type  n. " ] ; 

I f [ newcurve===y , { 

lowerlimit=Input[ "Input  new  lower  limit."]; 
lowerlimit=lowerlimit  /.  sqrt->Sqrt; 
lower limit=lower limit  /.  pi->Pi; 

While [lowerlimit  <-N[Pi],  { lower limit=Input[ "The  lower  limit 
can  not  be  less  than  -Pi,  input  lower  limit  again."]; 
lowerlimit=lowerlimit  /.  sqrt->Sqrt; 
lowerlimit=lower limit  /.  pi->Pi; 

>]; 

upperlimit=Input[" Input  new  upper  limit."]; 
upper limit=upper limit  /.  sqrt->Sqrt; 
upperlimit=upperlimit  /.  pi->Pi; 

While [upperlimit>2  N[Pi],  {upperlimit=Input["The  upperlimit 
can  not  be  greater  than  2Pi,  input  upper  limit  again."]; 
upper limit=upper limit  /.  sqrt->Sqrt; 
upperlimit=upperlimit  /.  pi->Pi; 

}]; 

If  [N  [upperlimit  ]<N[  lower  limit]  ,lowerlimit=N[lowerlimit-2Pi]  ] ; 
rmtranscend;  plotroutine; 

)]; 

}]? 

complexcheck;  negcheck; 

If [Length [complexf unci] <=2  &&  Length [complexfunc2]<=2  && 
endprogram==0 ,  f indarea2 ] ; 

If [ (Length [complexf unci]  >2  jj  Length [complexfunc2]  >2)  && 
endprogram==0 ,  f indareal ] ; runopt ion«l 
}]; 

control=0; 

If [Length [theroots ]>=2  &&  limitsknown===n  &&  runoption==0, { 
Print ["The  two  polar  curves  intersect  at  Q:  "," 

",  Map [topi,  theroots]  ]; 

Pause[8] ; 

newcurve= Input [ "The  polar  curves  are  plotted  from  Q=0 
to  Q=2Pi.  If  you  want  to  change  the  limits  type  y,  otherwise 
type  n."]; 

I f [ newcurve===y , ( 

lowerlimit=Input[" Input  new  lower  limit."]; 
lowerlimit=lower limit  /.  sqrt->Sqrt; 
lowerlimit=lower limit  /.  pi->Pi; 

While [lowerlimit  <-N[Pi],  { lower limit=Input[ "The  lower  limit 
can  not  be  less  than  -Pi,  input  lower  limit  again."]; 
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lowerlimit-lowerlimit  /.  sqrt->Sqrt; 
lower limit*lower limit  /.  pi->Pi; 

)]? 

upperl imit-Input[" Input  new  upper  limit."]; 
upperl imit-upperl imit  / .  sqrt->Sqrt ; 
upperl imit-upperl imit  /.  pi->Pi; 

While [ upperl imit>2  H[Pi]f  { upperl imit-Input[ "The  upperlimit 
can  not  be  greater  than  2Pi,  input  upper  limit  again."]; 
upperl imit-upperl imit  /.  sqrt->Sqrt; 
upperl imit=upper limit  /.  pi->Pi; 

>]; 

If  [N  [upperlimit  ]<N[  lower  limit] ,  lowerlimit=N[lowerlimit-2Pi]  ] ; 

>]; 

badlimits=Select[N[theroots] , 

(#>N[ lower limit]  &&  #<N [upperl imit ])&] ; 

While [ Length [badlimits] >0,  { 

If [control==l,  { 

Print ["For  Q=-Pi  to  Q=5/2  Pi  the  two  polar  curves  intersect 
at  Q:","",  Mapftopi,  testroots]  ];  Pause[8]}]; 
control»control+l ; 

lowerlimit=Input["A  point  of  intersection  can  not  lie 

between  the  two  limits.  Input  new  lower  limit."]; 

lowerlimit=lowerlimit  /.  sqrt->Sqrt; 

lowerlimit=lowerlimit  /.  pi->Pi; 

upperl imit=Input[" Input  new  upper  limit."]; 

upperl imi t=upperl imit  / .  sgrt->Sqrt ; 

upperl imit=upper limit  /.  pi->Pi; 

If  [ N [upperl imit ]<N[lowerl imit] ,  lowerlimit=N[lowerlimit-2Pi]  ] ; 
badlimits=Select[N[testroots] , 

(#>N[ lower limit]  &&  # <N[ upperl imit ])&] ; 

>]; 

rmtranscend;  plotroutine;  complexcheck;  negcheck; 

If [ (Length [complexfuncl] <=2  &&  Length [complexfunc2]<=2)  && 
endprogram==0 , f indarea2 ] ; 

If [ (Length [complexfuncl]  >2  jj  Length [complexfunc2]  >2)  && 
endprogram==0, f indareal] ;  runoption=l 
>]; 

onceprint=0 ; 

If [Length [theroots] >=1  &&  1 imitsknown===y  &&  runoption==0 , { 
badlimits=Select[N[theroots] , 

(#>N[lowerlimit]  &&  #<N [upperl imit] ) &] ; 

While[Length[badlimits]>0,  { 

If [onceprint==0, { 

Print ["The  two  polar  curves  intersect  at  Q;  "," 

",  Map [topi, theroots]  ];  onceprint=l } ] ; 

Pause [ 8 ] ; 

If [control==l,  { 

Print ["For  Q=-Pi  to  Q=5/2  Pi  the  two  polar  curves  intersect 
at  Q:","",  Map[topi,  testroots]  ];  Pause[8]}]; 
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control=control+i ? 

lowerlimit=Input["A  point  of  intersection  can  not  lie 
between  the  two  limits.  Input  new  lower  limit."]; 
lowerlimit=lower limit  /.  sqrt->Sqrt; 
lowerlimit=lowerlimit  /.  pi->Pi; 
upperlimit=Input[" Input  new  upper  limit."]; 
upperlimit=upper limit  /.  sqrt->Sqrt; 
upperlimit»upper limit  /.  pi->Pi; 

If  [N[upperlimit]<N[  lower  limit] ,  lowerlimit=N[lowerlimit-2Pi]  ] ; 
badlimits=*Select[N[  test  roots] , 

(#>N[lowerlimit]  &&  #<H[ upper limit] ) 6] ; 

If [Length [badlimits]==0,  { 
rmtranscend;  plotroutine) ] ; 

}]? 

complexcheck;  negcheck; 

If  [  (Length [complexf unci] <-2  &&  Length [complexfunc2] <=*2 )  && 
endprogram==0, findarea2] ; 

If [ (Length [complexf unci]  >2  J J  Length [ complexf unc2]  >2)  && 
endprogram==0, findareal] ;  runoption=l 
)]; 

]; 


complexcheck: =Module[ ( } , 

complexf uncl=DeleteCases [ DeleteCases [Table [ N [ fund  / . 
Q->values],  {values,  N[ lower limit] ,  N [ upper 1 imit ] ,  .1}], 
r_Integer ] , r_Real ] ; 

complexf unc2=DeleteCases [ DeleteCases [ Table [ N [ f unc2  / . 
Q->values],  { values, N[lowerl imit] ,N[upperl imit] ,  .1}], 
r_Integer],  r_Real]  ]; 

diff check: =Module[ { } , 

dif ference=DeleteCases [Table [funcl-func2  / .Q->values, 
{values,  lowerlimit+.Ol,  upperlimit-. 01, . 25} ] ,  x_Complex] ; 
If [Max[difference]<=0  &&  Min[dif ference]<0  J | 
Max[difference]>0  &&  Min[difference]>=0, 
testdiff=yes,  testdiff=no]  ]; 

findareal:=Module[ { } , 

For[tt=l,  tt<=2,  {tval=tt;  ++tt; 

If [allroots [ [tval] ]  =!=  {{Q,r}>,  { 
otherroot s=Transpose [  all roots [ [ tva 1 ] ]  ] ; 
otherroots=otherroots[ [2] ] ; 

otherroots=Sort[Union[ReplaceAll [Join[ {N[lowerlimit] } , 
otherroots,  {N [upperlimit] } ] ,  0.->0],{}]]» 
otherroots=Select [otherroots , (#>=N[ lower limit]  && 

#<=N[ upperlimit] ) &] ; 

} , {other root s=ReplaceAll[ {N[ lower limit] ,N[upperlimit] } , 
0.->0] }]  ; 

rmtranscend;  area=. ;  totalarea={ } ; 

For [w=i ,  w<=Length[ otherroots ]-l,  { wal=w;  ++w; 

Q=  (otherroots [  [wal]  ]  +  otherroots [  [wal  +l]])/2; 

Iff  N[bothfunc[ [tval] ]  ]>0,  {Q=. ; 
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area*  NIntegratef .5*  bothfunc[ [tval] ] A2,  (Q, 

topi [otherroots[  [wal]  ]  ],  topi[otherroots[  [wal+1]  ]  ]}]; 

AppendTo [ totalarea ,  Abs [ area ]  ] } ] ; 

)J; 

area*  Apply [Plus,  totalarea]; 

If[tval*=l,  areal=area]; 

If[tval*=2,  area2=area]; 

}]; 

area*  area2  -  areal;  cltranscend;Q*. ; 
printfuncl  ] ; 

f indarea2 :*Module[ { } , 

Q=. ;  rmtranscend;  area=NIntegrate[ . 5* (func2A2  -  funclA2) , 
{Q/  topi[lowerlimit] ,  topi[ upper 1 imit] }] ;  printfunc2]; 

printfuncl; “Module [ { } , 

lf[area<=0,  (temp*funcl;  funcl=func2;  func2=temp) ] ; 

Print ["  "]; 

Print [ "  "  ] ; 

Print ["The  outer  polar  curve  is  r  =  ",  func2]; 

Print  ["The  inner  polar  curve  is  r  *  ",  fund]; 

Print [ "The  area  of  the  region  bounded  by  the  two"]; 

Print ["polar  curves  is  ", Abs [area]]  ]; 

printf unc2 ; “Module [ { } , 
cl transcend; 

lf[area<=0,  (temp=funcl;  funcl*func2;  func2*temp} ] ; 

Print ["  "]; 

Print [ "  " ] ; 

Print ["The  outer  polar  curve  is  r  *  ",  func2]; 

Print  ["The  inner  polar  curve  is  r  *  ",  fund] ; 

Print ["The  limits  of  integration  are  a=", 
topi[N[Round[lowerlimit*100]/100] ] , 

"  and  b=",  topi[N[Round[upperlimit*100]/100] ]  ]; 

Print ["The  area  of  the  region  bounded  by  the  two"]; 

Print ["polar  curves  is  ",  Abs [area]]  ]; 

negcheck:*  Module [ { } , 

negfuncl=Min[DeleteCases [Table [fund  /.  Q->values, 

(values,  lowerlimit+. 01,  upperlimit-. 01,  .25}],  x_Complex] ] 
negfunc2=Min[DeleteCases [Table [func2  /.  Q->values, 

(values,  lowerlimit+. 01,  upperlimit-. 01,  .25}],  x_Complex] ] 
lf[negfuncl<0  Jj  negfunc2<0, (cl transcend; 
lf[negfuncl<0  &&  negfunc2<0, 

(Print["Both  functions  have  negative  values  for  r."]; 
Print ["This  program  can  not  solve  such  problems."]}, 

If  [negfuncKO, 

(Print ["The  function  r  =  ",funcl,"  has  negative  values."] 
Print ["This  program  can  not  solve  such  problems."]  }, 
(Print ["The  function  r  =  ",func2,"  has  negative  values."] 
Print ["This  program  can  not  solve  such  problems."]}] 

]  ; 
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endprogram*l 

)] 

]; 


plotroutine:* 

Module [ { } , 

plotchunk=PolarPlot[  { fund,  func2),  {Q,  N[ lower limit] , 

N [ upper 1 imit] } ,  Display Function->Identity ,  AxesLabel-> 

{ "x" , "y" } ,  AspectRatio->Automat ic ,  PlotRange-> 

{{-1.02  max,  1.02  max), {-1.02  max,  1.02  max})]; 

rayl=Re[N[  (fund  /.  Q->  lowerlimit)  ]  ] ; 
ray2=Re[N[  (fund  /.  Q->  upper  limit)  ]  ] ; 
ray3=Re[N[ (func2  /.  Q->  lowerlimit)]]; 
ray4=Re[N[ (func2  /.  Q->  upper 1 imit) ]] ; 
llequivalent=N[Mod[lowerlimit,  2Pi] ] ; 
ulequivalent=N[Mod[upperlimit,  2Pi] ] ; 

If [lowerlimit  <  0,  llequivalent=  2  N[Pi]+lowerlimit, 
llequivalent=lowerlimit] ; 

If [upper limit  <  0,  ulequivalent=  2  N [ Pi ]+upper limit, 
ulequivalent=upperlimit] ; 

Q=.  ; 

If [ (rayl=-ray2  &&  llequivalent®*ulequivalent  && 
ray3==ray4)  |  |  (rayl==0  &&  ray2==0  &&  ray3==0  &&  ray4~0) , 
{Show[plotchunk,  DisplayFunction->  $Display Function] } , 
{Show[plotchunk,  Graphics[{ 

Line[ { { Cos [ lowerlimit] *rayl,  Sin [lowerl imit ] *rayl) , 

{Cos[lowerlimit]*ray3,  Sin [lowerl imit ] *ray3 } }] , 

Line[ { { Cos [ upper 1 imit ]*ray 2,  Sin [upper limit] *ray2 ) , 

{ Cos [ upperl imit ]*ray4,  Sin[upperlimit]*ray4} } ] , 

} ] ,  DisplayFunction->  $DisplayFunction] 

}]] 

ptranscend:=Module[ { } , 

sin[x_] ;=  x  -  XA3/6  +  xA5/120  -  xA7/5040  +  xA9/362880  - 

XA 11/399 16800+XA 13/62 2 702 08 00-XA 15/13 07674 3 68000+ 
XA 17/355687428096000-xA 19/12 1645100408832000; 
cos[x_] ;=  1  -  xA2/2  +  XA4/24  -  xA6/720  +  xA8/40320  - 
XA 10/3628800  +  XA 12/47 9001600 ; 
tan[X_] :=  x  +  XA3/3  +  2xA5/15  +  17xA7/315  +  62xA9/2835  + 
1382XA1 1/15592 5; 

cot[x_]:=  1/x  -x/3  -xA3/45  -  2xA5/945  -  xA7/4725  - 
2xA9/93555 ; 

csc[x_] ;=  1/x  +  X/6  +  7xA 3/360  +  31xA5/15120  + 

127xA7/604800  +  73XA9/3421440 ; 
sec[X_]:=  1  +  xA2/2  +  5xA4/24  +  61XA6/720  +  277xA8/8064  + 
5052 lx A 10/3628800; 
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exp[x_]:«  1  +  X  +XA2/2  +  XA3/6  +  XA4/24  +  XA5/120  + 
XA7/5040  +  XA8/40320  +  XA9/362880] 


XA6/720  + 
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VIII.  FILB  POLPLOTl.MA 


(*  Copyright  Dennis  A.  Polaski,  May  4,1993  *) 

polplotl:* 

Module [ { 

func,  limitsknown,  theroots,  lowerlimit, 
upperlimit,  nevcurve,  curvesknown,  temp 

>, 

Q=  •  > 

cltranscend; 

func=Input[ "Enter  the  polar  curve  which  bounds  the  region. 


r  =  "]; 

func=func  / .  sqrt->Sqrt ; 
func=func  /.  pi->Pi; 
func=func  /.  q->Q; 

revfunc=Replace[func,  Sqrt[x_]  ->  x] ; 
limitsknown=Input["If  the  limits  are  known  type  y, 
otherwise  type  n."]; 

I f [ 1 imi tsknown===y , 

{lowerlimit=Input[ "Input  lower  limit."]; 
lowerlimit=lowerlimit  /.  sqrt->Sqrt; 
lowerlimit=lower limit  /.  pi->Pi? 

While [lowerlimit  <-N[Pi],  {lowerlimit=Input["The  lowerlimit 
can  not  be  less  than  -Pi,  input  lower  limit  again."]; 
lowerlimit=lowerlimit  /.  sqrt->Sqrt; 
lowerlimit=lowerlimit  /.  pi->Pi; 

>]; 

upperlimit=Input[" Input  upper  limit."]; 
upperl imi t=upper limit  /.  sqrt->Sqrt; 
upperl imi t=upper limit  /.  pi->Pi; 

While [ upperl imit>2  N[Pi],  { upperl imit=lnput[ "The  upperlimit 
can  not  be  greater  than  2Pi,  input  upper  limit  again."]; 
upperl imit=upperlimit  /.  sqrt->Sqrt; 
upperl imit=up^er limit  /.  pi->Pi; 

)] 

}, 

{ lower 1 imit=0. ;  upperlimit*  2Pi}]; 
rmtranscend ; 

max=Max[DeleteCases[ Table [func  /.  Q  ->  values, 

{values,  0,  2Pi,  . 25} ] ,x_Complex] ] ; 

rayl=N[(func  /.  Q->lowerlimit) ] ; 
ray2=N[(func  /.  Q->upperlimit) ] ; 


llequivalent=N[ Mod [lower limit,  2Pi]]; 
ulequivalent»N[  Mod  [upper  limit,  2Pi]]  ,* 

If [lower limit  <  0,  llequivalent=  2  N [ Pi ]+lower limit, 

11 equivalent*  lowerlimit] ; 

If [upperlimit  <  0,  ulequivalent*  2  N [ Pi ]+upper limit, 
ulequivalent-  upperlimit] ; 

Q=  •  * 

If [ (rayl==ray2)  &&  (llequivalent—ulequivalent)  J  | 

(rayl==0  &&  ray2=*0) , { 

PolarPlot [ f unc , { Q ,  lowerl imit ,  upper 1 imit } , 

AxesLabel  ->  { "x" , "y" } ,  AspectRatio  ->  Automatic, 

PlotRange->  {{-1.02  max,  1.02  max), {-1.02  max,  1.02  max})]}, 
{ 

Show [ PolarPlot [f unc, {Q,  lowerlimit,  upperlimit}, 
DisplayFunction  ->  Identity,  AxesLabel  ->  {"x","y"), 
AspectRatio  ->  Automatic,  PlotRange->  {{-1.02  max, 1.02  max}, 
(-1.02  max,  1.02  max}}], 

Graphics [ { 

Line[ { {0,0} , {Cos [lowerl imit ]*rayl,  Sin[ lowerlimit] *rayl) } ] , 
Line[ { {0,0}, { Cos [ upper 1 imit ] *ray2,  Sin [upperlimit] *ray2 } } ] , 

)], 

DisplayFunction->  $OisplayFunction] } ] ; 
ptranscend ; 

theroots=NSol ve  [  { r==0 ,  r~revf  unc } ,  { r ,  Q }  ] ; 
theroots=N [ { r , Q }  / .  theroots ] ; 

If [theroots=l={r,Q}  &&  theroots  !={},{ 
theroots=DeleteCases [ theroots , { r_, Q_Complex  / ; 

Abs[Im[Q] ]>. 1} ] ; 

theroots=ReplaceAll[ theroots, {r_,Q_}  ->  {r,Re[Q] } ] ; 
rmtranscend ; 

Q=  • » 

theroots=Table [ FindRoot [ revf unc==0 , 

{Q, theroots [ [m,2] ] } ] , {m, Length [theroots] } ] ; 
theroots= { 0 , Q }  / .  theroots ; 
theroots=N[ Union [theroots, { } ] ] ; 

theroots=Union [ theroots ,  theroots+  N[Table[ {0,2Pi} , 

{ i , Length [ theroots ]}]]]; 

theroots=Select [theroots, (# [ [2] ]<=N[ upperlimit]  && 

#  [  [ 2 ]  ]  >=  N[ lowerlimit] )&] ; 

theroots=Union[N[ (Round [ 100*theroots] )/100] , { }] ; 

)]; 

rmtranscend ; 

If [limitsknown=!=y,  { 

If [theroots=!={r,Q}  &&  theroots!={}  &a 
Length  [  theroots  ]  >=*1,  { 

Print ["  "]? 


152 


Print ["The  polar  curve  r  is  equal  to  zero  when  Q  is: 

", Map [ topi, Transpose [theroots J [ [2] ]  ]  ] 

>]; 

Pause [ 4 ] ; 

newcurve=Input[“The  polar  curve  is  plotted  from  Q=0 

to  Q*2  Pi.  If  you  want  to  change  the  limits  type  y,  otherwise 

type  n. "3 ; 

I f [ newcurve===y , { 

lower limit* Input [“Input  new  lower  limit. “3* 
lowerlimit=lowerlimit  /.  sqrt->Sqrt; 
lower limit=lower limit  /.  pi->Pi; 

While [lowerlimit  <-N[Pi3,  (lowerlimit*Input[wThe  lower  limit 
can  not  be  less  than  -Pi,  input  lower  limit  again. "3; 
lowerlimit=lowerlimit  /.  sqrt->Sqrt; 
lower limit=lower limit  /.  pi->Pi; 

)]; 

upperlimit=Input[" Input  new  upper  limit."]; 
upperlimit=upperlimit  /.  sqrt->Sqrt; 
upper limit=upper limit  /.  pi->Pi; 

While [upperlimit>2  N[Pi],  {upperlimit=Input["The  upperlimit 
can  not  be  greater  than  2Pi,  input  upper  limit  again. "3; 
upperlimit=upper limit  /.  sqrt->Sqrt; 
upperlimit=upperlimit  /.  pi->Pi; 

}]  ; 

rayl*N[(func  /.  Q->lowerlimit) ] ; 
ray2=N[ (func  /.  Q->upperlimit) 3 ? 

llequivalent*N[Mod[lowerlimit,  2Pi} ] ? 
ulequivalent=N[Mod[upperlimit,  2Pi3  3  * 

If [lowerlimit  <  0,  llequivalent*  2  N[Pi]+lowerlimit, 
llequivalent*  lowerlimit] ; 

If [upperlimit  <  0,  ulequivalent*  2  N [ Pi ]+upper limit, 
ulequivalent*  upperlimit] ; 

Q=.  ; 

If [ (rayl==ray2)  &&  (llequivalent==ulequivalent) \  J 
(rayl==0  &&  ray2*=0) , { 

PolarPlot [ func , { Q ,  lowerlimit,  upperlimit), 

AxesLabel  ->  { "x” , "y" ) ,  AspectRatio  ->  Automatic, 

PlotRange->  {{-1.02  max,  1.02  max), {-1.02  max,  1.02  max})]}, 

{ 

Show [PolarPlot [func, {Q,  lowerlimit,  upperlimit), 
DisplayFunction  ->  Identity,  AxesLabel  ->  { "x" , "y" ) , 
AspectRatio  ->  Automatic,  PlotRange->  {{-1.02  max, 1.02  max), 
{-1.02  max,  1.02  max}}], 

Graphics [ { 

Line[ {{0,0}, {Cos [lowerlimit] *rayl,  Sin [ lowerlimit ] *rayl) } ] , 
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Line[ { {0,0}, {Cos[upperlimit]*ray2,  Sin [upper limit] *ray2 } } ] , 


)  3  > 

DisplayFunction->  $Display Function] }] 


>3? 

)]; 

3 
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IX.  FILE  POLFLOT2 . MX 


(*  Copyright  Dennis  A.  Polaski,  May  4,1993  *) 

polplot2 : “Module [ {  }, 

cl transcend; 

funcl=Input[ "Enter  one  of  the  polar  curves  which  bounds  the 
region. 

r  =  "  ] ; 

funcl=funcl  /.  sqrt->Sqrt; 
funcl=funcl  /.  pi->Pi; 
funcl=funcl  /.  q->Q; 

func2=Input[ "Enter  the  other  polar  curve  which  bounds  the 
region. 

r  =  "  ] ; 

func2=func2  /.  sqrt->Sqrt; 
func2=func2  /.  pi->Pi; 
func2=func2  /.  q->Q; 

limitsknown=Input [ "If  the  limits  are  known  type  y, 
otherwise  type  n."]; 

If [limitsknown===y,  { 

lowerlimit“Input[ "Input  lower  limit."]; 
lowerlimit==lowerlimit  /.  sqrt->Sqrt; 
lowerlimit=lowerlimit  /.  pi->Pi; 
upper 1 imit=Input [ " Input  upper  limit."]; 
upperlimit=upperlimit  /.  sqrt->Sqrt; 
upperlimit=upper limit  /.  pi->Pi), 

(lowerlimit=0;  upperlimit=2  N[Pi]}]; 
bothf unc=  { fund ,  func2 } ; 
ptranscend ; 

If  [{fund,  func2  }  =  !  =  ({  fund,  func2}  /.  Sqrt[x _ ]->x), 

theroots=NRoots[funclA2==func2A2 ,Q] , 
theroots=NRoots [ f uncl==func2 ,  Q] ; 

]; 


If [theroots=! “False,  { 

theroots=  Table [ theroots [ [ i , 2 ] ] , { i , Length [ theroots ] } ] ; 
theroots=DeleteCases[theroots,  Q_Complex  /;  Abs[Im[Q] ]>. 1] ; 
theroots=ReplaceAll [theroots,  Q_  ->  Re[Q]  ]; 
rmtranscend ; 

Q=  • ' 

theroots=Table[FindRoot[funcl==func2 ,  {Q, theroots [ [m] ] }  ], 

{ m ,  Length [ theroots ] } ] ; 
theroots=  Q  /.  theroots; 
theroots=N [Union [theroots,  { }] ] ; 

theroots=Un ion [theroots,  theroots  +  N [Table [2Pi, 
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{ j ,  Length [ theroots ] } ] ] ] ; 

testroots=Union [ Select [ N [ Round [ theroots* 100 ] / 100 ] , 

(#>=N[-Pi]  44  #<®N[ (5/2)  Pi] ) 4] , { }  ] ; 
theroots=Select [theroots, (#<=  N[ upper limit]  46 
#>=  N[ lower limit] ) 4] ; 

theroots=Union[N[ (Round [ 100*theroots ] )/100] , { } ] ; 

)]; 

theroots=Transpose [Select [Table [ { theroots [ [i] ] , 
(Abs[funcl-func2]  /.  Q->  theroots [ [i] ] ) <. 01 > , 

{ i , 1 , Length [ theroots ]}]/(#[[ 2 ] ] ===True) 4 ] ] [ i 1 ] ] ? 

I f [ theroots**® { } [ [ 1 ] ] , theroots® { } ] ; 

max® {0,0}; 

For[rr=l,  rr<=2,  { rval=rr ?  ++rr; 

max[ [rval] ]=Max[DeleteCases[ Table [bothfunc[ [rval] ] 

/.  Q->values,  (values,  0,  2 Pi,  .25}],  x_Complex] ]  }]; 

max=Max [ max ] ; 

plotroutine; 

I f [ 1 imitsknown===n ,  { 

I f [ Length [ theroots ] == 1 , 

Print [H 

The  two  polar  curves  intersect  at  Q  =  ", 
topi [theroots [ [1] ] ]  ]]; 

If [Length [theroots]  >1, 

Print ["The  two  polar  curves  intersect  at  Q:  "," 

" ,  Map [ topi , theroots ]  ] ] ; 

Pause [ 1 ] ; Pause [ 4 ] ; 

newcurve® Input [ "The  polar  curves  are  plotted  from  Q=0 
to  Q=2Pi.  If  you  want  to  change  the  limits  type  y,  otherwise 
type  n . " ] ; 

I f [ newcurve===y , { 

lowerlimit=Input[" Input  new  lower  limit."]; 
lowerlimit=lowerlimit  /.  sqrt->Sqrt; 
lowerlimit=lowerlimit  /.  pi->Pi; 

while [lowerlimit  <-N[Pi],  { lowerlimit=Input[ "The  lower  limit 
can  not  be  less  than  -Pi,  input  lower  limit  again."]; 
lowerlimit=lowerlimit  /.  sqrt->Sqrt; 
lowerlimit=lowerlimit  /.  pi->Pi; 

}]; 

upperlimit=Input[ "Input  new  upper  limit."]; 
upperlimit=upper'1  imit  /.  sqrt->Sqrt; 
upperlimit®upper limit  /.  pi->Pi; 

While [upperlimit>2  N[Pi],  (upperlimit=Input["The  upperlimit 
can  not  be  greater  than  2Pi,  input  upper  limit  again."]; 
upper limit=upper limit  /.  sqrt->Sqrt; 
upperlimit=upper limit  /.  pi->Pi; 

)]; 

rmtranscend;  plotroutine 
}]? 

)]; 

]? 
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